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'^This  report  pments  the  results  of  a study  of  the  effects  of  gear  pattern  on  airfield 
pavement  performance.  The  report  presents  the  methodology  required  to  express  air- 
craft traffic  in  terms  of  passes  rather  than  the  current  coverage  concept.  The  method 
b capable  of  considering  wheel  interaction  rather  than  only  surface  geometry, 
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THE  EFFECTS  OF  GEAR  PATTERN  ON 
SYSTEMS  PERFORMANCE 


1 INTRODUCTION 

Probitm.  A major  problem  in  airfield  pavement 
evaluation  and  design  is  assessing  the  damage  which 
will  he  laused  by  new  landing  gear  loads  and  com 
ple»  geometries  In  addition,  a need  exists  to  deter 
mine  the  design  life  and  life  remaining  in  an  existing 
airfield  pavement  This  problem  is  compounded  due 
to  the  effects  of  mixed  traffic  and  the  lack  of  a 
suitable  rationale  for  handling  mixed  traffic 

Background.  The  two  new  jumho  )et  aircraft  mm 
in  operation,  the  ( -SA  and  the  Boeing  74’’.  may 
require  a new  trend  in  development  of  future  aircraft 
pavement.  Sufficient  ground  flotation  or  load  distri- 
bution in  the  landing  gear  design  was  required  for 
the  C-5A  so  that  the  atreraft  could  operate  from  a 
medium  load  pavement  as  defined  in  TM  5-824-3.‘ 
The  Federal  Aviation  Administration  (FAA).  prior 
to  the  introduction  of  the  747,  placed  a requirement 
on  aircraft  designers  which  limited  the  pavement 
flexural  stress  to  that  which  was  developed  by  the 
stretched  venion  of  the  DC-8  aircref!.  If  this  trend  is 
lonlinucd  in  the  future,  pavement  requirements  for 
new  aircraft  may  level  off  at  some  bounding  limit.  If 
these  requirements  do  approach  a bounding  limit, 
the  pavement  engineer  will  no  longer  have  the  luxury 
ol  designing  for  the  heaviest  aircraft  and  ignoring 
the  lighter  iraffic  because  many  aircraft  types  will 
have  ntarb  the  same  pavement  requirements  As 

continue  to  grow  in  size  and  weight,  an 
economic  breaking  point  is  certain  to  be  reached 
where  it  will  be  more  economical  to  upgrade  the 
p.ivemeiits  than  to  provide  more  flotation  on  the  air 
craft  In  order  to  provide  the  necessary  flotation,  air 
craft  designers  will  probably  introduce  radically  dif 
fercni  geometries  and  varying  wheel  loads  all  of 
which  will  require  a more  precise  method  of  analysts 

Scope  anti  Objaetiva.  The  objective  of  this  study  is 
to  Identify  a technique  which  will  permit  OCF  to 
provide  an  improved  method  for  specifying  the  rela 
live  effects  of  traffic  loads  on  rigid  pavements  to 
reflect  variations  in  gear  patterns,  with  particular 
emphasis  on  multiwhecl  aircraft  including  the  C .SA 

Pavtmtnit  fnr  Airfklds  Othtr  than  Army  TM 
5«Z<-3(1970I 


and  747  aircraft.  The  scope  of  this  study  is  limited  to 
provi'Se  the  necessary  methodology  lo  satisfy  that 
objective 

Pravious  Attampts  to  Solva.  Other  investigators' 
have  approached  this  problem  in  different  ways 
Attempts  have  been  based  on  an  equivalence  m 
stress  developed  in  I he  pavement  fhe  relative 
seventy  of  two  landing  gear  systems  is  judged  by 
comparing  the  critical  flexural  stresses  generated  in 
the  pavemeni  slab  Fatigue  effects  are  accounted  for 
through  the  use  of  Miner's  hypothesis  which  assumes 
that  damage  is  relative  and  linearly  proportional  to 
Ihe  ratio  of  actual  irafTu  to  traffic  at  failure  for  each 
wheel  load  configuration  The  new  multiwhecl  land- 
ing gear  systems,  m addition  to  developing  critical 
stress,  also  dision  pavements  over  large  areas  at 
stresses  below  the  critical  stress  Prcsious  methods 
have  Ignored  the  suhcritical  stresses  and  ,rtvc  some 
limited  success  because  the  subcntical  stresses  were 
substantially  smaller  than  the  critical 

Th«or»tical  ConikkNVtions.  Pavement  response 
to  static  loading  can  be  predicted  with  reasonable 
accuracy  using  the  Westergaard  algorithm  Using 
Ihe  same  assumptions.  Westergaard.  Slel/cr  and 
Hudson’ have  developed  a finite  element  representa 
lion  of  a loaded  pa-  cmcnl  slab  which  allows  a more 
flexible  nuxleling  of  boundary  condiiions  Such  a 
linile  clemcnl  analysis  has  hern  adapted  lo  the  use  of 
complex  gear  patterns  as  a part  of  this  study 
(Appendix  A)  It  is  jHissible  lo  construe!  influence 
lines  for  a jminl  on  the  pavemeni  which  can  be  used 
to  approximate  the  stress,  strain,  and  or  deflection 
hiMory  lor  the  point  pnxfuced  by  a load  at  other 
jiomts  on  the  pavement  The  stress  history  at  a point 
in  a rii;id  pavement  which  results  from  the  move 
ment  ol  a multiple  wheel  landing  gear  system  is 
likely  lo  contain  a number  of  sir.  .s  excursions  of 
varying  magnitude 

Various  investigatois  have  had  success  in  desenb- 
ing  fatigue  failure  by  jierforming  a wcightexi  summa 
lion  of  Ihe  excursions  ol  stress,  strain,  or  dcfleition 
which  a specimen  cn>ounlers  Ihe  excursions  have 
lo  be  weighted  as  to  severity  for  obvious  reasons  for 

‘A  AfriAf'ri  for  h'tlimaltng  the  Lift  of  fttgtd  Paotmrnii 
OHDl.  Tschnwai  Report  Mo  4 2,1  119621 

‘C  K StvsLtw  »n<t  W K Hudiwti  A Ihrrrt  f'ompulrr 
'\»lutum  for  Ptalrr  atui  Paormtn!  Slabt  HcMioBrrh  Report 
•S6  9 l('«m«v  for  Ml^hway  Reseeirch  tjnivnr'iil)'  of  loxsjt 
1867) 


I'xample  a specimen  mav  be  capable  of  suftaintng  a 
'.irair.  t ji..-  • ol  ’<<'  iukto  m.  hc\  per  in  only  10 

MPii-  ‘nil  I van  siisiain  KKK)  entuMons  o(  micro 

pc  111  l!  a strain  ssile  relationship  exprcs'iid  as 
.U)  micro  inches  per  iiuh  x 10  cycles  2000  strain 
cveles  were  to  be  used,  as  a basis  of  performance,  the 
pasemrnt  could  sustain  fiO  micro  inches  per  inch  ct 
'UiO  V)  40  cycles  instead  of  the  ,1000  cycles 
measured  Obviously  a weighting  technique  must  be 
applied  to  the  magnitude  of  the  excursion  This  is 
,u  cimplished  bs  assuming  that  damage  is  propor- 
'ii'n.il  In  some  function  of  the  number  of  operations 

1 i>f  the  purjMwe  of  this  study  a log  function  has  been 
seieciixl  It  sfuuikl  be  emphasi/ed  however  that  this 
rvliiiionship  requires  further  study  as  additional  data 
t'ecomr  available 

2 BASIS  OF  ANALYSIS 

A number  of  factors  mu«t  be  considered  in  an 
anabsis  of  lompicx  loading  history.  Generally 
speaking  only  two  response  parameters  are  available 
'r;im  existing  analytical  procedures  used  in  the 
evaluation  of  pavements  deflection  and  bending 
moment  The  single  most  important  factor  to  be 
4 onsidered  for  pavements  loaded  near  the  center  of 
the  slab  is  gear  configuration,  i e..  single  wheel,  twin 
tandem,  or  multi-wheelcd  landing  gears. 

The  two  pavement  response  parameters,  dcflcc 
luin  and  bending  moment,  provide  a basts  for  deter 
inming  the  amount  or  degree  of  distortion  intro- 
duced in  a pavement  by  a particular  gear  loading 
I rom  an  analysis  view  point,  bending  moment  is  of 
more  value  than  deflection  because  a pavement  can 
conceivably  be  subjected  to  Urge  deOcction  with  a 
small  bending  moment  and  vice  versa  Bending 
moment  can  be  converted  into  stresses  and  strain 
.issuming  linearly  clastic  component  behavior 
( urrenl  design  procedures  arc  based  on  flexural 
stress  considerations  which  arc  directly  obtained 
from  bending  motnenl  utilizing  elastic  plate  bending 
•heory 

I’erformancc  data  for  pavements  subjjectcd  to 
various  gear  configuralion  loadings  arc  rather  dsffi 
lull  14'  obtain  Several  controlled  traffic  tests  have 
tieen  pel hirnied  by  the  Corps  of  Engineers  m past 
o .fs  I hesc  lesis  are  4»f  the  test  track  type  where  air 
■ lull  loading  IS  simulate*!  bv  a loading  apparatus 
'latfiiking  .1  'cncj  i>i  lest  slabs  >n  which  several 


parameters  are  varied  such  as  slab  thickness,  sub 
grade  strength,  etc  fhese  tests  have  been  conducted 
at  low  speed  and  typically  the  loading  apparatus  is 
capable  of  simulating  one  mam  gear  of  the  aircraft 
Performance  of  the  test  items  is  assessed  in  accord 
ance  with  current  Corps  of  Engineers  performance 
cTiteria.  Three  degrees  of  failure  have  been  cstab 
fished  by  the  Corps  of  Hngiiicers  as  follows 

1 Intiial  h'ailutr — That  point  when  a crack 
develops  m a pavement  and  extends  through  the 
entire  thickness  of  the  pavement  as  a result  ul  traffic 
loading  The  crack  must  be  due  to  an  exicrnallv 
applied  loading  and  not  to  shrinkage  or  other  non 
traffic  effects 

2 Skattrrtd  Slab — That  point  veherr  traffic 
induced  cracking  has  si  bdtvidcd  the  pavement  slab 
into  SIX  visible  pieces  As  with  the  initial  failure 
critenon.  the  cracks  must  be  due  to  extcrnaliv 
applied  loading  and  must  extend  through  the  tull 
depth  of  the  pavement  'lab 

.!  Compieie  Failure — That  point  where  the  pave 
ment  slab  has  been  cracked  into  about  pieces, 
each  piece  having  an  area  of  15  to  20  square  feel  Ai 
this  point,  the  pavement  would  be  considered  non 
operational  due  to  excessive  roughness  and  the  like 
lihncxl  of  debris  formation 

Due  to  the  inherent  assumptions  asscKiated  with 
the  available  pavement  analysis  prcKedurcs.  prim.in 
emphasis  has  been  placed  on  the  initial  failure  ton 
dition  The  shattered  slab  and  complete  failure  con 
ditions  represent  situations  too  unwieldy  for  e lalyti 
cal  modeling  at  this  time 

Actual  m scrvK'c  pavements  arc  subjected  a 
mixture  of  traffic.  Air  Force  bases  which  arc  dedi 
calcd  to  bomber  wings  will  experience  not  onlv 
bomber  traffic  but  also  support  type  traffic  such  as 
tankers,  cargo  aircraft,  etc  Some  dala  are  .ivadahle 
on  traffic  operations  at  selected  Air  Force  bases 
! hesc  data  should  be  analyzed  to  determine  prob 
-ole  traffic  miaturcs  for  selected  Air  Force  bases 

3 WHEEL  LOADING  PATTERNS 

The  iixuprint  ptoduf  cd  by  a C4>mpk»  ianibiig 
gear  can  be  examined  bv  cithei  4leflei!u>iri  ci;nu>urs 
Of  vifam  contours  Ibe  vietatl  < i ih»'  44>nt4vur  jia!t4  ru 
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u’lll  depend  on  the  pavement  thickness,  subgrade 
modulus,  the  elastic  properties  of  the  concrete  pave- 
ment slab,  wheel  loads  and  spacing,  in  addition  to 
the  slab  geometry  in  relation  to  the  wheel  loads. 

It  IS  possible  therefore  to  define  the  cyclic  strain 
or  deflection  history  on  a pavement  by  using  the 
footprint  as  an  influence  diagram  in  which  the  wheel 
loading  paths  relative  to  a point  in  the  pavement  are 
traced  on  the  contour. 

The  basts  for  this  technique  is  outlined  as  follows 
A plot  of  either  pavement  deflection  or  strain  re- 
sponse in  the  form  of  contours  can  be  used  to  predict 
pavement  response  due  to  traffic  loading.  Since  the 
concrete  and  foundation  materials  are  assumed  to 
behave  elastically,  superposition  can  be  used  to 
determine  response  histones  for  each  of  the  traffic 
lanes  I'he  contour  plot  moves  along  the  pavement 
coincident  with  the  wheel  gear  system.  Therefore, 
the  deflection  or  strain  history  at  a point  in  the  pave- 
ment can  be  constnicted  from  the  contour  plot 
simply  by  tracing  the  location  of  the  point  in  the 
pavement  on  the  contour  as  the  wheel  gear  passes 
over. 

An  example  of  the  use  of  such  an  approach 
follows:  I'igure  1 is  a contour  plot  of  deflection  for  a 
single  wheel  loading  on  a btekboume  test  slab.  It 
represents  the  ocflection  hisiory  of  a test  pavement 
subjected  to  a 37,000  pound  wheel  loading  resulting 
from  traffic  by  an  A-3  Tournapull  with  a Model  NU 
scraper.  The  traffic  lines  are  designated  on  the  con- 
tour plot  and  represent  the  path  followed  by  each  of 
the  loading  wheels.  These  plots  were  used  to  develop 
response  histories.  The  deflection  history  of  the 
pavement  at  the  transverse  joint  shown  on  Figure  1 
produced  by  single  wheels  traversing  each  of  the 
traffic  lines  would  then  be  computed  as  follows: 

Lane  1 Front  Wheel  - .023  in. 

Rear  Wheel  - .033  in. 

Lane  2 Front  Wheel  - .007  in. 

Rear  Wheel  - .013  in. 

initial  failure  occurred  after  90  operations.  The 
area  of  the  histogram  representing  the  deflection 
history  is 

.033  In  90  -b  .023  In  90  + .C07  in  90  4-  .013  In  90  » .342 

The  procedure  desciihed  in  the  above  example 
was  applied  to  other  test  pavements  subjected  to 
single,  twin  tandem,  and  12-wheel  assembly  loading. 


Pertinent  physical  properties  from  test  items  traf- 
ficked in  the  Lockfvourne  and  Multiple  Wheel  Heavy 
Gear  Load  (MWHGL)  tests  were  used  in  a finite 
clement  computer  program  for  pavement  analysis 
developed  by  Austin  Research  F.ngineers.  Inc.  This 
program  computes  both  deflection  and  principle 
stress  or  momcni  at  c.uh  nodal  nt  (Appendix  A) 
Thus,  a single  run  of  the  prograoi  rovides  a rather 
complete  picture  of  the  response  of  an  entire  pave- 
nicnl  slab  Weigbicd  liisiogrnni  areas  were  calcu- 
lated tor  the  initial  crack  failure  condition.  The  sum- 
mation of  the  weighted  area  histogram  is  indicative 
of  the  distortion  of  the  slab  under  load.  Histogra.'^s 
were  prixliiced  lor  single,  twin  tandem,  and  12-whi*^l 
landing  gear  arrays  with  .lO.fXX)  lb  wheel  loads.  The 
average  values  lor  the  weighted  areas  were  318  for 
single  wheel.  .52b  lor  twin  tandem  and  1.43  for 
12-whccl  assemblies  Using  the  single  wheel  loading 
as  a datum,  this  would  indicate  that  twin  (andem 
assemblies  distort  the  slab  1 66  times  as  much  as 
single  wheel.  These  values  are  indicative  of  the  strain 
relief  and  reversals  occurring  between  loading 
wheels. 


4 PAVEMENT  LOADING  HISTOGRAMS 

Previous  test  track  data  were  analyred  in  view  of 
the  strain-deflection  history  approach.  The  current 
Corps  cf  Hngincers  design  method  is  based  on 
critical  flexural  stresses  generated  at  a jointed  edge. 
An  attempt  to  produce  strain-deflection  histories 
was  based  on  ernditions  which  exist  at  a jointed 
edge.  This  initial  attempt  was  performed  using  data 
collected  from  the  Lockbourne*  and  MWHGL*  lest 
tracks. 

The  analysis  of  data  from  these  tests  presented 
problems  due  to  jointed  edge  behcvior.  For  design 
purposes  a jointed  edge  is  assumed  to  transfer  2S 
percent  of  the  load  across  the  joint  over  the  design 
lifc.’Tfic  25  percent  load  transfer  factor  was  estab- 
lished based  on  performance  data  from  many  experi- 
mental and  in-service  pavements.  While  the  load 
transfer  value  is  probably  valid  as  an  average,  some 


'lAKlibourne  No  I Tru  Track,  Final  Report  (U  S Anny 
<Hik)  River  Divuiion  Ijiboratones.  19461 

Multiple  WkrtI  Heavy  (!<far  Load  Pavement  Test*. 
WFS.TR-571  t7,  Vol  MV  (Waterways  Fxperiment  Station. 
19711 

•R  1.  Hutchinson.  lia$is  for  Ri^id  Pavement  Design  for 
Military  Airfields.  Miscellaneous  Paper  No  5-7  (Waterways 
Experiment  Station,  1966! 
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Fi(«re2.  Instrumentalion  layout  for  12  wheel  trafnc.  rigid  p.»vement  test  vclson 


variations  undoubtedly  occur  during  the  course  of 
the  design  life.  These  variations  have  a rather 
decided  impact  or  a strain  or  deHection  history 
analysis. 

MWHGL  Test  Track  Data.  The  data  collected 
trom  test  tracks  were  analyzed  using  the  strain 
history  approach.  Because  detailed  strain  history 
data  for  rigid  pavement  performance  are  available 
from  the  MWHGL  test  reports,  extensive  use  was 
made  of  this  data.  Figure  2 shows  a layout  of  the 
MWHGL  test  track  and  locations  of  instrumenta- 
♦'  )n.  TralTic  was  applied  to  the  test  track  pavement 
in  five  lines  as  shown  in  Figure  3 in  order  to  simulate 
the  wander  of  typical  operational  C-5A  aircraft.  A 
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f igure  Ea.  instrumentation  layouts  (or  the  twin  tandem 
assembly,  rigid  pavement  teat  section. 
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Figore  3.  Traffic  patterns  for  the  C 6 A wheel  assembly, 
rigid  pavement  test  section. 


similar  series  of  traffic  tests  was  applievi  to  the  test 
•^cms  to  simulate  747  traffic.  The  traffic  fines  for 
these  loadings  are  shown  in  Figure  4. 

The  response  data  for  each  test  slab  were  collated 
in  strain  history  format.  In  order  to  utilize  this  data, 
a procedure  for  weighting  the  strain  excursions  was 
necessary.  Several  weighting  functions  were  exam- 
ined during  the  initial  phases  of  the  study.  The 
method  chosen  to  weight  the  magnitude  of  strain 
excursions  which  are  measured  from  maxima  to 
minima  was  based  on  Saint  Venent's  maximum 
strain  theory  of  failure.  The  maximum  strain  theory 
assumes  that  failure  occurs  simultaneously  as  the 
elastic  limit  of  the  material  is  reached.  This  failure 
theory  is  primarily  applicable  to  materials  which  fail 
by  brittle  fracture  such  as  concrete.  The  maximum 
strain  for  conventional  paving  quality  concrete  is  on 
the  order  of  200  micro  inches  per  inch.  Thus,  strain 
exemsions  were  divided  by  200  micro  inches  per  inch 
to  provide  a ratio  of  the  magnitude  of  each  excursion 
to  the  maximum  strain  which  can  be  tolerated  by  the 
concrete.  A somewhat  similar  technique  was 
employed  to  weight  the  magnitude  of  the  deflection 
excursions.  Some  limited  data  were  collected  to 
determine  an  approximate  value  of  elastic  deflection 
normally  associated  with  initial  failure  of  rigid  air- 
field pavements  over  a normal  design  life.  The  value 
selected  u typical  was  0.05  inches.  It  should  be 
noted  that  the  selection  of  a base  deflection  or  datum 
value  can  be  rather  arbitrary  as  the  value  is  merely 
used  to  establish  a relative  weighting  index.  All 
excursions  are  expressed  as  a ratio  of  the  datum 
deflection. 

The  relative  merits  of  using  a strain  or  deflection 
history  were  examined.  More  meaningful  results 
were  obtained  using  the  strain  history  approach  as 
compared  to  deflection  because  of  the  large  number 
of  factors  affecting  deflection.  This  would  appear 
reasonable  as  pavement  design  and  analysis  pro- 
cedures are  pr^icted  on  stress  rather  than  deflec- 
tion. Another  consideration  which  tends  to  favor  the 
strain  history  approach  is  the  fact  that  equal  deflec- 
tions can  be  vastly  different  in  damaging  effects  due 
to  difference  in  radius  of  curvature.  Deflections  are 
less  sensitive  to  interaction  between  wheels  and 
gears. 

Measurement  of  strain  in  the  MWHGL  pavement 
tests  was  accomplished  at  the  venter  and  edges  of 
Test  Items  I.  2 and  3 in  the  longitudinal  and  trans- 


verse directions  for  the  total  test  program.  An 
example  of  a typical  strain  history  for  lane  2 is  shown 
in  Figure  5.  Histograms  were  prepared  showing  the 
number  of  strain  excursions  per  pattern  versus  the 
magnitude  of  the  excursion.  The  interval  for  the 
histogram  is  S micro  inches  per  inch.  Histograms  for 
strain  measured  in  test  items  I to  3 are  shown  in 
Figures  6 through  8.  Each  strain  excursion  was 
weighted  by  dividing  by  200  micro  inches  per  inch 
and  the  number  of  occurrences  of  thb  weighted 
strain  was  considered  by  multiplying  the  quotient  by 
the  natural  logarithmic  of  the  number  of  occur- 
rences. This  calculation  was  performed  for  each 
interval  of  the  histogram  and  summed.  Table  1 
shows  the  value  of  each  of  the  weighted  histograms 
for  the  C-5A  test  program,  it  should  be  noted  that 
these  weighted  histogram  areas  must  be  evaluated  in 
light  of  the  performance  of  the  MWHGL  test  items. 
Item  1 had  a crack  prior  to  the  application  of  traffic, 
lest  item  3 experienced  a first  crack  prematurely 
due  to  a subgrade  pumping  failure.  Test  item  2 
performed  in  a normal  manner  and  can  be  used  as  a 
basis  for  normalization  of  the  data. 
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32.13 
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50.23 

55.76 
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69.84 
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2SCL 

56.16 

61.24 

2SCT 

14.27 

16.82 

3 

38WJT 

21.56 

48.99 

3SCL 

4.56 

10.33 

38CT 

13.96 

32.81 

The  use  of  jointed  edge  histories  was  later 
abandoned  due  to  difficuities  in  correlation  brought 
on  by  mathematically  modeling  the  jointed  edge. 
The  assignment  of  a value  for  load  transfer  at  a 
jointed  edge  is  at  best  a tenuous  estimate.  Joint  effi- 
ciency changes  with  time  and  increased  traffic 
volume  and  has  a large  impact  on  strain  and  deflec- 
tion history  computations.  For  design  purposes  a 
load  transfer  value  of  25  percent  is  assumed  for  the 
entire  design  life  of  a pavement.  The  25  percent 
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Pigwn  4.  TrsISs  for  twin-UiHkm  whMl  UMnMy,  rigid  (MVMnent  tMt  MctioR. 
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FREE  EDGE 


pav{ 


a«Miinp(lkm  has  been  demonstrated  to  be  reasonable 
for  design  puipoees  thnngh  test  track  csperiments 
and  inspecthms  of  in>serrke  pavements.  The  strain 
iiHtory  approach  b much  too  sensitive  to  the  load 
transfer  vaHm  to  use  an  average  vahte  for  the  life  of 
the  pavement.  Emphasis  was  shHfod  from  the  use  of 
a jointed  edge  loadHif  to  an  interior  loading  condi* 
tkm.  The  use  of  an  inUrior  loading  condition  etimi* 
nates  a n:'  .^ber  of  boundaiy  conditions  which  arc 
difficult  to  represent  mathematkalljr. 

The  twin  tandem  747  traffic  which  was  applied  to 
the  test  items  provides  a mechanism  for  evaluating 
the  interior  slab  loading  toehnique  in  as  much  as  thk 
traffic  was  applied  at  the  center  of  the  slabs  as  shown 
in  Figure  4.  Histograms  for  the  strain  eicuriions  in 
Test  Item  2 prodwead  bjr  this  traffic  are  shown  in 
Figure  9.  Table  2 shows  the  weighted  areas  for  these 
histograms  for  test  hems  2 and  3.  The  areas  were 
adjusted  to  the  loHial  failure  condhioa  and  IS  traffic 
patterns  for  test  Hems  2 and  3,  respectively. 
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S MIXED  TRAFFIC 

Chapter  3 discussed  a procedure  for  generating 
load  histories  from  gear  pattern  contours  for 
particular  aircraft  coafiguratioM  and  pavement 
designs.  Chapter  4 presented  a hypothesis  which 
establishes  a totat  weighted  area  traffic  histogram, 
Hj.  which  is  a unique  characteristic  of  the  fatigue 
lifr  of  each  pavement  system. 

In  order  to  esaasine  the  rsiative  severity  of  mind 
traffic  loading  fsrwhkh  test  track  is  net  available,  a 
series  ofcalculattohe  can  bn  made  using  the  program 
described  in  AnsMndfac  A to  caleulaSs  waiahlad  hlslo> 
ifsni  sfsas  The  nsnder  nattern  fnr  traffir  im  a 


way  is  as  shown  in  Figure  10.  Traffic  distribution 
which  is  typically  contained  in  the  central  40  inches 
of  the  tasiway  wander  pattern  accounts  for  75 
percent  of  all  traffic.' 

The  technique  for  generating  histograms  is  based 
on  the  calculated  stress  distribution.  The  weighted 
area  histogram  (Hg)  is  obtained  as  follows: 


where  Aej » sum  of  stress  excursions  in  lane  j 

Ef  « Young’s  modulus  for  the  concrete  pave- 
ment 

(f  * maximum  strain  at  failure 
Aj  * proportional  area  of  the  wander  pattern 
in  lane  j 

M the  number  of  lanes  in  the  distribution 
pattern. 


H5  represents  the  area  of  the  weighted  histogram  for 
a single  pass.  For  multiple  passes  represented  by  the 
wander  distribution  Aj  the  traffic  histogram  Hm  b 
obtained  from: 


|Eq2) 

If  we  wbh  to  compare  complex  traffic  patterns,  we 
must  solve  equatiM  2 for  Np. 


EctfHM-Z(AojlnAj) 

"p- 


lEflJ! 


Equation  3 can  be  used  to  estimate  the  impact  of  a 
compbx  mix  of  traffic  simply  by  summing  the  indi- 
vidual dbtribution  patterns  as  follows; 


EcifHT-II(AosjlnA|jNp 


(Eq4] 


where  Ht  » the  design  weighted  hbtogram  area  for 
the  life  of  the  pavement 
0)}  » the  stress  excursion  for  aberaft  i in 
lanej 

N|  "ratio of  i type  aircraft  operatioos  to  total 
aircraft  operations  occurring  on  a 
periodic  basb. 


’UebMt  F.  Bahw,  JUport  of  Photogrommotrk  Motkod* 
for  MootoHof  La*onU  MeewnMt  of  Aimnfi,  LockbwrM 
AFB  lOWe  State  Uaivmby,  IIM). 


19 


I § u 1 8 1 8 1 s 1*8  n 1 1 n I i 8 ° 


I- 

i- 

i-'. 


Fifwt  It.  Wander  pattern  for  runway  loading. 


TaUt3 

Wrighted  Strain  HlilegreM  Aimm  Hgi  far  I. ON  ftaaM 


Gear  CiaWgaratlaa 

TatalAraa 

26,0001  Single  Wbail 

3.52 

B-62 

57.49 

C-UI 

27.90 

B-47 

24.79 

KCI36 

31.27 

Table  3 shows  accumulated  stress  excursion  levels 
for  the  25,000  lb  single  wheel  load,  the  B-S2.  B-47. 
KC- 135  and  0141  ain  vafi  systems.  The  pavement  is 
a 12'inch  thick  slab  with  a r^ius  of  relative  stiffness 
of  SO  in.  Weighted  histograms  were  then  developed 
by  summing  the  individual  stress  excursions  aimg 
traffic  lanes  spaced  for  each  12-inch  segment  of  the 
wander  pattern  which  represents  the  appropriate 
trafRc  density. 

Table  4 is  a report  of  aircraft  traffic  compiled  by 
base  operations  of  a typical  heavy  load  Air  Force 
base  sdected  at  random.  Data  of  this  type  have  been 
compiled  for  some  24  heavy  load  Air  Force  bases 


over  a four-year  time  period.  On  the  average,  the 
typical  heavy  load  Air  Force  base  experienced  100 
cycles  of  B-52  traffic  per  month  and  95  cycles  of  KC 
135  traffic.  A considerable  amount  of  variation 
occurred  in  the  traffic  operations,  however,  as  the 
standard  deviation  for  B-52  traffic  was  64  and  for 
KC  135  the  standard  deviation  was  51. 

Some  of  the  variation  is  probably  due  to  differ- 
ences in  the  numbers  of  aircraft  stationed  at  a base. 
In  some  instances  different  numbers  of  squadrons, 
groups  or  wings  are  stationed  at  different  Air  Force 
bases.  Variations  in  numbers  of  aircraft  are  reflected 
in  the  above  standard  deviations. 

The  distribution  of  aircraft  traffic  on  taxiway  and 
runways  has  been  studied  to  a limited  extent*  and 
has  been  shown  to  be  approximately  normally  dis- 
tributed. These  studies  indicated  that  75  percent  of 
aircraft  traffic  is  concentrated  with  a wander  of  7.5 
feet  on  primary  taxiways  and  37.5  feet  on  runways  at 


'Robart  F.  Bakar,  Rtport  of  Phologrommttrie  MtlHodt 
for  UtoiuriHg  Lctml  Macamanr  of  Aircrofi,  Lockbouma 
AFB  (Ohio  SUta  University,  1»M). 


21 


about  the  runway  midpoint.  The  Corps  of  Engineers 
design  method  for  pavements  is  based  on  75  percent 
of  the  traflic  occurring  within  the  limits  of  the  traffic 
area  width.  The  data  generated  above  from  past 
operational  records  would  indicate  that  over  a 
20-year  design  iife  at  an  average  heavy  load  base. 
24.000  cycles  of  B-S2  and  22.800  cycles  of  KC  135 
traffic  will  be  encountered.  Applying  the  75  percent 
adjustment  for  traffic  design  yields — 18,000  cycles  of 
B-52  and  17,100  cycles  of  KC  i 35  traffic.  Using  the 
existing  Corps  of  Engineers  method  for  describing 
traffic,  the  above  traffic  volumes  would  be  converted 
to  9000  coverages  of  B-52  traffic  for  taxiways  and 
4500  coverages  for  runways.  These  figures  agree  well 
with  10,000  coverages  and  5000  coverages  used  for 
design.  In  the  present  design  method  for  heavy  load 
pavements  the  KC  135  traffic  would  be  ignored.  This 
is  in  contrast  to  the  relative  strain  history  areas 
shown  in  Table  3 which  indicate  that  the  KC  135 
traffic  is  54  percent  as  severe  as  the  B-52  traffic. 
Thus,  the  total  traffic  would  be  computed  as  9000  -f 
17.1/18x.S4  x 9000»  13,617  coverages  for  the  taxi- 
ways  and  implies  that  the  traffic  volume  used  in  the 
design  of  taxiways  and  runways  b understated  by  36 
percent. 


6 ANALYSIS  AND  DISCUSSION 

Several  problems  arose  during  the  course  of  this 
study  which  tend  to  limit  the  degree  of  validation 
which  is  possible.  No  data  are  available  on  the 
behavior  of  pavements  under  mixed  traffic  loadings 
under  ideal  test  conditions.  Data  from  various  test 
tracks  with  similar  properties  and  performance 
records  were  compared.  A prime  consideration  in 
studying  the  effects  of  wheel  interaction  U the  radius 
of  relative  stiffness  which  controb  wheel  interaction. 
Thus  past  test  tiack  data  had  to  be  analyzed  on 
approximately  the  same  basis  to  achieve  meaningful 
comparisons.  Thb  requirement  limited  the  data 
available  for  comparative  analysis. 

Computer  anaiysb  of  pavement  respon'.e  is  essen- 
tial when  using  the  strain  hbtory  approach.  No  con- 
venient manual  method  could  be  devised  which 
would  yield  meaningful  results.  The  output  required 
to  produce  a significant  strain  hbtory  Is  too  detailed 
for  manual  computation  except  fmr  single  wheel 
loadings.  Wheel  interactions  on  multi>«vheeled  gears 
are  too  complex  to  produce  without  the  aid  of  a 
digital  computer. 


It  has  been  necessary  to  make  a number  of 
implicit  assumptions  relative  to  the  hypothesis 
advanced  in  this  study  which  defines  a weighted  area 
histogram  as  a unique  characteristic  cf  each 
pavement  systems.  These  assumptions  are  listed 
below: 

1 . A lower  threshold  for  fatigue  failure  produced 
by  strain  excursions  in  a pavement  does  not  limit  the 
level  at  which  a significant  portion  of  the  strain 
excursions  are  accumulated  during  the  pavement 
life. 

2.  The  relative  severity  of  strain  excursions  is  a 
function  of  the  log  of  the  number  of  loading  cycles. 

3.  The  weighted  histogram  area,  Hj,  at  which 
fatigue  failure  occurs  is  a function  of  the  strain 
excursion  history  only  and  is  unique  for  each  pave- 
ment system  design. 

The  strain  hbtory  technique  yields  reasonable 
results  and  can  provide  a basis  for  designing  airfield 
pavements  for  a given  level  of  mixed  traffic.  The 
method  provides  a much  more  realistic  assessment  of 
traffic  on  rigid  pavements  than  the  coverage  concept 
as  wheel  interaction  b taken  into  consideration.  TIk 
coverage  concept  merely ' isiders  surface  geometry; 
interactive  effects  are  lefl  to  judgment  and 
experience. 


7 CONCLUSIONS  AND  RECOMMENDATIONS 

Conclusions.  Based  on  the  results  of  the  study 
reported  herein,  several  conclusions  can  be  drawn 
concerning  the  strain  history  method  of  assessing 
field  traffic  loadings. 

The  strain  hbtory  approach  appears  to  be  a 
viable  technique  to  assess  multi-wheeled  aircraft 
traffic  effects  on  pavements. 

The  strain  history  approach  should  be  based  on 
interior  loading  conditions  due  to  the  difficulties 
encountered  in  modeling  jointed  edge  boundary  con- 
ditions. If  sufficient  assumptions  are  made  concern- 
ing the  jointed  edge  boundary  conditions  to  permit  a 
comparative  analysis  from  the  same  reference,  the 
anaiysb  may  as  well  be  performed  using  interior 
loading  conditions.  The  analysis  is  intended  to 
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demonstrate  relative  efTects  rather  than  absolute 

values. 


the  stresS'Ktrain  diairam  with  a minor  influence  up 
to  that  point. 


Tile  deflection  hbtory  approach  was  abandoned 
early  in  the  study  as  the  deflexion  approach  was  not 
sensitive  enough  to  wheel  interaction  to  provide 
meaningful  data. 

The  strain  history  approach  is  a reasonable  tech- 
nique for  analyzing  the  effects  of  mized  traffic  on 
airfield  pavement  performance.  The  weighting  func- 
tions applied  to  strains  and  repetitions  wilt  permit 
mixed  traflic  to  be  assessed  fro.n  a common  basis. 
Insufficient  data  is  available  to  vc.ify  the  nature  of 
the  Fatigue-Cyclic  Excursions  (F-N)  Curve  for  con- 
crete pavements.  The  assumption  that  the  relation- 
ship is  a log  ^utKtion  needs  to  be  verified.  An 
assumption  that  fatigue  strain  is  linearly  propor- 
tional to  the  ratio  of  the  strain  excursion  experienced 
to  that  at  failure  is  undoubtedly  an  oversimplifica- 
tion. The  relationship  is  probably  influenced  by  the 
strain  level  beyond  that  purely  linearly  position  of 


RnooniiiMndstloiw- 

1 . The  Corps  of  Engineers  rigid  pavement  design 
criteria  should  be  based  on  hbtograms  of  strain  dis- 
tribution rather  than  coverage  concept  for  describing 
traffic  volume. 

2.  Lateral  positioning  of  aircraft  on  taxiways  and 
runways  during  normal  operations  should  be  investi- 
gated. Previous  studies  are  rather  old  and  should  be 
updated  to  reflect  changes  in  aircraft  and  opera- 
tional characteristics  which  may  affect  aircraft 
placement  during  normal  operations. 

3.  The  nature  of  the  weighting  factor  to  be 
applied  to  strain  excursions  should  be  verified  by 
either  model  tests  or  traffic  data. 

4.  A program  for  the  collection  and  analysis  of 
field  data  for  corrtlation  should  be  undertaken 
particularly  for  new  pavements. 
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APPENDIX  A:  SLAB  30E— A COMPUTER 
PROGRAM  FOR  AIRFIELD  SLAB  ANALYSIS 

I'hii  appendix  presents  the  dtxrumcnfation  of  a 
computer  program  developed  by  Austin  Research 
Engineers  Inc.  for  the  U.S-  Army  Construction 
Engineering  Research  Laboratory  (CERL).  Cham- 
paign. Illinois.  The  work  was  conducted  under  Con 
tract  No  DACA  23-70-C-007b  by  Mr.  John  J. 
Panak.  Mr.  Frank  L.  Endres  was  the  programmer 
who  made  the  revisions  to  the  original  program  and 
performed  the  checkout  work  on  the  CERL  input 
facility. 

The  FORTRAN  Program.  Program  SLAB  30E  is 
an  extensive  u-daie  and  revision  of  a program  called 
SLAB  30  which  is  reported  by  Stelzer  and  Hudson.* 
it  has  features  which  allow  it  to  work  more  efficiently 
than  SLAB  30,  and  in  addition,  includes  three 
significant  new  developments  which  arc: 

1.  Addition  of  a printer  plot  routine  to  allow  a 
visual  display  of  selected  lines  of  output  values  of 
deflection,  bending  moments  and  principal  stresses. 

2.  Addition  of  a feature  which  will  allow  a 
particular  wheel  load  pattern,  such  as  that  associ- 
ated with  the  C-5A.  to  be  simply  input  to  the  com- 
puter model  by  reference  to  a single  input  coordi- 
nate. 

3.  Addition  of  a data  input  generation  routine 
which  will  allow  much  simpler  definition  of  the  pave- 
nent  slab  in  terms  with  which  the  investigator  is 
accustomed. 

The  program  is  written  to  operate  on  any  Control 
Data  Corporation  Model  20&  remote  user  terminal 
which  can  accers  a CDC  6600  computer.  All  pro- 
gramming is  in  ASA  FORTRAN  and  therefore  is 
compatible  with  other  computers  such  as  IBM  360 
and  UNIVAC  1108  systems. 

Stomg*  RaquIrmiMnts.The  storage  requirements 
arc  variable,  depending  upon  the  size  of  the  problem 
to  be  run.  Cards  which  must  be  changed  for  different 
sized  problems  are  specified  at  the  beginning  of  the 
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program  and  only  include  the  dimension  statements 
and  two  variables  which  define  the  number  of  mere 
ments  in  both  directions  It  is  not  necessary  to  malvh 
exactly  the  dimensioned  storage  to  the  problem,  any 
si/c  larger  is  also  acceptable  It  is  recommended  that 
several  dimension  packages  be  made  of  those  cards 
with  RE-DIMEN  in  columns  73  through  80  of  the 
main  driving  program  (refer  to  program  listing). 
These  packages  can  be  of  multiples  ol  width  and 
length  ot  10  or  20  stations  up  to  the  maximum 
acceptable  for  the  particular  computer  being  used,  A 
ploi  of  the  CDC  6600  storage  requireine^'  is  shown 
in  Figure  Al.  If  necessary  to  gain  more  storage 
space,  certain  variables  could  be  placed  in  common. 
In  addition,  others  may  be  set  equal  to  each  other  in 
storage  by  an  equivalence  declaration.  (.  ommon  and 
equivalence  statements  have  not  been  added  in  the 
program  for  normal  operation,  to  avoid  initial  con 
fusion  when  converting  to  other  computer  systems. 
All  subroutines  in  the  program  arc  variably  dimen 
sioned  as  functi-ins  c the  short  (X)  and  long  (Y) 
lengths,  which  arc  specified  in  the  mam  driving 
program. 

Input  ol  Probl»ms.  The  procedures  that  are 
followed  in  input  ol  problems  are  outlined  in  the 
Guide  for  Data  Input.  A parallel  study  of  the  guide 
will  help  to  understand  the  lollowing  discussion.  Any 
number  of  problems  may  be  run  at  the  same  time. 

The  first  two  cards  of  a problem  series  arc  for 
identification  purposes.  Any  alphanumenc  desenp 
tivc  information  desired  can  be  entered  it  is  sug 
gested  that  the  date  of  the  run  and  the  user’s  name 
always  be  entered  within  these  two  cards.  The  next 
card  is  the  problem  number  card  wuh  a brief 
description  of  the  particular  problem  The  problem 
number  itsell  may  contain  alphabetical  characters  if 
desired.  The  problem  serifs  terminate'  when  a blank 
problem  number  is  enci  untered 

TabiM  of  Data  Input.  Table  3I  is  used  to  inpuf  s,'*; 
problem  control  data  and  is  comprised  of  -i  agie 
datacard  that  specilks  the  multiple  Sv.,jd  rtOd 

the  number  of  caros  in  th?  following  data  fabks  8 ‘ 
through  B7.  The  Multiple  Load  Optico  m c luinn 
ol  Table  Bl  is  left  blank  if  each  sutv*' .si.e  p'-.-b-sn*.  is 
independent  of  the  preceding  p'obten.  if  a toiiowiog 
problem  is  lor  the  same  s:^!'  out  i.uly  the  kiad 
pattern  and  placement  a'e  to  change  the  iirst  prob 
lem  in  the  loading  senes  is  sr » cified  wuh  a + i for  the 
option  This  will  he  itc*  ' orfrtot”  problem  Each 
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successive  loading  must  have  a - 1 for  the  option  and 
these  will  be  termed  “offspring"  problems.  When  a 
blank  option  or  another  +1  is  encountered,  that 
problem  is  then  another  independent  problem  or  a 
new  parent  problem. 

The  optional  stiffness  input  switches  in  columns 
55.  60.  and  65  may  be  used  to  automatically  com- 
pute conventional  slab  or  plate  stiffness  values  by 
entering  a one  for  the  options  selected.  Appropriate 
values  must  then  be  made  available  in  Table  B2.  If 
D*  or  oy  or  C are  to  be  computed,  then  there  must 
be  real  values  available  for  E and  t in  Table  B2.  if  a 
uniform  value  of  S is  to  be  computed,  then  a sub- 
grade  modulus  k must  be  available.  The  stiffness 
options  must  be  left  blank  for  offspring  problems. 

Two  Output  Options  may  be  specified  in  columns 
75  and  80.  All  of  the  detailed  output  of  Table  B8 
may  be  omitted  if  desired.  The  user  is  cautioned 
when  using  this  option  that  he  could  omit  what 
might  later  be  desirable  information.  Deletion  of 
output  reduces  central  computer  time  only  a small 
amount.  It  may,  however,  decrease  turn-around 
time  due  to  reductions  in  peripheral  processor  time 
and  printer  time.  If  the  output  is  suppressed,  a 
single  referr'ce  station  is  printed  which  is  at  or  nep,r 
the  center  of  the  slab.  The  other  output  option  is 
exercised  if  the  user  desires  the  printout  of  principal 
stresses  in  place  of  principal  moments.  If  the  stress 
option  is  exercised,  an  appropriate  value  of  stab 
thickness  must  then  be  available  in  Table  B2.  The 
principal  moment  is  converted  to  a stress  having  the 
same  sign  by  multiplication  of  the  moment  by  the 
plate  section  modulus  which  is  internally  computed 
on  the  assumption  of  a uniform  thickness  of  plate. 
The  use  of  the  stress  option  is  only  appropriate  for 
plates  of  a constant  thickness.  At  specify  discon- 
tinuities in  the  slab,  such  as  a crack  or  joint  which 
might  be  modeled  by  means  of  a reduc^  bending 
stiffness,  the  output  value  of  “stress"  at  that  location 
may  be  misleading.  A better  estimate  of  stress  at  dis- 
continuities. may  be  obtainc ''  by  inspecting  the  vari- 
ation in  computed  stress  at  several  stations  adjacent 
to  the  crack  while  incrementally  moving  the  con- 
centrated loads  away  from  the  crack. 

Table  B2  is  used  to  specify  the  constants  for  the 
problem.  These  are  the  number  of  increments  in  the 
X and  Y directions,  the  increment  lengths  in  both 
directions,  and  Poisson's  ratio.  For  efficient  solution 
of  the  program,  the  number  of  Y increments  must  be 


equal  to  or  greater  than  tk<;  number  of  X increments, 
fable  B2  is  omitted  for  offspring  problems  since  the 
constants  must  be  the  same  as  in  the  parent  prob- 
lem. The  constants  specified  for  the  parent  are 
retained  and  used  by  all  successive  offspring 
problems. 

Data  for  optional  stiffness  constants  are  also 
input  in  Table  B2.  The  Modulus  of  Elasticity  E and 
the  slab  thickness  t must  be  specified  if  D’^,  Dy.  or  C 
are  to  be  internally  computed  by  the  options  in  Table 
Bt.  A uniform  subgrade  modulus  k must  be 
specified  if  S is  to  be  computed.  All  of  these  three 
constants  may  be  input  if  desired  and  thus  printed 
for  reference  without  the  options  exercised  in  Table 
Bl . The  options  cause  the  computation  of  the  equiv- 
alent data  which  could  be  input  in  Table  B4  based 
on  these  customary  plate  stiffness  relationships; 
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The  computed  stifl>;esses  are  constant  over  the  X 
by  Y area  of  the  slab  and  appropriate  quarter  and 
half-values  are  generated  at  comers  r.nd  ^ges.  They 
are  printed  as  an  extension  to  Table  B4,  and  thus 
any  additional  stiffness  values  entered  in  Table  B4 
are  superimposed  algebraically.  By  the  addition  of 
appropriate  reduced  lines  of  stift'nf-';^  in  Table  B4, 
any  degree  of  discontinuity  may  be  modeled.  Multi- 
valued subgrade  moduli,  or  foundations  with  voids 
may  also  be  modeled  by  entering  the  appropriate 
additive  spring  stiffnesses  in  Table  B4. 

Table  B3  is  used  to  define  the  lines  or  areas  of 
selected  output  for  deflection,  bending  moments  in 
the  X and  Y directions,  and  either  the  maximum 
principal  moment  or  stress  depending  on  the  print 
option  in  Table  Bl.  The  number  of  cards  are  as 
specified  in  Table  Bl  and  may  include  up  to  a maxi- 
mum of  10  cards.  Each  card  may  encompass  up  to  a 
m>*ximum  of  300  points.  That  is,  coordinates  from 
10.  10  through  20.  40  or  from  0,  0 through  12,  25 
could  be  specified.  If  a larger  area  is  required, 
another  card  covering  the  adjacent  area  could  be 
added. 
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The  major  advantage  of  Table  B3  is  that  a crude 
printer  plot  display  is  also  obtained  for  each  area 
specified.  This  feature  will  be  discussed  more  com- 
pletely under  the  heading  Computed  Results.  Table 
B3  is  especially  useful  when  studying  the  deRection 
and  moment  variations  along  a line  or  over  a local 
area.  Table  B3  may  be  omitted  if  desired,  since  all 
selected  output  value:  are  duplicated  in  the  complete 
printout  of  result*  in  Table  B8.  The  selected  output 
is  controlled  by  the  same  input  joint  coordinate 
system  described  below  for  Table  B4.  The  user  is 
cautioned  when  omitting  Table  B3.  since  he  may 
have  elected  to  suppress  the  complete  printout  of 
results  by  the  option  in  Table  Bl.  It  is  thus  possible 
to  use  a significant  amount  of  computer  time  and 
print  no  results. 

Table  B4  has  the  number  of  cards  as  specified  in 
Table  Bl.  Card  counts  should  be  carefully  checked. 
It  is  recommended  that  a listing  of  the  data  cards  be 
checked  by  the  user  prior  to  submission  of  the  pro- 
gram for  a run. 

Stiffness  and  load  data  are  entered  by  a coordi- 
nate system  notation.  The  coordinates  refer  to  the 
discrete-element  model  of  the  slab.  A joint  is  defined 
as  occurring  at  the  intersection  of  the  station  lines  in 
each  X and  y-direction.  A mesh  is  defined  as  that 
area  surrounded  by  four  joints.  A bar  b defined  as 
the  dbcrete-element  length  between  adjacent  joints. 
Figure  A2  of  the  Guide  for  Data  Input  summarizes 
this  notation.  Note  that  mesh  data  cannot  have 
either  a zero  x or  zero  y-coordinate;  x-bar  data 
cannot  have  a zero  x-coordinate,  and,  similarly,  the 
y-bar  cannot  have  a zero  y-coordinate.  If  the  data 
occur  only  at  one  location  (such  as  a concentrated 
load),  the  From  and  Through  coordinate  b simply 
repeated.  If  the  data  occur  along  a line,  the  coordi- 
nates will  reflect  this  by  having  either  both  x or  both 
y-coordinate$  the  same.  Data  dbtributed  over  a 
rectangular  area  are  specified  by  entering  the  lower 
left  and  upper  right  coordinates  of  the  area. 

Values  of  stiffness  entered  in  Table  B4  are  added 
algebraically  to  the  automatically  created  optional 
stiffnesses  computed  by  Tables  Bl  and  B2.  Thus,  th* 
user  has  the  ability  to  freely  model  any  dbcontinui 
ties  that  might  exist  in  the  actual  pavement. 

The  orthogonal  bending  stiffnesses  D^^  and  DX 
are  entered  in  each  direction  and  are  specified  on  a 
per  unit  width  basis.  If  the  edge  of  the  slab  coincides 


with  a station  line,  a half-value  of  stiffness  should  be 
input  for  both  D*  and  Df  along  the  edge.  If  the  edge 
of  the  real  slab  b not  on  a station  line,  a proportion- 
ate value  of  full  stiffness  is  entered.  This  is  demon- 
strated by  a sample  input  in  Figure  A3  of  the  Guide 
for  Data  Input.  The  stiffness  proportlonment  may  be 
thought  of  as  a direct  function  of  the  plan  area  of 
real  slab  surrounding  each  joint. 

Load  0 b concentrated  on  a per  joint  basis  and 
may  be  apportioned  at  each  joint  by  the  contributory 
area  load^  around  each  joint.  Positive  loads  act 
upward.  Loads  that  occur  between  joints  may  be 
fractionally  proportioned  to  the  adjacent  joints. 
Support  springs  S are  concentrated  values  input  and 
apportioned  exactly  like  loads.  A rigid  support  may 
be  specified  by  introducing  a large  value  of  support 
spring.  A maximum  value  of  1 x 1(P*  b suggested  to 
avoid  computational  difficulties  for  some  computers. 
Loads  may  also  be  applied  to  the  pavement  by  means 
of  specifM  load  pattern  in  Tables  B6  and  B7 
discussed  below. 

The  twisting  stiffness  C is  input  on  a per  unit 
width  basb  for  each  mesh  surrounded  by  four  joints. 
When  the  geometric  edges  of  the  actual  slab  do  not 
fall  on  a station  line,  proportionate  values  of  unit 
twbting  stiffness  may  be  input  similar  to  bending 
stiffheu  propoctlonment.  Computations  of  twisting 
stiffnesses  for  slabs  or  plates  are  at  best  still  approxi- 
mate procedures.  This  b due  to  uncertainty  in  the 
defining  of  the  shearing  modulus  of  rigidity.  The 
best  procedure  is  to  ascertain  the  twbting  stiffness 
experimentally  as  outlined  by  Reference  10.  The 
formulas  shown  above  in  the  Table  B2  discussion  are 
correct  for  uniformly  thick  isotropic  plates.  An 
approximate  value  of  twbting  stiffness  for  ortho- 
tropic slabs  or  stiffened  plates  may  be  obtained  by 
using  procedures  outlined  by  Huffington**  or  com- 
putati^  summarized  by  Troibky  in  a recent  publi- 
cation.** Fortunately,  precise  values  of  twisting  stiff- 
ness are  unnecessary  to  model  a slab.  The  main  load 


"W.  R.  Hudson  and  Hudson  Matlock,  DiMcontinuous 
Orthotropk  PUtf  mi  favtmtnt  Sisks,  Rsssaich  Rspoit 
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"N.  J.  Huffiagton,  Jr.,  Tkeontktl  Dttwmiitaiion  of 
Rigidity  of  Orihogomdiy  StifftHod  Pious 

(Appitsd  Msehanics  Division,  ASME,  1966). 
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JOINT  2,4 


JOINT  data:  D*,d’o,S 

(0‘anD  d’  are  per  unit  width,  0 AND  S ARE  CONCENTRATED) 

MESH  data:  C 

(C  1$  PER  UNIT  WIDTH) 

BAR  data;  p“,p^ 

(P*  AND  P’  ARE  CONCENTRATED) 


I1g«r«  Al.  0«U  eoordinato  aumlwriag  •ystem. 


Ji 


i n u 1 1 u 11  i 1 1' 


Vioh-.| 


UNIFORM  AXLE  COMPRESSION 
Y DIRECTION  OF  600  lb  ptr 
INCH  OF  WIDTH 


IM.  I 


• hi.  •«( 


r 

1 

L 

'■  1 

1 

I 

1 

\ 

1 

1 

1 

I 

1 

■4 

1 

1 

•••hi. 


UNIFORM  LOAD  0>20pti 
EQUIV.  PER  JOINT  LOAD 

0«  20-^06110)  in.‘ 
>1200  ibi 


PLATE  STIFFNESS 
D*  AND  0^«300-M^ 

ifl* 


•s-kip  CONCENTRATED  LOAD 


t ft  t ft  t t ft  t ftt 


FROM 

X y 


TNROUQH  0*  AND  0^  0 


7.500E«04 

7.5001404 

1.500B405 

9.000B404 

9.0002404 


-1.200E403 

-3.600»f02 

-8.000E403 


data  incoapiata  for  this  aaapla 


-3.0002403 

-6.00bE403 

-3.600E403 


FNpw  A8.  Baaph  Mil  hi|tat. 


X 


carrying  capabilities  of  a slab  are  due  to  its  bending 
stiffness  which  is  more  accurately  definable. 

By  using  a zero  value  of  twisting  stiffness  and  zero 
Poisson's  ratio,  a simple  grid  system  may  be 
modeled.  Each  beam  would  be  modeled  by  an 
appropriate  D*  or  Df  term  which  would  then  be  the 
per  unit  width  stiffness.  The  beam  stiffnesses 
entered  should  therefore  be  divided  by  the  increment 
width. 

Table  BS  is  for  input  of  in-plane  axial  tensions 
or  py  if  present.  These  might  be  generated  due  to 
temperature  differentials  or  traffic  braking  and 
acceleration  forces.  There  is  no  provision  in  the  pro- 
gram for  automatic  distribution  of  applied  axial 
forces  sinee  no  in-plane  supports  are  used  which 
would  restrain  them.  The  user  must  specify  the  dis- 
tribution of  the  axial  tensions  (■¥)  and  compressions 
( I in  each  x-bar  and  y-bar  of  the  model.  Since  these 
ate  oar  forces,  no  data  should  be  input  which  would 
represent  forces  outside  the  boundaries  of  the  actual 
slab.  A brief  sample  of  data  input  is  given  in  Figure 
A3  in  the  Guide  for  Data  Input. 

All  data  in  Tables  B4  and  B5  are  algebraically 
accumulated  and  values  therefore  may  be  added  or 
subtracted  regardless  of  other  values  specified. 

Table  B6  is  for  input  of  special  load  patteras  such 
as  are  found  on  the  C-SA  or  Boeing  747  aircraft. 
This  table  allows  the  user  to  specify  one  pattern 
which  can  then  be  repeated  or  placed  several  times 
on  the  same  pavement.  By  this  means,  the  placement 
of  wheel  load  patterns  by  individual  load  inputs  in 
Table  B4  is  unneoesury. 

A maximum  of  nine  different  patterns  may  be 
specified  with  a maximum  of  12  loads  in  each 
pattern.  Each  pattern  requires  two  input  data  cards. 
The  first  card  specifics  the  pattern  number  designa- 
tion. the  number  of  loads  in  the  pattern,  and  the 
pattern  coordinates.  The  pattern  has  its  own  local 
coordinate  system,  the  only  restriction  being  that  the 
increment  spacing  is  the  same  as  for  the  complete 
pavement  described  by  Table  B2.  The  second  card  of 
the  pattern  is  for  the  entry  of  the  indWiduai  pattern 
load  magnitudes.  For  loads  between  stations,  a geo- 
metric apportionment  of  load  could  be  used.  The 
first  pattern  coordinate  and  load  entered  in  columns 
9 through  14  b the  reference  load  for  the  pattern  and 
U the  load  which  b placed  by  Table  B7.  It  may  have 


any  coordinate  including  negative  values  and  may 
also  have  a zero  magnitude  of  load  which  is  some- 
times convenient  A typical  pattern  and  Table  B6 
input  is  shown  in  Figure  A4  of  the  Guide  for  Data 
Input. 

For  a problem  series  in  which  the  effects  of 
several  load  patterns  may  be  studied,  all  patterns 
could  be  entered  in  Table  6 of  the  first  problem.  The 
selected  pattern  numbers  are  then  applied  to  the  slab 
by  the  Table  B7  of  each  problem.  The  series  can  be 
parent  and  offspring  problems,  or  individual  prob- 
lems in  which  stiffness  might  be  changed  from  prob- 
lem to  problem.  New  patterns  can  be  entered  at  any 
problem  in  the  series  and  are  then  available  for  sub- 
sequent problems.  If  the  same  pattern  number  is 
used  as  was  previously  used,  then  the  new  pattern 
replaces  the  old  pattern. 

Table  B7  is  for  the  placement  of  any  or  all  of  the 
load  patterns  defined  by  Table  B6.  Any  number  of 
pattern-placement  cards  may  be  used.  Placements  of 
the  reference  load  are  by  coordinates  which  are 
related  to  the  overall  system  defined  by  Table  B2. 
Placements  outside  the  boundaries  of  the  actual  slab 
arc  permissable.  Any  loads  of  the  pattern  that  thus 
fall  off  the  slab  are  not  included  in  the  solution.  The 
same  pattern  number  may  be  repeated  on  several 
cards  if  desired.  They  need  not  be  in  numerical 
ascending  'sequence. 

For  offspring  problems,  only  Table  Bl  and  Table 
B4  with  added  lowls  are  required.  Any  loads  entered 
in  the  offspring  Table  B4  replace  all  other  previous 
Table  B4  loads.  Offspring  problems  may  also  be  run 
with  only  Table  Bl  and  new  load  pattern  placements 
in  Table  B7.  Table  B2.  stiffnesses  in  Table  B4,  and 
Table  BS  cannot  be  used  in  offspring  problems. 
Tables  B3,  B6  and  B7  may  be  used  as  desired. 

Input  CtMCkllft.  An  input  checklist  is  provided  in 
Table  A I of  the  Guide  for  Data  Input.  A listing  of 
the  input  data  prior  to  problem  submission  will  also 
aid  in  detecting  coding  and  keypunch  errors. 

Otta  Errert,  All  data  are  checked  for  compatibility 
with  the  geometry  of  the  specified  slab  and  con- 
sbtency  of  coordinate  input.  A count  of  the  number 
of  data  erron  is  made  arid  if  any  errors  are  encount- 
ered, the  problem  is  terminated  and  a message 
showing  the  number  of  data  errors  made  is  printed. 
Typical  errors  ate;  I.  misuse  of  the  multiple-load 
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Table  BC  Input  for  this  pttttrn; 
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ri|«r«  A4.  Tjrptcd  iiwdN  kad  peUNR. 
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option,  such  «s  a - 1 following  a 0 in  the  preceding 
problem;  2.  the  number  of  increments  in  the  a- 
direction  exceeding  those  in  the  y-direction  (if  this 
were  allowed,  an  ineflicient  and  time  consuming 
computer  solution  would  result);  3.  a negative  or 
zero  increment  length  specified;  4.  a negative 
Poisson's  ratio  or  thickness  input;  5.  the  “Through" 
X or  y-coordinate  in  a data  specification  numerically 
less  than  the  “From"  coordinate;  6.  data  specific 
outside  the  geometric  limits  of  the  slab  (exceptions 
arc  Table  Bb  and  B7  which  allow  load  patterns  to  be 
outside  slab  boundaries);  7.  a zero  x or  y-coordinate 
specified  for  a twisting  stiffness;  8.  a zero  x- 
coordinate  used  for  x-bar  axial  tensions  or  a zero  y- 
coordinate  for  y-bar  tensions;  9.  the  number  of 
increments  specified  are  greater  than  the  dimen- 
sioned storage  the  program  can  operate  with;  and 
10.  misuse  of  the  selected  output  option, 

Computdr  fiMUitt.  All  input  data  cards  are  re- 
flected in  the  printout  of  results  exactly  as  they  were 
input.  It  is  good  practice  to  again  check  these  data 
for  possible  errors  prior  to  inspection  of  the  re- 
mainder of  the  results.  The  importance  of  using 
descriptive  alphanumeric  information  for  the 
problem  series  header  cards  and  for  each  problem 
number  card  cannot  be  over-emphasized.  *rhis  will 
avoid  confusion  when  running  a large  number  of 
problems. 

A continuation  of  Table  B4  is  printed  after  the 
normal  Table  B4  input  which  reflects  the  computed 
optional  pavement  stiffness  constants.  Appropriate 
quarter,  half,  and  full  values  are  generated  for  bend- 
ing stiffnesses  and  the  subgrade  spring  equivalent. 
Full  values  of  twbting  stiffness  are  compuM  for  all 
the  mesh  areas  of  the  slab.  Reference  to  Figure  A2  in 
the  Guide  for  Data  Input  explains  the  different  data 
types  and  their  associated  numbering  systems. 

Immediately  after  the  echo  print  of  load  place- 
ments in  Table  B7,  a summary  of  the  applied  load 
patterns  is  printed.  Thb  should  be  inspect^  to  see  if 
the  total  number  of  wheel  loads  and  their  algebraic 
sum  are  of  the  intended  amounts.  The  number  of 
pattern  wheels  which  fall  outside  the  slab  boundaries 
and  their  algebraic  sum  are  also  printed  for 
information. 

Table  B8  gives  the  complete  printout  of  final 
results.  If  the  option  to  delete  Table  B8  b exercised 
in  Table  Bl,  then  only  a single  reference  station  is 


printed  at  or  near  the  center  of  the  slab.  Again,  the 
user  b cautioned  to  use  this  option  with  discretion  to 
avoid  deleting  significant  results.  Table  B8  is 
arranged  to  give  the  x and  y-joint  coordinate,  the 
transverse  deflection  at  each  joint  (upward  deflec- 
tions are  positive),  the  bending  and  twisting 
moments,  the  support  reaction,  the  principal 
moment  or  stress  and  its  direction.  The  tabulations 
are  arranged  in  groups  for  each  y-station.  Output 
values  of  bending  and  twisting  moments  are  given  on 
a per  unit  width  basb.  Bending  moments  arc  positive 
for  compression  in  the  top  of  the  plate  or  slab.  The 
x-bending  moments  act  in  the  x-direction  and 
y-bending  moments  in  the  y-direction. 

The  per  unit  width  x twisting  moment  b tabu- 
lated and  is  exactly  equal  to  the  y twbting  moment 
with  opposite  sign.  The  x-twbting  moment  acts  in 
the  x-direction  and  is  about  the  y-axis.  Even  though 
the  input  values  of  twisting  stiffness  were  specified  at 
each  mesh,  the  output  value  of  twisting  moment  is 
the  average  of  four  adjacent  mesh  areas  and  is  there- 
fore given  at  the  joint.  The  user  b cautioned  that  the 
output  values  of  twisting  moment  along  the  edges  or 
other  discontinuities  of  a slab  or  plate  reflect  the 
average  and  therefore  may  be  a one-quarter, 
one-half,  or  some  other  prop^irtionate  value  of 
twbting  moment. 

The  support  reaction  is  the  concentrated  value  of 
resistance  to  displacement  offered  by  any  support 
springs  if  present.  A subgrade  modulus  spring  will 
reflect  the  concentrated  value  of  pressure  under  the 
slab.  The  value  of  support  reaction  is  different  than 
the  similar  value  tabulated  in  prior  versions  of 
DSLAB  and  SLAB  computer  programs.  In  those 
programs,  the  value  might  have  been  better  labeled 
Net  Reaction,  or  Net  Force,  since  it  was  the  summa- 
tion of  the  applied  load  and  support  reaction. 

Internally  in  the  program,  a Mohr's  circle 
analysb  b made  at  each  joint  using  the  orthogonal 
bending  moments  and  the  twisting  moments  to  yield 
the  larger  numeric  value  'positive  or  negative)  of 
principal  moment  per  unit  width  and  the  angle  from 
the  x-axb  of  the  coordinate  system  to  the  acting 
direction  of  this  larger  value.  Counterclockwise 
angles  are  positive.  If  the  stress  option  is  specified  in 
Table  Bl,  the  values  are  converted  to  the  larger 
numeric  value  of  principal  stress  instead  of  moment. 
A positive  stress  indicates  tension  in  the  bottom  of 
the  slab.  The  input  value  of  thickness  b property 
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used  only  for  slabs  of  constant  thickness.  For  slabs  of 
variable  stiffness  and  thickness,  a direct  conversiofi 
can  be  made  for  principal  stress  from  the  principal 
moment. 

As  both  a check  on  the  internal  computer  solu- 
tion and  a check  on  the  input  of  the  load  system  to  a 
problem,  a statics  check  is  printed  at  the  end  of 
Table  B8  which  is  the  algebraic  sum  of  all  the  re- 
action values  and  should  be  equal  to  the  sum  of  all 
the  applied  loads  from  Tables  M and  B7.  This  check 
should  always  be  inspected  to  verify  that  the  desired 
load  system  was  specified  and  that  the  problem  was 
properly  solved. 

Tabic  B9  has  been  found  to  be  the  most  immedi- 
ately useful  type  of  output  in  the  SLAB  30E  pro- 
gram. It  lists  the  output  parameters  which  were 
specifled  as  selected  output  in  Table  B3.  Each  area 
of  Table  B3  is  treated  as  a separate  output  array  for 
each  of  the  four  desired  output  values  of  deflection, 
bending  moments  in  the  X and  Y directions,  and 
principal  moment  or  stress  depending  on  the  Table 
Bi  streu  option. 

The  coordinate  references  are  printed  in  consecu- 
tive groups  associated  with  the  largest  number  of 
increments  in  the  rectangular  area.  For  instance,  if 
an  area  were  desired  from  10,  IS  to  12,  20  then  there 
would  be  three  groups  of  values  each  with  five 
values.  If  the  area  is  square,  then  the  groups  are  for 
each  Y value  included  which  is  the  same  as  the 
normal  arrangement  of  Table  B8  output.  Adjacent 
to  the  coordinate  is  the  numerical  value  of  the  ^fkc- 
tion,  montent,  etc.  To  the  right  of  the  .output  values 
is  a series  of  uterisks  which  have  a relative  place- 
ment to  one  another  which  are  proportional  to  the 
numerical  output  value.  Thus,  a crude  plot  of  the 
output  values  is  obtained.  The  plot  has  a width  of  20 
printer  characters  and  a length  equal  to  the  number 
of  nodes  encompassed  by  the  area  specifled.  The  plot 
has  no  scale  and  no  zero,  the  values  are  relative  to 
one  anotiter.  increasing  positively  to  the  right.  The 
20  character  width  is  based  on  the  minimum  and 
maximum  value  in  the  area  to  be  plotted.  The  user  is 
cautioned  to  not  misiaterprtt  apparant  changet  in 
plot  curvatures  which  might  have  been  generated 
due  to  very  slight  numerical  change. 

The  plots  have  been  found  to  bo  especially  valu- 
able in  understanding  slab  behavior  in  areas  adja- 
cent to  concentrated  wheel  loads.  All  deflection 


areas  are  printed  first  followed  by  bending  moment 
areas  which  are  printed  adjacent  to  each  other  if 
both  X and  y-moments  were  desired  in  the  same  area. 
The  final  selected  output  is  for  the  principal 
moments  or  stresses.  Plots  of  principal  stresses  along 
slab  lines  are  somewhat  misleading  since  the  direc- 
tion of  the  stress  usually  varies  along  the  Hne.  it  is 
valuable,  however,  in  pointing  out  maximum  values 
of  stress  which  might  overlooked  when  inspecting 
a mau  of  numbers  in  the  normal  Table  B8  output. 

The  final  printed  output  is  of  the  computer  time 
which  has  bMn  used  for  the  problem  plus  another 
time  which  is  the  total  elapsed  time.  The  user  should 
record  his  problem  run  times  for  parent  and  off- 
spring problems  for  each  different  problem  size  he 
runs.  By  this  means,  an  estimate  of  required  run 
times  can  be  made  for  future  problems.  For  s*  all 
problems,  the  offspring  times  will  be  from  20  to  50 
percent  of  the  patent  problem  times.  Fortunately, 
the  offspring  problem  time  decreases  to  a very  smaM 
proportion  of  the  parent  problem  time  as  the  prob- 
lem size  becomes  large.  A time  as  low  as  four  percent 
is  possible.  For  this  reason,  the  user  should,  for 
computer  economy,  run  as  many  offspring  problems 
for  each  parent  problem  as  he  thinks  he  might  need. 

ComfMtar  Timd  RtquIrMiMnU.  The  computer 
time  required  for  running  problems  with  program 
SLAB  30E  cannot  be  as  precisely  defined  as  the  stor- 
age requirements  indicated  in  Figure  Al.  As  pre- 
vioutly  stated,  it  would  be  good  practice  for  the  user 
to  record  run  times  as  he  uses  the  program  to  be  able 
to  make  better  run  time  predictims. 

Based  on  a small  sample  of  run  timet  using  SLAB 
30E  and  similar  programs,  the  following  time 
exprenions  were  derived: 

For  patent  and  independent  problems, 

Where  Tp  is  the  predicted  central  processor  time  in 
seconds,  and  X and  Y are  the  number  of  slab 
increments. 

Per  oflkpring  problems, 

To-Tp(5J^  (EqA3) 
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whi-n;  T„  h>  the  ofTsprinit  lime  and  Tn  is  the  ab«we 
piirent  problem  time.  A minimum  ot  UHir  percent  ttl' 
Tp  should  be  reserved  for  Tg  times  for  large 
problems. 

in  both  parent  and  offspring  problems,  a very 
rough  estimate  of  input-output  time  requited  may  be 
found  by  this  expression: 

Tl/0»^  (EqA41 

where  T|/q  is  the  estimate  of  input-output  time  in 
seconds. 

The  above  time  estimates  should  be  reasonably 
accurate  for  problems  with  X less  than  30.  For  larger 
problems,  the  estimated  times  will  probably  be 
somewhat  conservative. 

The  program  compile  time  must  be  added  to  the 
estimate  problem  run  times.  The  complete  pro- 
gram deck  wilt  compile  in  approximately  ten 
seconds.  If  stored  on  tape  in  a re-locatable  form,  the 
compile  time  will  be  nominal. 

Smnpln  Preblvins. 

Bridgt  Approach  Slab-Problam  601M,  This 
problem  is  the  ume  as  example  problem  601  of 
Stclxer  and  Hudson's  report'*  and  similar  to 
problem  601  of  Hudson  and  Matlock's  report.**  It  is 
discussed  here  to  demonstrate  the  use  of  Program 
SLAB  JOE  and  to  give  a reference  solution  for  pro- 
gram checkout  when  necessary, 

The  problem  is  illustrated  in  Figure  AS.  A 
lO-inch-thick,  reinforced-coocretc  bridge  approach 
slab  is  used.  It  is  supported  on  one  end  by  the  bridge 
abutment;  the  othw  end  rests  on  the  embankment. 
Because  of  poor  compaction,  which  often  results 
when  there  is  backfill,  the  soil  has  settled  under  the 
interior  of  the  slab  and  left  a section  unsupported. 
The  stab  has  a center-tine  joint  and  a crack  which 
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developed  from  a combination  of  shrinkage  and  pre- 
vious overstress.  For  .-i  nun-uniformly  supported  slab 
such  as  this,  the  dead  weight  of  the  slab  must  be  con- 
sidered when  evaluating  moment  and  stresses.  This 
weight  acts  as  a uniform  load  of  600  lb  per  station. 
Two  10-kip  wheel  loads  were  considered  in  this 
example.  An  axial  load  of  5000  lb  per  in.  has  been 
indut^  by  the  expansion  of  the  adjoining  pavement. 

The  problem  is  modeled  by  a 12  by  16  increment 
slab.  The  slab  stiffness  options  are  exercised  in 
Table  Bl.  The  subgrade  is  introduced  in  Table  B4 
over  the  proper  area  from  stations  0,  7 to  12,  IS.  The 
crack  and  joint  are  modeled  by  subtracting  the  full 
bending  stiffness  along  tlieir  lines  from  0,  7 to  12,  7 
and  6, 0 to  6.  16.  The  abutment  support  is  simulated 
by  a double  line  of  rigid  supports  along  that  edge 
from  0.  0 to  12,  I.  The  loads  are  introduced  by 
Tables  B6  and  B7  to  demonstrate  their  use,  although 
they  could  have  been  coded  in  Table  B4.  it  will  be 
noted  that  there  are  two  load  patterns  in  Table  B6 
which  simulate  a truck  and  trailer.  The  load  place- 
ment in  Table  B7  is  such  that  only  two  of  the  trailer 
wheels  fall  on  the  approach  slab  in  the  same  position 
as  shown  on  Figure  AS. 

The  reader  will  note  that  the  computed  results  are 
slightly  different  than  Stelzer  indicates.  This  is  due 
to  a coding  error  in  his  problem.  The  illustration  of 
Figure  AS  shows  an  axial  compression  of  5000  Ib/in. 
in  the  y-direction  of  the  slab,  but  the  problem  as 
coded  in  the  Stelzer  report  has  this  axial  compres- 
sion in  the  x-direction.  The  thrust  is  slight  and  has 
little  effect  on  the  results  in  either  case  however. 

The  results  for  problem  601 M are  also  somewhat 
misleading  near  the  face  of  the  abutment.  A double 
line  of  rigid  supports  was  used  to  represent  the  abut- 
ment support  which  effectively  fixes  the  slab  against 
rotation  at  that  end.  A more  realistic  modeling 
would  have  been  achieved  by  a single  line  of  supports 
at  the  abutment  edge,  thus  allowing  rotation  to 
occur.  The  amount  of  nominal  anchorage  that  is 
utuaHy  placed  between  slab  and  abutment  is  not 
nearly  enough  to  preclude  rotation. 

Taxiway  Slab  with  C-5A-Problem  JRS.  This 
and  the  next  problem  show  how  the  application  of 
Tablet  B6  and  B7  can  be  used  to  conveniently  code 
complex  load  patterns  and  their  placements.  The 
problem  it  lliuitraled  in  Figure  A6  with  the  heavy 
crosses  indicating  the  C-SA  wheel  loads,  and  the 
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darkened  circles  the  Boeing  747  loads  which  are  used 
in  the  following  problem. 

The  slab  area  under  study  has  four  ZSToot  square 
panels  with  cracks  between  the  panels  that  have  25 
percent  fleiural  stiffness  remaining.  The  slab  is  12 
inches  thick,  the  subgrade  modulus  is  125  lbs  per 
cubic  in.,  and  the  coiKrete  modulus  of  elasticity  is 
five  million  lbs  per  square  in.  with  a Poisson’s  ratio 
of  0.2. 

The  selected  increment  lengths  of  30  in.  and  21.5 
in.  are  coarse  for  a problem  of  this  type,  and  the  user 
would  probably  use  more  increments  to  effectively 
model  his  problem.  However,  a coarse  solution  when 
inspected  would  allow  the  user  to  then  possibly 
delete,  for  additional  finer  mesh  problems,  areas  of 
the  slab  which  obviously  do  not  affect  the  area  under 
study. 

All  the  input  and  output  options  are  specified  in 
Table  Bl.  A selected  output  area  which  includes  the 
loaded  area  of  the  slab  is  specified  in  Table  B3.  In 
addition,  three  lines  of  selected  output  are  specified 
in  the  x and  y-directions.  Seventy-five  percent  of  the 
slab  bending  stiffness  is  removed  by  Table  B4  input 
along  the  lines  of  tlie  joint.  It  is  noted  that  this  is 
accomplished  by  removing  37.5  percent  along  the 
complete  line  and  the  remaining  37.5  percent  at  the 
interior  stations  along  the  joint.  By  this  means,  the 
removed  stiffneu  at  the  edges  is  compatible  with  the 
automatically  generated  half-values  at  the  edges. 
One  six-wheel  gear  group  of  the  C-5A  is  described  in 
Table  B6  with  an  additional  reference  wheel  at  the 
center  of  the  pattern  with  a zero  load  value.  This 
pattern  is  then  placed  twice  by  Table  B7  to  create  the 
total  load  pattern  as  shown  on  Figure  A6. 

The  output  of  Table  B8  has  been  omitted.  The 
output  of  Table  B9  indicates  that  the  maximum 
deflection  occurs  at  station  10,  10  and  the  maximum 
streu  occurred  at  station  9,  11.  The  large  slab  area 
specified  in  Table  B3  from  3, 8 to  17,  21  allows  these 
maximums  to  be  immediately  discerned  by  the  rela- 
tive plot  values  adjacent  to  the  tabulated  values. 
Three  additional  lines  of  selected  output  were  speci- 
fied by  Table  B3  which  coincide  with  some  of  the 
wheel  loads  as  shown  on  Figure  A6.  These  lines  are 
valuable  fix  comparing  deflection  and  moment  vari- 
ations in  specific  directions  or  areas.  It  may  be  noted 
that  the  plotted  output  from  areas  can  be  thought  of 
as  plots  of  contiguous  strips  with  the  relative  values 


of  each  line  in  proper  proportion  to  adjacent  lines. 
Variations  of  deflection,  stress,  or  bending  moment 
may  thus  be  studieti  in  one  direction  along  consecu- 
tive strips  taken  in  the  other  direction,  but  still  allow 
tbc  maximums  to  lie  emphasized  by  the  plots. 

Taxiway  Slab  with  B-747— Problem  JR6.  This 
problem  is  the  same  slab  shown  in  Figure  A3  which 
was  used  for  Problem  JR5.  A different  load  pattern 
is  used  which  represents  two  gear  groups  of  the 
Boeing  747  aircraft  Due  to  the  modeling  coarseness 
of  this  demonstration  example,  two  wheels  of  each 
group  are  split  to  adjacent  node  points,  thus  forming 
a six-wheel  group  for  each  pattern.  This  apportion- 
ment must  be  taken  into  consideration  when 
interpreting  the  results. 

This  problem  is  run  as  an  offspring  problem  to 
the  previous  problem.  Only  Tables  Bl.  B3,  Bb  and 
B7  were  input.  The  pattern  could  have  been 
described  as  pattern  No.  2 in  the  previous  problem 
and  then  simply  placed  by  the  Table  B7  of  this  prob- 
lem. The  computer  efficiency  of  multiple  loadings 
via  offspring  problems  may  be  seen  by  comparison  of 
tfte  run  time  of  23  seconds  in  Problem  JR5  to  7 
seconds  in  this  problem  which  is  30  percent  of  the 
parent  problem  time.  For  extremely  large  problems 
(which  might  be  a maximum  of  68  by  68  or  20  by  342 
as  seen  in  Figure  Al)  the  offspring  time  might  be  as 
little  as  3 to  4 percent  of  the  parent  problem's  time. 
A procedure  for  estimating  computer  times  was 
given  in  the  discussion  of  the  program. 

The  output  for  this  problem  may  be  compared  to 
that  of  Problem  JR5.  The  major  bending  moments 
with  the  B-747  loads,  even  though  the  individual 
wheels  are  of  almost  SO  percent  greater  magnitude, 
are  not  50  percent  greater.  When  the  stresses  (or 
principal  mon.'ents)  arc  compared,  however,  the 
B-747  maximum  is  significantly  higher  than  the 
stress  induced  by  che  C-5A.  This  is  due  to  the  much 
higher  twisting  mon>ents  present  which  are  evidently 
caused  by  the  influei.'ce  of  the  two  gear  patterns  on 
each  other.  The  two  C-5A  patterns  are  orthogonal  to 
the  major  slab  directions,  but  the  two  B-747  patterns 
are  oriented  almost  at  45  degrees  to  each  other  in 
relation  to  the  major  slab  directions,  thus  causing  a 
significant  tw'isting  effect. 

Simply  Supported  Plata— Problems  A C 1 . This 
problem  is  added  for  reference  and  to  roughly 
demonstrate  the  accuracy  of  the  program.  The  prob- 
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Icni  IS  a simply  supp4ir1cd  isotropic  plate  48  imhes 
square.  ('.9788  inches  thick,  with  a concentrated 
load  ot  10‘  pounds  at  the  center  The  optional  stiff- 
ness  input  was  not  used,  although  it  could  have 
l>cen.  Ihc  stillnesses  entered  in  lable  B4  arc  com- 
puted from  the  above  plate  thickness  for  a modulus 
ofclasticity  of  3.0x  lOMbs/in^  and  a Poisson's  ratio 
of  0.25.  The  simple  supports  arc  created  by  very  stiff 
support  springs  around  all  edges. 

The  results  indicate  a maximum  deflection  at  the 
center  of  1 . l.)8  inches.  Fven  though  this  is  a coarse 


mesh  ol  only  8 by  8 increments,  the  dcncction  is  seen 
to  agree  closely  with  a closcd-lorrn  solution  of  1 .07 
ini-b*^  given  by  Timoshenko. 

This  problem  is  recommended  to  be  run  lor  quick 
program  checkout  whenever  it  becomes  necess?*-“ 
An  8 by  8 driver  for  the  program  will  only  require 
approximately  24,(XX)  (56.(XXJ  octal)  words  of  storage 
and  the  problem  will  run  in  about  12  (M  octal) 
seconds  which  includes  10  seconds  for  a complete 
recompilation. 
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mPFENDIXB  GUiOE  for  data  input 


MAHWII  (.(Itl)l  KOKOAIA  INI’Ur  CAROFORMS 
II  )1  N I II  H A I ION  ()^  Rl  N (iwi.  .i!|ili.iniinirru'  I'atris  (x’f  mot 


I ' ^ 


IDi'N  III  l(  A I ION  or  I’ROBI.l'.M  (urn- card  lor  each  problem:  program  Mop«  if  PROB  NUM  is  IcIt  blank) 


I’ROIINUM 

1 ^ ..L.  . 

pescriptioti  ol  problem  (alphanumerie) 

' ■■■  1 

1 5 

II 

80 

fable  Ml.  CONTROl,  DATA  (one card (r reach pmbEcml. 


Multiple  load  Option 

Number  of  Cards 

Optional  Stillness 

Input  Constants 

Options  1 1 
Delete 

Yes) 

Print 

* I lot  Parent  probs 

(or  Table 

Enter  1 to  compute 

Dcl.iilcd 

Principal 

1 lorOllsprings  2 

1 4 5 (>  7 

D**!))'  C S 

Out  pill 

Stress 

1 1 rm 

■ TT  " T 1 

••  1 **  ”1 

5 20 

25  .30  35  40  45 

55  W bS 

75  80 

Table  B2.  CONSTANTS  (one  card,  o«H  forodiiHrini  ptabiMit). 


Increment  Lengths  Person's 

Ntim  Iners  X-Direction  Y-Direction  Ratio 

X Y h,  hy  V 

t— — Data  for  Optional  Stiffness  Constani  —a 

Modulus  of  Slab  Subgrade 

Elasticity  Thickness  Modulus 

F.  t k 

. .nj 

] 

[ __L.  - J 

' :3 

s 10  20  Ml  40  50  ()0  70 


I I,  - F/I.J  I.  hVl.J 

f Korcc  Onus.  I.  lx-n)!lh  Units.  All  units  nuisi  1h’  consistent 


Table  B3.  SPKCIFIEU  AREAS  FOR  SEl.ECTEO  OUTPUT  (namberofrarde  aiihown  InTaMt  BJ  lOmaihaum) 


from  IhrouRh  PRINT  and  SPI.OT  (I  Yes) 

X Y X Y UF.FL  X MOMY  MOM  PRIN  STRF.SS/MOMENT 

I ■ T""xiri:i:zn::izT~ri 

5 in  15  20  25  JO  35  40 


Table  84.  STIFFTIESS  AND  LOAD  DATA  (number  of  caedt  aa  ahasm  In  Table  B1  enter  onl;  toad  for  offtprini  ptobleina) 


from  Ibroutih 

X 3 X Y 

Cl 

5 10  15  20 


Bending  Stiffness 
U»  


Load  Spring  resisting  Stiffness 

o s c 


.JO 


40 


50 


W) 


70 


’for  otlspii  g problems  J abirs  B2  and  BS  must  be  omitted,  amt  only  load  may  l>c  entered  m I able  B4  I ables  B.l  Bh  .iiid  B7  mas  Ik- 
iisi-d  4-^  desmd 

••Opinnal  Siillm-ss  Saiuhes  nuisi  be  blank  lor  olTspring  problems 


4i 


f 


mmmm 


Tabt*  ns.  AMAl.  THRUST  DATA  (namber  •feanb  m rfMwn  In  TaMc  11  Mnk  ft  nlfayriwc  ptoblww). 


I'r  nil 


Tluougll 
X Y 


Anial  ^^lru^l 
,.x  (.y 


T ■ mm 

5 10  IS  20 


izzm 

70 


III 


Ml 


TfWt  BA.  SPf.n  AT  IXJAO  PATTERNS  lanMlMr  nf  cm*  m tbtni  In  Tnblt  Bt  a|»  to  9 pnltonw  nuy  b*  uted). 


I’.it  Ntini  of  These  coordinates  refer  to  the  load  pattern 

^.•rn  loads  m ^ Reference  Load 

Nimi  Pattern  ,yx)*y«y*y*y«y«y*y»y*y') 


i.a. i~i-r-n-i  i i i i i i i i i i i i i i ii~m 

I h It  I4  I7  20  2.1  2h  20  J2  2S  38  4 I 44  47  SO  S3  56  59  62  65  68  7|  74  77  M) 


Magnitude  ol  Loads 


I' m"  i :::i: 


rzT  rm 


I4 


20 


26 


32 


38 


50 


56 


62 


68 


74  Ml 


S(Hxi.il  P.'illerns  arc  placed  ass(Kei(icd  in  Table  B7  Up  loOdifTercnl  patterns  (18  cards)  may  be  used.  The  first  load  of  each  pattern  is  Ihc 
ivtcrence  lor  I able  B7  ii  may  hare  a aero  load  magnitude  Negative  pattern  coordinates  arc  pcrmissabic 


Table  B7.  PLACEMENTS  OF  TABLE  B6  LOAD  PATTERNS  Innnbct  nf  etda  m ahnwn  In  TnUt  Bt), 


Pattern  Numot 
Ntini  Place 

mcnis  XYXYXVXYXVX 

L 1--1  I i--r— I I I I 'I T-~r 

5 10  IS  20  25  30  35  40  45  SO  55  60 


X Y 


im 

70  75 


80 


Ibssc  coordinates  refer  to  the  coordinate  system  defined  by  Table  B2.  Negative  coordinate  placements  arc  pcrmissabic. 
< .1  NP.K At,  PROGRAM  NO  I ES 


The  data  cards  must  be  assembled  in  proper  order  for  the  program  to  run. 

A consistent  system  of  units  must  be  used  for  all  Input  data,  for  example,  kips  and  feel. 

All  2 to  S space  words  arc  undcrsi  id  to  be  right-justified  integers  j ^ 1 1 

All  lO-space  words  are  floating-point  numbers  (EI0.3I I — 4 . 3 2 I fc  T oT] 

All  6-space  words  are  noalmg-poim  numbers  (F6.0) | — 3 0 . 0 [ 

Table  Bl . CONTROL  DATA. 

Ihc  multiple  load  option  is  ciercised  for  problem  series  In  which  only  the  load  positions  and  magnitudes  of  Tables  84  or  86  and  will 
vary  The  llrsi  problem  m a series  is  the  Parent  and  is  specified  by  cnterini  + i.  successive  loadings  are  the  Offspring  and  are  specified 
by  entering  —I.  If  the  option  is  left  blank.  Ihc  problem  is  complete  within  itself. 

1 he  number  ol  cards  input  for  Table  B2  through  Table  B7  should  be  carefully  checked  after  coding  is  completed.  A checklist  for  allowable 
I ables  and  numbers  of  cards  lor  each  Table  is  given  in  Table  At. 

< ..loeniumal  plate  or  ilab  input  stiffness  may  be  selected  by  entering  1 for  those  to  be  computed  Appropriate  vsluw  must  then  lie  made 
avaitahir  in  I able  B2 

tbiipiii  ..piuiio  may  ai«i  be  selected  d desired  All  detailed  output  of  Tabic  B8  may  be  suppressed  except  a single  reference  location  at  the 
u-nt.i  ol  iln  slab  The  over  ma>  also  print  principal  betiding  stresses  computed  from  a constant  thickness  I in  place  ol  prinupal 
niontrni  I tio  stress  tias  liii  same  sign  as  the  principal  moment.  ,s  value  of  thickness  must  then  be  available  in  Table  B2 
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TaUcK.  CONSTANTS. 


Vart.iblcs 

K 

h 

V 

F. 

» 

k 

lvpic.ll  Input  Units: 

in. 

3 

tn. 

none 

Ibs/ini 

in. 

lbs/!n' 

Only  one  card  h needed  for  Parent  or  Independent  pntbkint.  Ihh  table  it  omitted  for  Offspring  problems. 

Poisson's  ratio  will  be  taken  as  mo  aniess  specified  (always  positive). 

Data  for  computation  of  optional  stiffnesses  need  be  entered  only  if  the  options  are  exercised  in  Table  Bl . However,  they  may  be  entered 
and  thus  pinted  for  reference  if  desired.  The  slIfTnetses  arc  computed  from  these  customary  relationships; 


D*=  0)'  = 


Et» 

12(1— v*l 


C=^ 


Et» 

12(1  +v) 


S = kh,hy 


fhe  computed  stiffnesses  are  constant  over  the  X by  Y area  of  the  slab  and  appropriate  quarter  and  half-values  are  generated  at 
corners  and  edges.  They  are  printed  as  an  extension  to  Table  B2  and  thus  any  additional  stilTncss  values  entered  in  Table  B4  arc 
superimposed  algebraically. 


TaMe  B3.  SPCanED  AREAS  FOR  SELECTED  OtrmJT. 

A maximum  of  ten  cards  may  be  used,  with  a maximun  of  XJO  points  specified  by  each  card. 

The  selected  results  are  printed  if  any  of  the  four  print  options  are  exercised. 

Table  B3  may  be  omitted  if  desired  since  all  Klected  output  values  are  duplicated  in  the  complete  printout  of  the  results.  The  user  is 
cautioned,  however,  that  the  complete  printout  may  have  been  suppressed  by  the  option  in  Table  Bt. 

The  selected  output  is  controlled  by  the  same  joint  coordinate  system  shown  in  Fig  A2  and  described  below  for  Table  B4. 

TaMe  B4.  STIFFNESS  AND  LOAD  DATA. 


VarUbles; 

0* 

Df 

Q 

S 

C 

Typical  Input  Units: 

lb-in* 

lb-in* 

tb 

Jb 

lb-in* 

in. 

in. 

in. 

in/rad 

All  data  are  described  with  a coordinate  system  natation  which  is  related  to  the  discrete-element  model  of  the  slab.  This  b shown  in 
Figure  A2. 

Todisiributedata  over  a rectangular  area,  the  lower  teft-hand  and  the  upper  right-hand  coordinates  must  be  specified.  Figure  A3  illu- 
strates a sample  data  input. 


All  data  entered  In  this  (able  arc  algebraically  added  to  (he  optional  sliffneu  values  generated  from  Tables  Bt  and  B2. 

To  specify  data  at  a single  location,  (he  same  coordinates  must  be  specified  in  both  the  "From"  and  "Through"  columns. 

Ihc  "Through"  coordinates  must  always  be  equal  to  or  numcricalty  greater  than  the  "From"  coordinates. 

fhe  user  may  input  values  on  (he  edges  of  the  slab  and  the  comers  to  represent  the  proportionate  area  desired  as  illustrated  in  Figure  A 3 

I'hert  ate  no  reslrktions  on  the  order  of  cards  in  Table  B4.  Cumulative  input  is  used,  with  full  values  at  each  coordinate. 

Unit  stiffness  values  D*  and  Df  are  input  at  all  joints.  The  values  are  reduced  proportionately  for  edges. 

Load  values,  Q,  and  support  springs  S for  any  joint  an  determined  by  multiplying  the  unit  load  or  unit  support  value  by  the  appropriate 
area  of  the  real  slab  assigned  to  that  joint.  Hinged  supports  are  provided  by  using  large  S values.  Concentrated  loads  that  occur 
between  joints  can  be  proportioned  gaomctrically  to  adjacent  joints. 

Unit  twisting  stiffness  C is  dcfloed  for  the  mesh  of  the  plate  or  slab  surrounded  by  four  rigid  bars  and  four  joints.  The  mesh  is  numbered 
according  to  the  joint  number  at  the  upper  right  comer  of  the  mesh  ai  shown  In  Figure  A2. 
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TaMt  M.  AXIAI.  THRUST  DATA. 


VariabU-v 

P’ 

|.V 

1 yptcal  Input  Units 

lb 

lb 

All  dait  in  thii  labk  arc  concentrated . Diitriboted  data  must  be  Mmmed  over  the  width  of  the  incrcmcnl  involved.  Proportionair  values 
can  be  used  atong  ed|cs. 

Axial  tension  1 1 I or  compression  (— ) values  P arc  specified  for  each  a-bar  or  jr*bar.  There  is  no  mechanism  in  the  program  to  automati- 
cally distribute  the  internal  elTects  of  any  eilemally  applied  axial  loads. 

The  axial  thrust  P*  refers  to  the  force  In  the  x-bar  in  ‘he  x-direction.  Since  H is  a bar  value,  no  coordinate  should  be  uitd  which  would 

specil)  a P*  value  in  a bar  outside  the  real  plate  m slab.  The  bars  arc  numbered  according  to  the  joint  number  as  shown  in  Figure  A2 

TaMa  M.  SPCOAL  LOAD  PATTCRNS. 

A maximum  of  4 different  patterns  118  cards)  may  be  used  with  up  to  12  loads  in  each  pattern. 

The  first  load  entered  in  columns  4 through  14  Is  designated  w the  reference  load  for  placement  of  the  pattern  by  table  B7.  It  may  have 
a xcro  value  for  convenience  if  desired. 

Tfse  coordinates  entered  in  this  table  refer  to  the  relative  positiuns  of  the  toads  within  the  pattern.  Ibis  is  illustrated  by  Fig  A4  ( nordi 
nates  may  also  have  negative  values  if  desired. 

Special  Load  patterns  need  not  be  entered  if  they  ate  the  same  as  those  In  previous  probicmt.  however,  additional  patterns  may  be  used 
Any  new  patterns  entered  replace  those  with  the  same  pattern  number.  All  load  patterns  defined  in  previous  problems,  including 
previous  Parent  and  Independent  problemt.  are  retain^,  and  may  then  be  used  by  specifying  their  placement  in  fable  B7. 

TaMe  IT.  PUCEMENTS  OF  TAILE  W LOAD  PATnCRNS. 

Each  pattern  defined  by  Table  Bb  may  be  placed  at  a number  of  different  locations  on  the  slab.  Any  patterns  specified  In  previous  prob- 
lems may  be  also  placed  by  this  table:  they  need  not  be  re-entered  in  this  problem's  Table  B6.  Pattern  numben  must  be  unique 
throughout  the  scries,  however. 

Placements  ol  the  telerencc  load  arc  by  coordinates  which  are  related  to  the  overall  system  defined  by  Table  B2.  Placements  outside  the 
boundaries  of  the  actual  slab  are  pcrmissable.  Any  loads  of  the  pattern  that  thus  faN  off  the  slab  are  not  included  In  the  solution.  For 
leferencc.  however,  a summation  of  all  loads  placed  on  and  off  the  slab  is  given  in  the  output  of  Table  87. 

Ifonly  new  patterns  and  placements  are  specified,  and  no  stiffnesa  data  is  changed  In  Tables  82,  84.  or  8S.  then  this  ptublem  is  an  nil- 
spring  p^lcm. 


TdMt 

Na. 

Matfranni 

Ne.a( 

Faraalerii 

Fn« 

OApetag  FreWama 

Owla 

Raqrired 

OptiMal 

Raqtrirad  Optlmml  DMbwtd 

81 

1 

X 

X 

82 

1 

X 

X 

83 

to 

X 

X 

84 

no  limit 

X 

X* 

8S 

no  limit 

X 

X 

Bb 

IB 

X 

X 

87 

no  limit 

X 

X 

*LoadOnly 

Tbs 

user  should  also  ensure  that  if  Optional  Siiffhau  Options  are  estreised  in  Tabic 

81.  that  appropriate  data  is  available  In  Table  B2.  He  should  aleo  verify  that  if  the 
Option  to  Delete  Drisiled  Output  Is  used,  that  sufficient  output  trill  be  printed  tor 


areas  specllied  In  Table  83. 

Figsw«  Bl.  ClMeUidi  iapvA  Uhltst. 
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GLOSSARY  OF  NOTATION 


NOT At  ION  FCM  SLAI  SCC 


A(  ill 

AAI  til 
ALF 

mMK  til 
Ai42(  tl) 

ANK  It  AN2(  I 


RCCUNStON  COCFFICIENT 

COCFF  IN  STIFFNESS  HATH IX 

ANGLE  ON  HOHN'S  CIKCLE 

NECURStON  CCCFFICICNT  A(  tl)  AT  J*1 

RECURSION  COEFFICIENT  AC  ,1)  AT  J-2 

lOENTIFICATION  AND  HEHARKS  (ALPHA  - NUlM 


ATC  t ) 
NC  t I 
BBC  tS) 
BETA 
BHA 
BMO 


BMT 

BHXC  t I 
BMYI  t I 
bHlC  t I 
bH2(  t I 
Cl  t I 
CCC  tS) 

CHC  t I 
CNN 

CHIC  t I 
CM2I  t I 
CROC  I 

0(  t I 

UOC  tS) 

OXC  t I t DYC 
OXN  t OVN 
EC  I 
EEC  I 
ER 

FFC  I 
MX 

MXOHY 

HXOHYS 

HY 

HYOHX 

MYUHXS 

I 

INI 

1N2 

CNISC  I 

IN2S(  I 

lOPC 

lOPO 

lOPPS 

lOPS 

lOPS 

lOPP 

ISTA 

ITESTC  I 


TEMP  STORAGE  FOR  AC  tl)  RECURSION  COEFF 
RECURSION  COEFFICIENT 
COEFFS  IN  STIFFNESS  KATRIX 
half  theta  ICOUNTER  CLOCK-jtSE  IS  4-1 
AVERAGE  OF  X AND  Y BENDING  MOMENTS 
FIRST  principal  BENDING  MOMENT 

tHX  ~ SNA 

AOIUS  OF  MOHRi  CIRCLE 
SECOND  PRINCIPAL  BENDING  MOMENT 
•ending  moment  in  the  X DIRECTION 
•ENDING  MOMENT  IN  THE  Y DIRECTION 
RECURSION  COEFFICIENT  Bl  t I AT  J-1 
RECURSION  COEFFICIENT  Bl  t I At  J-2 
RECURSION  COEFFICIENT 
COEFFS  IN  STIFFNESS  MATRIX 
TMISTING  STIFFNESS  PER  UNIT  WIDTH 
INPUT  VALUE  OF  TWISTING  STIFFNESS 
RECURSION  COEFFICIENT  Cl  t I AT  J-1 
RECURSION  COEFFICIENT  CC  t I AT  J-2 
CROSS  bending  STIFFNESS  FOR  PR  EFFECTS 
RECURSION  MULTIPLIER 
COEFFS  IN  STIFFNESS  MATRIX 
•ending  STIFFNESSES  PER  UNIT  WIDTH 
INPUT  VALUES  OF  PENDING  STIFFNESSES 
RECURSION  MUUIPLIER 
COCFF  IN  STIFFNESS  MATRIX 
A TEST  TO  ELIMINATE  REMNANT  REACTIONS 
COEFF  IN  LOAD  VECTOR 
INCREMENT  LENGTH  IN  X DIRECTION 
HX  DIVIDED  OV  MY 
HX  DIVIDED  BY  HY  CUBED 
increment  length  IN  Y DIRECTION 
NY  DIVIOEO  BY  HX 
HY  DIVIOEO  BY  NX  CUBED 
station  number  in  X DIRECTION 
INITIAL  EXTERNAL  X COORDINATE 
FINAL  EXTERNAL  X COORDINATE 
initial  external  X COORDINATE  IN  TADLE  1 

final  external  X coordinate  in  table  3 

TWISTING  STIFFNESS I optional  I SWJTCH 
•ending  STIFFNESS  (OPTIONAL)  SWITCH 
OPTION  FOR  PRIN  STRESS  INSTEAD  OF  MOM 
SUBCRADE  MODULAS  (OPTIONAL)  SWITCH 
OPTION  TO  SUPRESS  DETAILED  OUTPUT  (TAB  I 
internal  plot  option  SWITCH 

external  X coordinate  NUMAcr 
alpmanumchic  Blanks,  used  to  terminate 

THE  PROGRAM 
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l?NOf 

19APR 

17N09 

12NOD 

12N09 

280C7 

12N09 

12N09 

03JA9 

17N09 

280E7 

09JAR 

280E7 

280E7 

03JAR 

03JAB 

03JAB 

I2N09 

12N09 

12N09 

03JAR 

05JE" 

OSJEB 

12N09 

12N09 

30JLR 

12N09 

03JAB 

2BOE7 

03JAA 

12N09 

09JAB 

28DE.7 

03JAA 

03JAA 

280F7 

28DE7 

280E7 

280E7 

28DE7 

280r7 

03JA8 

28DE7 

I2N09 

12N09 

07SEC 

O7Sf0 

07SEP 

07Sf0 

07SE0 

07SEC 

12N07 

12N07 

UNOY 


r*rtn/>r»rtftnr*Ar\AArtAnr»nr>nnnnnr»nr»Ar»r»Anr*nr>Anr»r»nr»nr»r»nnnr\onAAf%r»r»r» 


II  * 12 

J 

JJ 

JNl 

^N2 

JNISI  ) 

JN23<  ) 

J$TA 
J1  t J2 
It 

KASEXI  I 
KASEYt  ) 

ALONG 

KSHORT 

AML 

ARROBI  I 

LL 

t2 

LI 

HU 

MX 

NXPl  THRU  HXP5 
MY 

MYPl  THRU  KYP5 
H 

MCT2 

ACTS 

NOE 

NPROBI  I 
OOHX 
OOHXMY 
. UOHX2 
OOHY 
OOHY2 
POMXHY 
PKMAX 
PR 

PStGO 

PX(  • ) * PYI  t I 
PXN  t PYN 
0(  * ) 

UBMX  t OBMY 
UN 

wPX  I UPY 
wTMX  t CTMY 
REACT 
$i  t I 
$012 
$ICO 
' $AA 
$N 
SOM 
SUHR 

theta 

THA 


INI  AND  IN2  PLUS  2 
STATION  NUMBER  IN  Y DIRECTION 
J FOR  SUOROUTINE  MATRIX 
INITIAL  external  Y COORDINATE 
FINAL  external  Y COORDINATE 

initial  external  y. coordinate  in  table  1 

final  external  y coordinate  in  table  3 

external  y coordinate  number 

JNl  AND  JN2  plus  2 

00  LOOP  INDEX  USED  INSTEAD  OF  I 

OPTION  FOR  SELECTED  PRINT  OF  BMX 

OPTION  FOR  SELECTED  PRINT  OF  BMY 

MAXIMUFt  REAL  Y DIMENSION  SIZE 

MAXIMUM  REAL  X DIMENSION  SIZE 

KEEP  MULTIPLE  LOADING  FOR  ERROR  CHECKS 

PROBLEM  NUMBER  FROM  PARENT 

DO  LOOP  INDEX  FOR  REVERSED  J 

RL0NG43*  USED  FOR  VARIABLE  OIMfNSIONING 

KSHORTF'S*  used  for  variable  OIMfNSIONING 

MULTIPLE  LOADING  SWITCH 

NUMBER  OF  INCREMENTS  IN  X DIRECTION 

MX  ♦ 1 THRU  MX  ♦ 5 

NUMBER  OF  INCREMENTS  IN  Y DIRECTION 

MY  ♦ I THRU  my  ♦ 5 

INDEX  FOR  READING  CARDS 

NUMBER  OF  CARDS  IN  TABLE  2 

NUMBER  OF  CARDS  IN  TABLE  3 

NUMBER  OF  DATA  ERRORS 

PROBLEM  NUMBER  (PROG  STOPS  IF  BlANK) 

ONE  DIVIDED  BY  MX 

ONE  DIVIDED  BY  MX  TIMFS  HY 

ONE  DIVIDED  BY  HX  SQUARED 

ONE  DIVIDED  BY  HY 

ONE  DIVIDED  BY  HY  SQUARED 

POISSONS  RATIO  DIVIDED  BY  HX  TIMFS  HY 

LARGEST  PRINCIPAL  MOMENT 

POISSONS  RATIO 

PLOTTING  ARRAY  FOR  PRIN  STRESS  OR  MOMENTS 
AXIAL  TENSIONS  IN  X AND  Y DIRECTIONS 
INPUT  VALUES  OF  X AND  Y AXIAL  TFNSIONS 
TRANSVERSE  LOAD  PER  JOINT 
LOAD  ABSORBED  IN  BENDING 
INPUT  VALUE  OF  TRANSVERSE  LOAD 
LAOO  ABSORBED  DUE  TO  AXIAL  TENSIONS 
LOAD  ABSORBED  IN  TWISTING 
SUPPORT  REACTION  PFR  JOINT 

spring  support*  value  per  joint 

MOMENT  MULTIPLIER  FOR  PLATE.  STRESS 
LARGEST  PRINCIPAL  MOMENT  OR  STRESS 
SUBGRAOE  MODULUS 
INPUT  VALUE  OF  SUPPORT  SPRINGS 
VARIABLE  LABEL  (-STRESS  OR  -MOMENT  I 
SUMMATION  OF  REACTIONS  FOR  STATICS  CHECK 
MOHRS  CIRCLE  ANGLE  BETWEEN  X AND 
PRINCIPAL  MOMENT 

THICKNESS  OF  SlAB  FOR  STRESS  CALCULATIONS 


280ET 

280E7 

28DE7 

28DF7 

280E7 

12N09 

12N09 

12N09 

280E7 

03JA8 

12N09 

12N09 

01AG8 

01AG8 

280E7 

12N09 

03JA8 

19AP8 

19AP8 

03JAS 

28DE7 

01AG8 

28DE7 

OlAGR 

03JA8 

03JAR 

OlAGfl 

03JAR 

12N09 

28Dr7 

28DE7 

03JAB 

280F7 

03JAB 

28DF7 

03JAB 

28DE7 

07SE(' 

OlAGB 

OlAGR 

OlAGB 

OlAGR 

03JAB 

OlAGB 

OlAGB 

OlAGB 

OlAGB 

03JA8 

03JAB 

07SEP 

2eDF7 

07SEC 

03JAB 

17N09 

17N09 

C3JA8 
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TGKXt  TCHy 

TE^•>P  STORAGE  ARRAYS  FOR  PLOT  ROUTINE 

oysec 

XX 

IAOEPENOENT  variable  FOR  PLOT  ROUTINE 

07SF? 

••UTATION 

TSEO 

SUBROUTINE  SPLOT  3 AND  4 

07SE0 

TSEO 

t»ETA 

SCALED  DATA  1 FLT  PT  > 

07SEC 

1 

general  index 

07SCf 

IOTA 

FIXED  BETA 

07SEP 

IR 

POINT  positioner  FOR  PLOTTING  STRING 

07SFP 

IS 

INDEX  OVER  X-STA  TO  OE  PLOTTED 

07SFC 

ISKP 

length  OF  SPACE  array  TO  PF  USED 

OYSE" 

JS 

INDEX  OVER  J-STA  TO  BE  PLOTTED 

OYSFf 

L 

INDEX  OVER  ISKP 

07SEC 

OMEGA 

MAXIMUM  VALUr  IN  X ARRAY 

OYSEP 

StCMA 

scaler 

OYSEP 

SRACd 

INITIAL  PART  OF  PLOTTING  STRING 

OYSEP 

SYMUI  J 

LAST  4 characters  OF  PLOTTING  STRING 

OYSEP 

SYMO 

SYMBIIRI  - TEMP 

OYSEP 

THETA 

MINIMUM  VALUE  IN  X ARRAY 

OYSEP 

hlUTH 

NIDTM  OF  PLOT  IN  CHARACTERS  (COLUMNS) 

OYSEP 

X 

array  to  BE  PLOTTED 

OYSEP 
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rtf>r>rtnnrtr\f>rtn/^r>nrtr\nnrtr>nrtrtnr\nt^r>rtr>r>nr>r> 


USTiNQ  OF  PROQRAM  DECK 


PRO0HA.M  SLA030CJ  IKPUT.  OUTPUT.  TAPEl,  TAPf?  » 

■THIS  PKOCKAH  IS  A VtRSlOX  EXTETiOED  FROM  SLAB  30  FOR  CERl  BY 
AUSTir4  RESLaRCM  ENGIHLERS  INC*  IT  HAS  THE  ADDED  CAPABILITIES 

selected  Plot  output.  pave.'.ent  oriented  input,  and  the 

ABILITY  TU  INPUT  A SPCCIEIEO  PATTERN  OF  LOADS  AT  ANY  DESIRED 
LOCATION  . 

■FOR  DIFFERENT  SIZED  PROBLEwS  ONLY  THE  DlKCNSION  CAROS  OF  THIS 
URIVLR  NEED  oE  CHANGED.  FOR  EXAMPLE*  AAIS^^S.I)  « BBiS^43*3t  * 
ALSU  CHANGE  ASHCRT  AND  XLOi.G  TO  EUUAL  THE  REAL  X AND  Y STATIONS. 

CCIS«3.Si  t eiS>3.S*3l  * BMXIS43.L4^3>  WHERE  S AMD  L RFFER 
TO  THE  SHORT  AND  LONG  LENGTHS  OF  THE  REAL  PROBLELi, 

■THIS  PROGRAM  IS  NOW  DIMENSIONED  TO  SOLVE  A 20  BY  40  GRID. 


■THIS  PROGRAM  wiLL  OPERATE  ON  EITHER  COC6600  OR  IBM360/50  SYSTEMS. 
THOSE  CAROS  NEEDED  TO  OPERATE  UN  THE  IQM360/S0  ARE  INCLUDED  AS 
FOLLOWING  companion  CARDS  TC  ThE  CDC  CARDS  AND  HAVE  A C IN  COLUMN 
ONE  .'NU  THL  SYMBOLS  IQK  IN  COLUMNS  76  THRU  60*  OTHER  ADDITIONAL 
CAROS  SUCH  AS  the  selective  OOuDLE  PRLCISION  STATEf^EMTS  ARE  ALSO 
TAGGED  WITH  IBM  ANO  NULlCO  WITH  A C*  WHEN  CONVERTING  TO  THE 
tBM360/S0  SYSTEM.  THE  COMPANION  C0C6C00  CAROS  SHOULD  BE  RETAINED 
ANO  NULLED  WITH  AN  ADDED  C. 


DOUBLE 

PRECISION 

AA. 

6b  t 

CCt 

OOt  EEt 

FFt 

Wt 

A. 

AMlt 

AM2t 

Bt  BHlt 

6M2t 

WPlt 

Ct 

CMlt 

CM2t 

Of  Et 

ATt 

WP2t 

ALFt 

BETAt 

BMAt 

RMOt  BMP. 

BMRt 

BMT, 

CROt 

ERt 

HXi 

HXDHYt 

HX0HV3. 

HYt 

HYDHXtHVDHX3.  OOHX.OOHXHY. 

90HX2t 

POHXHY . 

Pit 

PMKAXt 

PR. 

QBHXt 

OBMY. 

REACT  t 

SDT2t 

SICO. 

SUMRt 

THETA. 

THKt 

THAt 

TKXt 

TMY.WSUMl. 

WSUM2t 

WSUM3. 

OPXt 

OPYt 

OTMX,  OTMY 

dimension 

AAI  29 

t 1 1. 

EEt  23 

t 1 I. 

FF( 

23 

1 1 

It 

Al  29 

t 1 It 

AMlI  23 

t 1 I. 

AM2( 

23 

t 1 

)t 

WPK  29 

t I It 

WP2I  23 

t 1 i • 

BB( 

29 

t 3 

It 

ODI  29 

t 3 It 

CC  ( 

29 

t S 

)t 

Bl 

29 

t 23 

It 

BHII  29 

t 23  It 

BM2I  23 

t 23  It 

Cl 

29 

f 29 

It 

CMII  29 

t 23  It 

CH2(  23 

t 29  1. 

01 

29 

t 29 

It 

Cl  29 

t 29  It 

BMX  I 

29 

t 43 

It 

SMYI  23 

1 49  It 

01  23 

• 43  I . 

SI 

29 

t 43 

It 

CHI  29 

t 49  It 

0X1 

29 

. 49 

It 

OVI  29 

t 49  It 

PXI  23 

. 43  I. 

B 

PY( 

23 

t 43 

It 

Ml  23 

t 49  1. 

ATI  29 

. 43  I 

DIMENSION  PSIGO  (23.43) 
XSHORT  . 20 
KLONG  r AO 
LI  f XSHORT  4 3 
L2  • KLONG  4 3 


CALL  SB30ES 

1 

2 

9 

END 


AA*  BB*  CC|  00. 
AH2t  Bt  BHl.  BM2t 
NPl.  WP2«  BMX.  BHY. 
St  PX.  PY.  W. 


EE»  FF'.  A* 
Ct  CMl.  CM2. 
CM.  PSIGO. OX. 
LI.  L2 


AT.  AMI. 

0.  E. 

OY.  0. 


26JF0 


RE-OIMFN 


26SE9IP'- 

26SF6IFM 

26SE9IP*' 

265E9I«»-' 

26SE9tPM 

26SE9IPM 

26SE9IPV 

26SC3IPr.’ 

26SE9IP^* 

26SE9I»'^' 

re-oi-mfn 

rf-dimfn 

RF-OP'PN 

RE-DIMFN 

RE-DIVEN 

RE-01Mr*J 

RE-OIMEN 

RE-or’FN 

RE-OIMfn 

RE-DIKEN 

RE-D1^*FN 

RE-OIMF.i 

RE-OI*' 

Rp-ni v^^ 

RE-nr'rv 

OS  JAB 

0SJA8 

26JEO 

OSJAB 

23JF0 

19APA 

OSJAB 


SUORwUTIKE 

SaSCES 

1 AA, 

SB,  CC* 

DO, 

EE» 

FF,  A, 

AT,  AMI, 

26Jrp 

1 

AM2, 

e,  SMI, 

l)M2* 

C» 

CM  I,  CM2, 

D,  E, 

rsJAfl 

2 

WPl, 

WP2,  BMX, 

BMY, 

CH, 

PSIGO,OX, 

OY,  0, 

23JFO 

3 

S* 

Px,  PV» 

M, 

LI, 

L2 

1 

19APA 

1 FuHhnT  1 

S2H 

PROGRAM  SLAS30E 

CERL 

SPECIAL  DECX-i 

ARE-JJP,FLE 

26JE0 

20H  REVISION  DATE  07  SEP  70 


IREVISED 


C 

OOUbLE 

PRECISION 

AA, 

B6*  CC, 

00, 

EC,  FF, 

M 

26SE9to/ 

C 

1 

A , 

AMI,  AM2, 

B, 

BMl,  BM2, 

WPl 

26SE9I"V 

c 

2 

C, 

CMl,  CH2* 

0, 

E,  AT. 

WP2 

26SE91PV 

c 

3 

ALF, 

BETA,  6MA, 

BKO. 

BHPv  BMR, 

BMT 

26SC9ln" 

c 

A 

CRD, 

ER,  HX, 

MXOHY 

26SE9|i>-^ 

c 

S 

HX0HY3, 

HY,  HY0HX,HY0HX3,  ODHX.ODHXHY, 

0DHX2 

26SE9ref/ 

c 

4 

POHXHY, 

PI,  PMMAX, 

PR, 

OBMX 

?6SE9I'»'' 

c 

7 

OBMY, 

REACT,  S0T2* 

SIGO, 

SUMR, 

THETA 

26SE9IPM 

c 

i( 

THX, 

THA,  TMX* 

TMY,WSUH1,  WSUM2, 

WSUM3 

26SE9IPN' 

c 

9 

OPX,  QPY, 

OTHX,  OTMY 

26SE91PV 

DIMENSION 

' 

AAI  LI 

♦ 1 ), 

EEI  LI 

, 1 

>. 

OSJAR 

I 

FF< 

LI 

, 1 

I,  Al  LI 

,1  >* 

AHII  LI 

. 1 

>, 

OSJAA 

2 

AH2( 

LI 

, I 

i,  WPlI  LI 

* 1 ), 

WP2I  H 

, I 

) • 

C5JA9 

3 

eei 

LI 

* 3 

1,  001  LI 

• 3 >, 

OSJAS 

A 

CC  ( 

LI 

* S 

I, 

osjaa 

S 

ei 

LI 

*L1 

I,  BMII  LI 

,L1  1, 

BM2I  LI 

,L1 

) , 

I9APR 

A 

Cl 

LI 

,L1 

1,  CMII  LI 

,tl  ), 

CM2I  LI 

, Ll 

) , 

19AP». 

7 

■ 01 

LI, 

♦ .LI 

1,  El  LI 

, LI  ), 

19APA 

I 

BMXI 

LI 

♦ L2 

),  BMYI  LI 

♦ L2  ), 

01  LI 

, L2 

), 

19APR 

9 

SI 

Ll 

♦ L2 

),  CHI  LI 

, L2  ), 

07MYR 

A 

DXI 

LI 

, L2 

) , OYI  LI 

, L2  )* 

PXI  Ll 

, L2 

) , 

19AP8 

• 

PYI 

LI 

* L2 

>,  Ml  LI 

, L2  >, 

ATI  Ll 

, L2 

> 

3 IOC  9 

dimension 

IMISI  10 

>• 

JN13I  10 

1 , 

03NO9 

1 

IN23I 

10 

1,  JN23I  10 

I, 

KASEXI  10 

) , 

03N09 

2 

AASEYI 

10 

)•  NPROBI  2 

), 

XPR06 1 2 

), 

03N09 

3 

ITESTI 

2 

),  ANII  AO 

), 

AN2I  18 

>, 

03N09 

A,  . 

CROI 

S 

» 

03N09 

OIKEi,SION  PS  ICO 

1 LI 

, L2  t 

•28JE0FN 

COt-IMuN 
« FORMAT 

10  FORMAT 

11  FORMAT 

12  FORMAT 
IS  FORMAT 
lA  FORMhT 
IS  FORMAT 

IS  format 

If  FOHM»,T 
1 

2 / 5X 

/ SX 
/ SX 
/ SX 


//IfX  lOHOY 


SPLT  / ,MXt  my,  TH<,  KCT3»  NOC,  PR,  SOM 
Cups  / EM,  SKK,  lOPO,  lOPC,'  I0PS»  Hx,.  my 

» 

SH  ,''IOX,  iOHI TRIM  » 

SHI  , SOX,  lOMi TRIM  > 

2CAA  ) 

SX,  20AA  I 

Al,  AA,  SX,  17AA,  A2  » 

(///lOH  PROS  , /SX,  Al,  AA,  SX,  17AA,  A2  ) 

(///17H  prod  ICONTO),  /SX,  Al,  AA,  SXt  17AA,  A2  > 
t AI/ltAAH  •••THIS  SERIES  CF  PROBtEKS  WAS  RuM  USING  A 
AOH  PROGRAM  OEVELOPEO  FOR  THE 
S2MTRANSP0RTAT10:.  FACItITIES  BRANCH, 

52H0EPARTMtNT  OF  THE  ARMY, 

S2HC0NSTRUCT1UN  ENGINEERING  RESEARCH  lAQORATORY 
S2HCHAKPAIGN,  ILLINOIS 


// 


/ SX  2AHAUSTIN  RESEARCH  CNCINEERSi 
/ SX  20H  AUSTIN,  TEXAS  > 


31  AGO 
17AG0 
AMY3 
03FEA 
03FEA 
170C9 
170C9 
170C9 
170C9 
170C9 
18AG0 
31  AGO 
31  AGO 
31A60 
31  AGO 
SIAGO 
07SE0 
31  AGO 
31AG0 


SU 


4*  ruKMAT  iSX*  2(lX«l2«lX*t3l»22X.Cll*3> 

4*  (SXt  2(lX.t2«lx*t3t*22X«lPC11.3  I 

47  FUKMAT  ( 3X.  12.  1X«  13.  6E11.3.  F6*l 

47  format  ( 5X.  12.  IX.  13.  IP6EU.3*  0PF6.1  > 


2SOC7COC 

ITOCtl*'; 

izoetcdc 

170C*H»-* 


4t  format  I 
1 
2 

3 

4 
3 


X MOVENT  AND  X TWISTING  MOMENT  ACT  IN04N09 


32H  TME  X DIRECTION  (ABOUT  Y AXIS). 


04NO9 


Y twisting  moment  • -X  TWISTING  MOM  04M09 


3SHENT*  counterclockwise  BETA  ANGLES  •/• 


04N09 


ARE  POSITIVE  FROM  X AXIS  TO  THE  OlR  04NO9 


28HECT10N  OF  LARGEST  PRINCIPAL 


18AG0 


49  format  (///50H 


•**  THE  detailed  output  HAS  BEEN  DELETED  BY  T 18AC0 


21HHE  OPTION  IN  TABLE  1. 


18  AGO 


2 /»  9X 

3 


SOHA  SINGLE  SET  OF  VALUES  IS  PRINTED  AT  OR  NEAR  THE  S 07SEO 


29HLAa  CENTER  FOR  REFERENCE 


50  FORmmT  (///SOH 


55  FORMmT  « //44H 


SOH  STATICS  CHECK*  SUMMATION  OF  Rl 

«HS  ■ • E10.3  ) 

6HS  > * 1PE10.3  ) 

I.4H  table  3.  SPECIFIED  AREAS  FOR  SELECTED 


//  ) - 07SE0 

summation  of  reaction.  120E7 

25JA8COC 


1 I SHI 

2 ASH 

3 / 4SH 

A IX  I8H 

SO  format  < //3SH 

53  format  ( //40H 

54  format  I //47H 

1 32H 

7S  FORMAT  < //SOH 
91  format  (///30H 
1'  20H 


ISHPLOTTEO  OUTPUT  ♦// 

4SH  PLOT  II-YESI 

ASH  FROM  THRU  OEFL  X-MOMENT 

18H  Y-MOMENT  PR IN  * AS  > 

3SH  table  S.  special  LOAD  PATTERNS//) 

40H  TABLE  9.  SELECTED  OUTPUT 

47H  TABLE  9.  SELECTED  OUTPUT  — USING 

32HS7IFFNESS  DATA  FROM  PARENT  PROS  . Alt  A4  ) 


170C9IPM 

04N09 

2SJE0 

76JEP 

17A60 

17AG0 

2SUE0 

2SUF0 

2SJE0 

07SE0 


TABLE  7.  PLACEK3NTS  OF  SPECIAL  LOAD  PATTERNS  /)I7AG0 


91  format  I //SIH 


30H  ••••  PROBLEM  terminated- 

20H  DATA  ERRORS  •••*  ) 


I80E7 

I80E7 


CAUTION*  MULTIPLE  LOADING  OPTION  MISUSED  I80E7 


3SH  FOR  THIS  OR  PRIOR  PROBLEM  ***• 


93  format  I///38H 


••••  PROBLEM  WILL  BE  TERMINATED. 


1 40HTHE  DIMENSION  STORAGE  IS  TOO  SMALL  ) 

94  format  I (///47H  •••  NO  CAROS  ALLOWED  IN  TAD  2 OR  S FOR  ML»-l 

95  FoRM«T  I ///SOH  •••  PRIOR  PROBLEM  IS  NOT  A PARENT  *** 

9S  format  ( ///SCH  •••  My  should  not  be  less  than  MX 


1I0E7 
290C9 
290C9 
)17AC0 
) I TACO 
)I7AG0 


97  format  ( S^/XSOH  •••ILLEGAL  DATA  INPUT  ATTEMPTED  IN  OFFSPRING  PR0B)31AG0 


91  F0R.3mT  I///40H 


••••  UNOESIGNATED  ERROR  STOP  ••*•  ) 


-PROGRAM  AND  PROBLEM  IDENTIFICATION 


ITEST(l)  • IH 
ITESTdl  • 4H 
XML  « 0 

READ  12*  I ANl(N)*  N ■ l«  40  ) 

CALL  TIC  TOC  (1) 

1010  READ  14*  NPROB.  ( AN2(N).  N • 1*  18  > 

IF  « NPROOdI  - ITESTIl)  I 1020.  lOlS.  1020  • 
lOlS  IF  I NPR0BI2I  - ITCST(2)  ) 1020.  9990.  1020 

1020  PRINT  11 

1021  PRINT  1 

PRINT  13*  ( ANIIN).  N > 1.  40  ) 

PRINT  IS.  NPROB*  ( AN2(N)*  N • 1*  18  ) 


REWIND  1 
REWIND  2 


NOE  • 0 


180E7 


I70C9 

170C9 

1S0E7 

170C9 

2SSFS 

170C9 

020C9 

020C9 

2SAG3 

)9MRS 

170C9 

170C9 

200C7 

20OC7 

180C7 
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r%r\n  r\r\r>  r\r>r> 


■ INWTTAiif  1 

R£AO  2w.  .<tLi  NCT2t  KCT).  ««CT4*  NCT»«  NCT6*  NCTT*  lOROt  lORCt  02JL0 

1 lOPS*  lUPUt  IOPP5  17AGC 

PRIoiT  KLi  NCT2*  i\CT3.  hCTA*  NCT).  NCT6*  NCTTt  lOPDt  tOPCt  02JL0 

1 lOPSt  lOPttt  lOPPS  17AGO 

IF  I MU  i 62f  6U  A1  C70C9 

INPUT  TABLE  2 

A1  READ  21.  MX.  MY.  HX.HV.  PR.  EM.  THK.  SKK  17AG0 

«2  PRINT  22.  MX.  MY.  HX.HY.  PR.  EH.  THK.  SKK  17AG0 

IF  I HL.NE.-l  I CO  TO  69  17AGO 

IF  I NCT2  ♦ NCT»  .NE.  0 I PRINT  94  ''  17AG0 

IF  ( NCT2  4 NCT6  .NE.  0 I NUE  • NOE  4 1 17AG0 

69  IF  I KHL  ) 79.  7C.  72  17AGO 

TO  IF  C ML  » 71.  79.  79  19DE7 

71  NOE  > NOE  4 1 ieOE7 

PRINT  99  17AG'' 

CO  TO  79  19DE7 

72  IF  < ML  » 75.  73.  73  15DE7 

73  PRINT  92  190E7 

75  KML  • ML  190E7 

IF  I L1-3-MX  > 76.  77.  77  290C9 

76  PRINT  93  - 290C9 

NOE  « NOE  ♦ 1 • 290C9 

77  IF  I L2-3-MY  I 76.  79.  79  290C9 

76  PRINT  93  290C9 

NOE  • NOE  4 1 290C9 

79"  IF  ( MX  - MY  I •!»  «1»  10  290C9 

•10  NOE  ■ NOE  4 1 180E7 

PRINT  96  17AG0 

61  IF  I HX  • HY  I 62.  62*  83  - 29N07 

82  NOE  ■ NOE  4 1 180E7 

•3  IF  I PR  4 THK  ) 84.  89.  89  29N07 

84  NOE  ■ NOE  4 1 18DE7 

89  CONTINUE  29N07 

COMPUTE  FOR  CONVENIENCE 

I ML  I 136.  100.  ICO  . 29N07 

100  MXPl  • MX  4 1 29M07 

MYPl  » KY  4 1 29N07 

MXP2  ■ -tX  4 2 060E7 

MYP2  • MY  4 2 29N07 

MXP3  • HX  4 3 080E7 

MYP3  ■ MY  4 3 080E7 

MX?4  '■  fix  4 4 080F7 

MYP4  • my  4 4 080E7 

MXP9  ■ MX  4 9 080F7 

MYP5  » MY  4 9 080E7 

OOHX  • 1.0  / HX  14SF6 

OOHY  • 1.0  / MV  14SE6 

OCHXHY  > OOHX  • OOHY  14N07 

POHXHY  » PR  * OOHXHY  14N07 
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r\r\  Ck  non 


00HY2  « OOHY  ’ OOmY 

UN07 

00HX2  « OCHX  • OCXtX 

lANOT 

MXOMY  « HX  • OOHY 

14N07 

HYOHX  « HY  * OOHX 

UNO? 

HY0HX3  ■ HYOHX  • 00HX2 

1AN07 

HX0HY3  •^XDHY  • 00HY2 

14N07 

KPRuaat  ■ r<pRoa(i) 

170C9 

RPR0B(2I  * NPR0512> 

I70C9 

00  105  J ■ 1.  MYP3 

29N07 

00  103  I « It  KXP3 

29N07 

OXdiJl  • 0.0 

01N06 

OYI 1 tJI  ■ 0.0 

OINOG 

0(1 tJ)  ■ 0.0 

01N06 

S(ltJ>  > 0.0 

01N06 

CHtt  iJI  • 0.0 

07MyS 

pxdij)  ■ 0.0 

23JE7 

PYdtJJ  • 0.0 

23JF7 

BHX(ttJ)  > 0.0 

29N07 

BMY(ltJ)  • 0.0 

29N07 

tu3 

CONTINUE 

29N07 

105 

CONTINUE 

23JE7 

00  135  X « It  MXP3 

29N07 

00  130  1 « It  MXP3 

29M07 

• 

BdtXI  « 0.0 

29N07 

OHUUXt  > C.O 

2lf)F7 

CdtXI  ■ 0.0 

29N07 

CMldiX)  • OtO 

210E7 

150 

CONTINUE 

1SAC7 

155 

CONTINUE 

1SA67 

15« 

00  138  K « It  HXP3 

29N07 

A(Ktl)  > 0.0 

210E7 

AHKKtll  ■ 0.0 

29N07 

W(XtMYP3>  * O.C 

2I0E7 

WPKXtll  • C.O 

29N07 

WP2(Ktli  • C.O 

280C7 

15S 

CONTINUE 

29N07 

— WPUTTAOLea  

SOH  . 6HSTRESS 

17AG0 

IF  ( lOPPS.NE.l  ) SOH  • SMMOHENT 

17  AGO 

PRINT  55  . SOM 

17AG0 

'IP  1 NCT3  1 9990.  ISCt  150 

24SE9 

150 

IPOP  • 0 

2AJE0 

CAJ-L  sLPLT  ( M.  BHX.  BHY.  PSIGOt  LI.  l.2>  IPOPt  0 » 

2SJC0 

GO  TO  1«1 

24SE9 

110 

PRINT  3C 

24SE9 

Itl 

CONTINUE 

24SE9 

■ INPUT  TABLE  4 

190 

IF  ( ML  ) 200.  320.  320 

190F7 

200 

PRINT  34.  XPROB 

150E7 

DC  22u  J « 1.  MYP3 

190F7 

00  21^^  I > 1.  HXP3 

29N07 

Wd.JI  • 0.0 

240C7 

53 


U(  I tJI 

a O.C 

240C7 

210 

cuNi  n.uE, 

240C7 

220 

CONI INUL 

240C  1 

IF  i NCHi  1 

9983i  22<* 

. 223 

07Sf  0 

224> 

PRINT  3H 

. 

7SFP 

CO  TO  362 

07sro 

22S 

00  26v,  N 

« NCT4 

n2SE0 

READ  26  t INI. 

JNli  |N2« 

JN2*  OH 

20OC7 

PRINT  «6t  INlt 

JN1«  |M2i 

JN2*  ON 

20OC7 

11  • 

INI  4 2 

29N07 

J1  > 

JNl  > 2 

29N07 

12  « 

IN2  * 2 

29H07 

J2  « 

JN2  * 2 

29N07 

IF  ( INI  - 

IN2  1 232< 

232.  231 

29N07 

231 

NOE  * 

NOE  ♦ i 

180E7 

232 

IF  < JNl  - 

JN2  1 234 t 

234.  233 

29N07 

233 

NOE  ■ 

NOE  ♦ 1 

180F7 

234 

IF  I IN2  - 

MX  I 236* 

236.  2J3 

29N07 

233 

NOE  ■ 

NOE  * 1 

18DF7 

234 

IF  ( JN2  - 

MY  > 238* 

238.  237 

29N07 

237 

NOE  > 

NOE  + 1 

ieDF‘» 

231 

00  233  1 

< 11*  12 

29Nn7 

00  230  0 • 

• Jl*  J2 

20OC7 

0( ItJI 

1 • 0< 1 tJ) 

4 ON 

20OC7 

230 

CONTINUE 

200C7 

233 

CONTINUE 

2C0C7 

260 

CONTINUE 

• 

200C? 

CO  TO  386 

19DF7 

320 

PRINT  33 

13DE7 

IF 

( NCT4  ) 99BC.  323.  323 

025F'' 

323 

PRINT  38 

2SE0 

CO 

TO  362 

02SF0 

323 

00 

360  N « 1.  NCT4 

025F0 

READ  23.  IM.  JNl.  IN2.  JN2  • 

OXN. 

OYN.  ON.  SN. 

CMN 

OTMYB 

PRINT  63*  INI.  JNl.  1N2.  JN2. 

OXN. 

OYN.  ON.  SN, 

CHN 

07MY8 

IF 

( CHN  1 331*  333. 

331 

07MY8 

331 

IF 

1 INI  • |N2  • JNl  • JN2 

1 9980*  333*  333 

180E7 

333 

NOE  .NOE  4 1 

1B0E7 

333 

11  ■ INI  4 2 

150E7 

Jl  • JNl  4 2 

29H07 

12  • IN2  4 2 

29N07 

J2  • JN2  4 2 

6 

29N07 

JP 

( INl  - 1N2  1 342*  342. 

361 

29N07 

361 

NOE  - NOE  4 1 

180'^7 

362 

IF 

1 JNl  - JN2  ) 346*  366* 

363 

29N07 

363 

NOE  - NOE  4 1 

180E7 

364 

IF 

1 IN2  - MX  1 346.  346. 

345 

29N07 

363 

NOE  * NOE  4 1 

laor? 

366 

IF 

1 Ji'i2  - MY  J 348.  348. 

347 

29N07 

347 

NOE  . NOE  4 1 

• 

ieoF7 

368 

DO 

333  1 > 11.  12 

29N07 

00 

330  J > Jl.  J2 

01067 

PXn.JI  • OXd.JI  4 OXN 

6SE64 

OYII  .J)  • DYI I .J)  4 OYN 

6SE64 

0(1. Jl  • 0(1. Jl  4 

ON 

• 

13  API 

6(1. Jl  ■ Cd.JI  4 

SN 

13AP3 

S4 


n nr>  n n n 


CMd.J!  - CHtl.JJ  ♦ CHN  07MY« 

yt>0  CONTINUE  11AP3 

355  CONTINUE  22JE7 

360  CONTINUE  13AP3 

C GENERATE  EXTENCE'J  TABLE  A 

362  IF  « lOP0*l0PC*40PS  .£0.  0 > GO  TO  363  26JF0 

JF  ( ML.EO.-l  1 PRINT  97  31AG0 

IF  I ML.EQ.-l  > GO  TO  363  31AGC 

CALL  COPS  ( CH,  OX.  UY,  S.  LI,  L2  J 26JE0 

363  CONTINUE  26JE0 


INPUT  TABLE  5 


PRINT  37 

IF  I NCT5  ) 9980,  385.  36<, 
36A  00  362  N > 1,  NCT5 


READ  2A,  INI,  JNl.  IN2,  JN2,  PxN, 
PRINT  AA,  INI.  JNl.  tN2,  JN2,  PXN, 


IF 

( PXN 

1 

365.  367,  365 

365 

IF 

( INI 

1M2  1 9980,  366,  367 

366 

NOE 

NOE  ♦ 1 

367 

IF 

( PYN 

) 

366.  370 368 

368 

IF 

1 JNl 

JN2  ) 9980,  369,  370 

369 

NOE 

NOE  ♦ 1 

370 

11 

INI  ♦ 2 

J1 

JNl  <•  2 

12 

IN2  ♦ 2 

* 

J2 

JN7.  4 2 

IF 

( INI 

- 

1N2  I 372,  372,  371 

371  ^ 

NOE 

NOE  ♦ 1 

372 

IF 

( JNl 

- 

J'/I2  » 37a,  37A,  373 

373 

NOE 

• 

NOE  4 1 

37A 

IF 

1 IN2 

- 

MX  1 376,  376,  375 

375 

NOE 

n 

NOE  4 1 

376 

IF 

I JN2 

- 

MY  ) 378,  378.  377 

377 

NOE 

• 

NOE  4 1 

378 

DO 

38w 

I 

> 11,  12 

DO 

379 

J 

> Jl,  J2 

PXII 

9 

J)  » PXI I »J)  4 PXN 

PYII 

• 

Jl  ■ PYI I ,JI  4 PYN 

379 

CONTINUE 

380 

CONTINUE 

382 

CONTINUE 

00 

TO  386 

385  PRINT  38 
INr>UTTABLES8AN0  7 


22AP5 

2ASE9 

2ASE9 

PYN  110E7 

PYN  UDE7 

18PF7 

180E7 

18DE7 

180E7 

180E7 

180E7 

180E7 

29N07 

29N07 

29N07 

29N07 

18DE7 

29N07 

180E7 

29N07 

18DF7 

29N07 

18DE7 

I20F7 

125JE7 

20AP5 

20AP5 

29N07 

120C7 

120E7 

I90EV 

150E7 


386  PRINT  6v  26JE0 

IF  I NCT7  1 9980,  369,  388  31AG0 

388  CALL  PLOP  I C,  LI,  L2 , NCT6,  NCT7,  MX.  MY  I 26JE0 

- GO  TO  AOO  26JEC 

389  PRINT  38  17AG0 

IF  I NCT7  » 9980,  392,  AOO  17AG0 

392  PRINT  7A  I7AG0 

PRINT  3tt  I7AGC 


r\nc\of\r>ryr\r>  r>r» 


I 


1 


400  IF  t NOE  ) 998C.  ^lOl.  397 
397  PRINT  91 1 NOE 
GO  TO  9990 


ZfcJEO 

26JtO 

29N07 


-FORM  SUo-MATRlCcS 

-A  SPRING  IS  PLACEO  AT  PTS  uEYONO  BCUKOARIES  OF  THE  REAL  SLAB 
TO  MAF,E  solution  OF  NON-RE CT ANGULAR  SLABS  OR  SLABS  WIT 
HOLES  POSSIbLE.  THIS  IS  DONE  BY  TESTING  ON  THE  CClItS) 
TEKMSi  ANO  IF  ZERO,  SET  EOUAL  TO  1.0 

-Qd.Jl  IS  the  input  load  for  this  problem,  other  prior  problem 
LOADS  are  OISCAROEO 

00  600  J = 1 . HYP3 

00  LOA  I » 1,  MXP3 
FFtl.l)  ■ UII.J) 

IF  (J-1)  9980,  A02,  A03 
AAII  ,11  « 0.0 
GO  TO  40A 

AAtl  ,11  OYI  I ,J-11  • HX0HY3 
, CONTINUE 

IF  ( ML  1 501,  A05,  A05 
, DO  500  1 « 1»  WXP3 

IF  (J-11  9900,  AlC,  A07 
r IF  (MYP3-J1  9980,  A38,  A20 

—COEFFICIENTS  COMPUTED  AT  J • 1 


IF  (l-ll  9960,  All,  A12 
BBII  ,21  • 0.0 

BBIl ,3)  • 0*^ 

ceil, 3!  • 1,0 

CCtl.Al  « 0.0 

CC( I ,51  * 0.0 

0011 ,21  • 0.0 

00ll,31  • 0.0 

EEii,;:  “ 0.0 

GO  TO  500 

IF  IMXP3-I1  9980,  A19,  A13 
CR0(5I  ■ AMINI  ( 0X1 1 

60(1,11  • 0.0 

BC(1,2I  ■ 0.0 

B0(l,31  • 0,0 

CC(I  ,21  ' 0.0 


( CC(1,31 
CC(1  ,31 
CC(1,A1 
00(1,11 
00(1,21 


00(1,31 
EE(t  ,11 

IF  (I~2*  9980, 

ceil,  11 


INI  ( OX(l  ,J^l>  , OY(l  »J+1 

« 0.0 
■ 0.0 
• 0.0 

B 0.0 

• HX0HY3  • OY(I,J+ll 

OOHY  • PYI I ,J+l  I 
A15,  AlA,  415 

* 1.0 
. 0.0 

• POHXHY  * CR0(5) 

. -2.0  • ( HX0HY3  • OY(I,J+H 

4 POHXHY  • CRDt5>  I 
- OOHY  • PY(I,J4ll 
■ POHXHY  • CR0(51 

• HX0MY3  * OYII.J+ll 

17,  416 

> 0.0 


29N07 
08DE7 
300C9 
OBOE  7 
03JAR 
08DE7 
03JA8 
oeoE7 
080E7 
080E7 
130E7 
05N09 

080E7 

OaOE7 

08DE7 

170C9 

080E7 

080E7 

080E7 

080E7 

03JA8 

0e0E7 

130F7 

29JL0 

29JLB 

0flDE7 

08DE7 

080F.7 

080E7 

11DE7 

14DE7 

170C9 

130F7 

29JLe 

08DE7 

29JLB 

080E7 

29JL8 

03vlA8 

130E7 

130F.7 


i 

i 

i 
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417 

IF 

(HXP2-I>  99a0t  500«  4l« 

130f7 

4ia 

CC(Ii5l  • 0.0 

130F7 

CO 

TO  500 

OtOET 

419 

ao(i«i)  - 0.0 

130f7 

BB(1.2)  « 0.0 

oeoF7 

ccd.ii  ■ 0.0 

080E7 

CC(Ii3)  ■ 1.0 

170C9 

ecu  t2)  ■ 0.0 

080E7 

00(1  .21  ■ 0.0 

oeoE7 

00(1 .1)  • 0.0 

080E7 

ECdtl)  « 0.0 

Co 

TO  500 

080E7 

COEFFICIENTS  COMPUTED  FROM  J • 2 TO  MVP2 

42u 

IF 

(l>lt  9980.  421.  424 

05N09 

421 

CRD(3)  ■ AKINI  ( DXd4-l.J  ) . 0Y(l4^1.J  ) ) 

29JL8 

8B(1.2)  > 0.0 

29JL9 

BOd.3)  - POHXHY  • CRD(3) 

29JL8 

CCdt3)  ••  HY0HX3  • OX(Ul.J) 

080E7 

1 

♦ OOHX  * PXd-»l.J) 

08DE7 

IF 

( CCd  «3>  1 441.  422.  423 

05N09 

422 

CC(I«3>  ■ 1.0 

05N09 

423 

CCd. 4)  « -*2.0  • (HY0HX3  • DX(Ul.J) 

05N09 

1 

♦ POHXHY  * CR0(3>  » 

29JL8 

2 

- ODHX  • PXd  + l.JI 

oeoET 

CCd. 5)  • HY0HX3  • OXd  + l.JI 

110E7 

00(1.2)  * 0.0 

08DE7 

00(1.3)  ■ POHXHY  • CR0(3> 

29JL8 

££d.l>  >0.0 

03JAS 

. CO 

TO  500 

08DE7 

424 

IF 

(MXP3-I)  998C.  431.  425 

09N09 

425 

CROd)  • AhlNl  ( OXd-l.J  ) . OYd-l.J  ) ) 

05NO9 

CRD(2)  ■ AMINI  ( OXd  .0  ) » OYd  >J  > ) 

29JL8 

CR0(3)  « AHINI  ( 0X(U1.J  ) t 0Y(I41.J  ) ) 

29JL8 

CR0(4)  ■ AMINI  ( OXd  .J-1)  . OYd  .J-D  ) 

29JL8 

CR0(5)  • AMINI  ( OXd  .J4l>  • OYd  .J-f))  ) 

29JLR 

BBd.l)  • POHXHY  • ( CROd)  ♦ CR0(4)  ) ♦ 

30JL8 

1 

OOHXHY  • ( CHd.J)  ♦ CHd.J)  ) 

07MY8 

BOd.2)  ■ -2.0  • ( POHXHY  • ( CR0(2)  + CRD(4)  ) 

29JL8 

1 

♦ HX0HY3  * ( OYd. J-1)  4 DY(I.J)  ) ) 

14SF8 

2 

♦ OOHXHY  • ( - CHd.J)  - CH(1  + 1.J) 

07MY8 

3 

- CHd.J)  - CHd-d.J)  ) - OOHY  * PYd.J) 

OTMYfl 

BOdO)  • POHXHY  • ( CRD(3)  + CR0(4)  ) 

29JLA 

1 

♦ OOHXHY  • « CH(I*l»J)  + CH(I+1,J)  ) 

07MY8 

CCd. 2)  « -2.0  * ( HY0HX3  * ( DX(I-l»J)  ♦ OXd.J) 

) 29N07 

\ 

♦ POHXHY  • ( CROd)  ♦ CRD(2)  ) ) 

29JLB 

2 

♦ OOHXHY  • ( - CHd.J)  - CH(I.J41) 

07MY8 

3 

- CHd.J)  - CHd,J+l)  ) - ODHX  • PXd.J) 

07MY8 

CCd. 3)  ■ HY0HX3  « ( OX(l-l,J)  ♦ 4.0  • OXd.J) 

29N07 

1 

4 OXd  + l.J)  ) ♦ HX0HY3  * ( DYd.J-l)  ♦ 4 

.014SE6 

2 

• OYd.J)  ♦ OYd.J*!)  ) ♦ POHXHY  * 4.n 

14  5E6 

3 

• ( CR0(2)  ♦ CR0(2)  ) ♦ OOHXHY 

79JLB 

4 

* i CHd.J)  ♦ CH(1.J  + 1).-*  CHd-^l.cA) 

09MYB 

5 

4 CHd4l.J4l)  4 CHd.J)  4 CHd.l.J) 

07MY9 

i 

4 CH(I.J41)  4 CH(141,J41)  ) 4 OOHX 

07MY8 

7 

• ( PXd.J)  4 Pxd41,J)  ) 4 OOHY 

14SF6 

• 

• ( PYd.J^)  4 PYd.J4l)  ) 4 Sd.J) 

14SE8 

S7 


IF  ( CC(1 t3l  ) 

427*  426*  427 

09N09 

426 

CC(I«3) 

- 1.0 

05N09 

427 

CC<I t4* 

■ -2.0  • ( HYUHX3  • ( OX(l*J) 

4 0X(l4l,J) 

) 09N09 

1 

4 POHXHY  • ( CRD(3)  ♦ CR0(2)  ) ) 

29JL9 

2 

♦ 00((XHY  • ( - CH(1*1*J)  - 

CH(U1  ,J41 ) 

07MYR 

i 

- CH(l4l*J)  - CH(U1*J+1)  ) 

- ODHX 

OT^Yfi 

4 

• PX(U1*J) 

14SF6 

00(1*11 

• POHXHY  • ( CRO(l)  ♦ CR0(9) 

) 

29JLX 

1 

•»  OOHXHY  • ( CH(I*J+1)  + CH(1*J*1)  ) 

09MY4 

00(1 *21 

• -2.0  • ( HX0HY3  • ( DNIUJ)  ♦ 

DY(l*J4l)  ) 

29H07 

1 

•*  POHXHY  * ( CR0(2)  ♦ CRD(9)  ) ) 

29JLB 

2 

4 OOHXHY  • ( - CH( I *J41 ) - 

CH(i4i,j4n 

07MY8 

3 

- CH(t*J4l)  - CH(1+1,J41)  ) 

- OnnY 

07MYB 

4 

* PY(1,J41> 

14SE6 

00(1*3) 

« POHXHY  • ( CROO)  4 CRD(9) 

) 

29JLB 

1 

4 OOHXHY  * ( CH(l4l,j4l)  4 

CH(I+l,J4l)  ) 

C7MYfl 

EE(1  *11 

■ HXCHY3  • 0Y(l*J4l) 

03JAR 

IF  (1-2)  9930* 

429*  426 

09N09 

428 

CC(1*1) 

• 0X(I-1*J)  * HY0MX3 

0$N09 

429 

IF  (MxP2  - 1 ) 

9960*  900*  4)0 

OBNO" 

430 

CC(1  .5) 

« HY0MX3  * 0X(l4l,0) 

09N09 

GO  TO  900 

08DE7 

431 

CRO(l)  • 

AMlNl  ( 0X(t-l*0  ) * 0Y(I-1*J  ) 

) 

05N09 

aU(l *1) 

• CRO(l)  • POHXHY 

29JLe 

BB(t  *2) 

• 0.0 

08DE7 

CC(1*1) 

• 0X(1-1*J)  • HY0HX3 

08DE7 

CC(1*2) 

• -2.0  • ( HY0HX3  * OX(I-l.J) 

08DE7 

1 

4 POHXHY  * CRO(l)  ) 

29JL8 

CC(1 *3) 

■ HY0HX3  * OXd-l.J) 

0eDE7 

IF  ( CC(I.3)  ) 

436*  439*  436 

05Nn9 

439 

CC( 1 *3) 

■ 1.0 

09N09 

436 

00(1*1) 

• POHXHY  * CROd) 

0SNO9 

00(1*2) 

> 0.0 

080E7 

CE(1*1) 

• 0.0 

03JA8 

CO  TO  900 

060E7 

COEFFICIENTS  COr.PUTEO  AT  J ■ KvP3 

438 

IF  (l-l)  9980* 

439*  440 

05N09 

439 

60(1*2) 

• 0*0 

09N09 

e0( 1 *3) 

* 0.0 

080E7 

CC(1*3) 

■ 1,0 

170C9 

CC(1.4) 

• 0.0 

080E7 

CC(1*9). 

• 0*0 

080E7 

00(1*2) 

■ 0.0 

080F.7 

0D( 1 *3) 

• 0.0 

0R0F7 

EE(1*1) 

■ 0.0 

05JA8 

CO  TO  900 

080E7 

440 

IF  (MXP3-1)  9980*  447,  441 

09N09 

441 

CR0(4)  ■ 

AMINI  ( OX(I  ,J-1)  , 0Y(1  ,J-1) 

) 

09N09 

60(1*1) 

• CR0(4)  * POHXHY 

29JL8 

66(1*2) 

• -2.0  • ( POHXHY  • CR0(4) 

29JL8 

1 

4 HX0HY3  * 0Y(1.J-1>  ! 

080E7 

68(1*3) 

• POHXHY  » CR0(4) 

29JU8 

CC(i  ,2) 

••  O.C 

080E7 

CC(i *3) 

4 HX0HY3  * 0Y(1*J-1) 

080E7 

IF  ( CC(1,3)  ) 

9980*  442*  443 

09N09 

442 

CC(It3) 

• 1,0 

09N09 

58 


%J%J  uuuuuu  uwu 


44) 


11-2)  9980i  445.  444 


-RETAIN  RECURSIOH  COEFFICIENTS  TO  USE  AT  NEXT  J STEP 
AM2(|«U  « AHKItl) 

AHKItX)  « A(t*l) 

» IF  ( ML  ) SIS*  SOSt  SOS 

» 00  SlU  X c 1»  MXP3 

BM2(1»X)  • BMUItX) 

BHl(l.X)  • BUtX) 

CN2(ltKI  ■ CM1(I«K> 

CMlIl.Xl  • CU*K) 

> CONTINUE 

) CONTINUE 

-SOLVE  FUR  ALL  RECURSION  COEFFICIENTS  AND  RETAIN  THE  Adtll 
COEFFICIENT  AT  THIS  J STEP  IN  THE  ATII.J)  ARRAY 
JJ  • J 

CALL  MATRIX  ( LliJJtMXP3fKY.AA.BB«CCtDDtEE«FF*A.AMl> 

1 AM2»BtBMl*BH2»C>CMltCM2»D»E.«ML  > 

00  S20  1 • If  MXP3 

ATlIfJI  ■ Allfl) 

) CONTINUE 


esNOf 
OfOET 
OBOE  7 
080E7 
03JA8 
OSNOB 


444 

CCIlfl)-'  • 

0.0 

0SN09 

445 

• IF 

(MXP2-1I  9980, 

SOC,  446 

09N09 

446 

CCIlfSi  ■ 

0.0 

09N09 

CO 

TO  500 

0e0F7 

447 

BBIlfl) 

0.0 

0SN09 

B0(l,2)  > 

O.C 

0e0E7 

CCdfli  « 

0.0 

OBOE  7 

CC(i,2) 

0.0 

0B0E7 

CC(I,3I  • 

1.0 

170C9 

00(1,11  > 

0.0 

0BDE7 

00(1,2) 

0.0 

0BDE7 

EEllfl)  • 

0.0 

03JA8 

900 

CONTINUE 

04MV7 

-BEGIN  MAIN  SOLUTION 

SOI 

00 

SIS  1 « 1,  MXP3 

29N07 

710E7 

2I0C7 

29N07 

200C7 

210E7 

2I0E7 

210E7 

2I0F7 

04MY7 

30JE7 


0BAG7 

29N07 

19AP8 

200C7 

29N07 

200C7 


■TEST  FOR  MULTIPLE  LOADING  " 

IF  ZERO,  RETAIN  B AND  C COEFFICIENTS  ON  TAPE  1. 

IF  PARENT,  ALSO  RETAIN  0 AND  E MULTIPLIERS  ON  TAPE  2. 
IF  offspring,  MOVE  TAPE  I COMPLETELY  FORWARD  IN  STEPS. 


IF(  ml  ) S22,  S30,  S2S 
S22  READ  (1) 

GO  TO  600 

S2S  WRITE  (2)  11  On,AI,ClI,XI,  I«1,MXP3I,  X«1,MXP3  ) 

SSO  WRITE  (1)  i(  8 (l,X),  C (1,K),  I*l,KxP3l,  K>l,MxF^  I 
600  CONTINUE 


290C7 

2S0C7 

31AGC 

.200CT 

29N07 

4MY7 


COMPUTE  AND  PRINT  RESULTS 


C 


?9N07 

29N07 


DO  fciu  LL  * 1.  KYP3 
J « MYP4  - LL 

POSIT IO\  TAPt  1 FOR  READING  AMO  RCTRIlvE  THE  A RECURSION  COEFF 
BACKSPACE  1 16A07 

00  62S  I « 1.  MXP3  21JL7 

All. II  m ATII.JI  2eOE7 

S25  CONTINUE  22JE7 

C RETHlEVl  d AND  C RECURSION  COEFFICIENTS  AT  THIS  J STEP 

READ  in  (IB  (I.KI.  C II.K)*t>liMXP3l  . K>ltMXP3  I 29N07 

C REPOSITION  TAPE 

BACKSPACE  I I6AG7 

COMPUTE  deflections 

AMI  AND  AM2  ARC  NOW  TEMPS  USED  TO  REPRESENT  B»WP1  AND  C*WP2 
CALL  MATMPY  ILl,  MXP3.  1.  B*  WPl»  AMI  » ^ 26JAR 

CALL  MATMPY  ILl*  MXP3t  1.  C.  WP2.  AM2  ) ' 26JAR 

00  *30  I • 1.  MXP3  21JL7 

WIt.JI  ••  All. It  ♦ AMIII.II  ♦ AM2II.1)  ?9N07 

WP2II«lt  ■ WPlIltll  2SDE7 

UPlI  1 .1 1 ■ Wl I .J|  28DE7 

BSO  CONTINUE  AMY7 

«»0  CONTINUE  0AMY7 

SET  DEFLECTIONS  AT  CORNER  STATIONS  OUTSIDE  THE  BOUNDARIES 

«U*1I«  2.0  • Wll.21  - wa.3)  29N07 

WIHXP3.il  ■ 2.0*  Wl  MXP3.2I  -WIKXP3.3I  29N07 

WII.MYP3I  •2.0*  WII.MYP2I  -W  H.MYPD  29N07 

WIMXP3.MYP3I  ■ 2.0  • WIMXP3.MYP2I • -WCMXP3.MYP1 > 29N07 

COMPUTE  BENDING  MOMENTS.  REACTIONS  AND  TWISTING  MOMENTS 

00  73u  J ■2.  MYP2  2ASF9 

00  72U  I ■ 2.  HXP2  29N07 

CR0I2I  « AMINI  I OXII  .J  I * OYtl  »J  ) ) 29JLR 

WSUMl  « OOMX2  • I Wll-1.41  - 2.0  ♦ WII.Jl  ♦ WlUl.JI  I 06N07 

WSUM2  • OOHY2  • I WII.J-ll  - 2.0  * Wll.JI  ♦ Wll.J*!!  I P6N07 

BMXII.JI  « OXll.JI  • WSUMl  ♦ CR0I2I  • PR  • WSUM2  29JLA 

BMYII.JI  • OYll.JI  • WSUM2  ♦ CR0I2I  • PR  * WSUMl  29JLA 

720  . CONTINUE  OiSE7 

790  CONTINUE  08SE7 

—— OUTPUT  TABLE  8 

PRINT  11  1SSE6 

PRINT  1 2UL7 

PRINT  13.  I ANIINI.  N • 1.  AO  I 170C9 

PRINT  18.  NPROO.  I AN2INI.  N • 1.  18  I 170C9 

IF  I ML  I 888*  887*  887  08N09 

887  PRINT  39  8N09 

PRINT  Att*  SOM  IBACO 

GO  TO  889  08N09 

888  PRINT  A2.  KPROB  08N09 

PRINT  A8.  SOM  18AGC 

889  IF  IIOPPSI  9980.  892.  891  17AG0 

891  SDT2  • 8.0  / I THK  • THK  » I7AG0 

CO  TO  89S  08N09 
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oo<^r\r>  rt  r»  r»  n f\ 


•92  $0T2  ■ 1.0  17AC0 

•95  PRINT  AC.  SON  17ACC 

IF  ( iOPd.EO.l  > 17A00 

•PRINT  A9  17AG0 

CR  > l.OC-12  2IN07 

SUMR  •0.0  120C7 

00  960  J ■ 2.  MYP2  29N07 

IF  I lOPV.NE.T  > PRINT  6 7ACC 

C COMPUTE  ABiORDEO  LOADS 

JSTA  • J - 2 2fN07 

00  955  I « 2.  MXP2  29N07 

ISTA  • 1 - 2 29N07 

OBMX  ■ ( BMXU-l.J)  - 2.0  • BPXtl.JI  ♦ BMX(I+1,J»  > 02N06CDC 

C OBMX  « { OKX( !-l«JJ*1.000  - 2.CD0  * BMX(I.J>  ♦ BMX I I+l . J) JOOC9 J P- 

I • MYOHX  06N07Cf>r 

Cl  • 1.000  I • HYOKX  S0OC9IPM 

OBMY  - I BMY(I.J-1>  - 2.0  • BMY<1»J>  ♦ PKYd.J^l)  ) 02N06Cf»C 

C OBMY  • I 0MY(  1»J-U»1.0D0  - 2. COO  • BHYII.J)  ♦ RMyI  I . JaI  )3COC9lRr‘ 

1 • HXOKY  06NO7COC 

C l • 1.000  ) • HXOHY  ))D0C9lpy 

OTMX  • { Wd-l.J-n  * CMd.JJ  - W(I-l#J)  • ( CMd.Jj  07MY8 

1 ♦ CH(1.J+U  » ♦ wn-l.J+l»  * CHd.JtL)  07MY8 

2 - Wd«J-l>  * t CHd.J)  ♦ CHd  + l.JdA  WdtJl  07MY8 

S • ( CMd.JI  4 CHd.J+1)  ♦ CHn4l»J»  ♦ CHd  + l.J  07MYB 

A All  » - Wd.JAl)  * « CHdtJAl)  * CHClAUJAd  » 07MY8 

5 A NdAl»J-l)  * CHdAl.J)  - WdAl.J)  * I CHd  07MY8 


6 Al.jj  ♦ CHdAl.jAli  WdAl.jAlJ  • CHdAl.J  07MY8 


7 All  I • OOHXHY  . 06N07 

OTMY  • OTMX  07MY8 

OPX  • OOHX  • C PXd.JI  • Wd-l.J)  “ « PXCl.Jl  06NO7 
1 A PXdAl.JI  I • Wd,J)  A PXCIAI.J)  • WdAl.J)  I02N06COC 

c.  ^ A 0.000  A PXdAl.JI  I * Hd.Ji  A PXdAl.JI  * wd4:.jii3C0C9ipf 

OPY  • OOMY  .•  I PVd.J)  • «d.j-ll  - (PYdtJl  06N07 
1 A PYd.JAll  I • Wd.J)  A PYd.JAll  • WdiJAll  I02NO6COC 

1 A 0.000  A PYd.JAll  I * Wd.JI  A PYd.JAll  * Wd.JAll|30OC9I5t* 

COMPUTE  TWISTING  MOMENTS 

WSUM3  • ( Wd-l.J-ll  - Wd-l.JAl)  - WdAl.J-n  A 1A0E7 

1 MdAl.jAll  I * 0.0625  • OOHXHY  06N07 

TMX  ■ I CMd.JI  A CHd.JAl)  a'chIIaI.JI  ♦ 07MY8CftC 

TMX  • { CHd.JI  A CHd.JAl)  A CHdAl.JI  ♦ 0.0O0A30OC9lpy 

1 CHdAl.jAl)  ) * WSUM3  07MY8 

subtract  applied  load  FROM  SUM  OF  ABSORBEO  LOADS  TO  GET  REACT 

REACT  ■ OBMX  a QBHY  a OTMX  A OTMY  ~ OPX  - OPY  -Od.J)  06N07 

SUMwaioN  OF  REACTIONS  FOR  STATICS  CHECR 

SUMR  • SUMR  A react  17JAI 

IF!  REACT  • REACT  - ER  ) 905.  905.  906  10N07 

905  REACT  -0.0  10N07 


COMPUTE  PRINCIPAL  MOMENTS  OR  STRESSES 


1 


V 


906  BMA  • ( BHXd.J)  A BMYd.J)  I • 0.50  17JA8 

BMX  IS  BENDING  MOMENT  IN  X DIRECTION  (COMPRESSION.  IN  .TOP  I»>A  K 
TMX  is  twisting  MOMENT  IN  X DIRECTION  (ABOUT  Y AXIS)' 

SIGO  IS  THE  maximum  NUMERIC  VALUE  OF  PRINCIPAL  MOMENT  I a OR  - ) . 
iF  THICKNESS  SWITCH  IS  INPUT  IT  IS  THE  PLATE  STRESS  ANO  I S A 
FUR  TENSION  IN  BOTTOM  OF  PLATE 
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C CU0.4TCft  CLOCKWISE  BETA  ANCLES  ARE  POSITIVE  AND  ARE  MEASURED  FROM 

C the  X AXIS  TO  THE  DIRSCTIOK  OT  THE  LARGEST  PRINCIPAL  STRESS 

C oR  MOMENT  (POSITIVE  OR  NEGATIVE » 


TMY  • -TMX 

BMP  • BHXII tJl  > BMA 

BMR  • SORT  ( BMP  * BMP  ♦ TKV  • TKV  ) 

C BMR  • OSORT  < BMP  • BMP  ♦ TMY  • TMY  I 

BMO  ■ BMA  * BMR 

BMT  « BMA  BMR 

C TEST  TO  PRINT  ONLY  THE  MAXIMUM  VALUE 

IF  I BMA  ) 916«  918*  918 
916  PMMAX  • BMT 

IF  I BMP  I 940«  930*  920 

918  PMMAX  • BMO 

IF  ( BMP  I 920*  930*  940 

920  ALF  • TMY  / BMP 

ALF  > ATAN  ( ALF  ) • 97.29578 
C ALF  • OATAN  ( ALF  > • 97.29978 

IF  ( ALF  t 922*  924*  924 

922  THETA  »-ALF  - 180*0 

CO  TO  945 

924  THETA  « * 180.0  > ALF 

GO  TO  949 

93U  IF  ( THY  I 932*  934*  936 
932  theta  « ♦ 90*0 

GO  TO  949 

994  THETA  • 0.0 

GO  TO  945 

996  . THETA  « - 90.0 

GO  TO  949 

940  _ ALF  • TMY  / BMP 

ALF  • ATAN  I ALF  I • 97.29978 
C ALF  • OATAN  I ALF  > * 97.29978 

THETA  --ALF 

C—— CLOCKWISE  ANCLES  ARE  NEGATIVE 
949  BETA  • 0.9  • THETA 

SiCO  • PMMAX  • SOT2 
PSIGOIlfJ)  - SICO 
IF  ( IOP8.NE.1  I CO  TO  948 

IF  I ISTA.EO.IMX/21  .A.  JSTA.EO. (MYF2 » > GO  TO  948 
GO  TO  999 

948  PRINT  47  i ISTA.  JSTA*  W( 1 .Jl .BMXII *J» *BMYI I *J) *TMX  * REACT* 
1 SICO*  BETA 

999  CONTINUE 

960  CONTINUE 

PRINT  6 
PRINT  40.  SOH 
PRINT  9v.  SUMR 

IF  I NCT3  I 998C,  989.  970 
970  CONTINUE 

C 

C OUTPUTTABLE9 

C . 

PRINT  11 
PRINT  I 

PRINT  13*  I ANKII*  I*!*  40  t 


10SE9 

190F7 

10SF9CDC 

310C9IPM 

190E7 

190E7 

150E7 

loseo 

150E7 

19DE7 

150E7 

10SE9 

10SE9CDC 

170C9lP'-‘ 

19St9 

10SE9 

19Df7 

10SE9 

190E7 

10SE9 

10SE9 

I9DE7 

190E7 

150F7 

iosr9 

I90E7 

10SE9 

10SE9Cr>C 

170C9 I Rv 

10SC9 

190E7 

190E7 

26JE0  ■ 

17AGP 

17AGC 

17AG0 

oeoE7 

140F7 

nSE7 

09A07 

17AC0 

17AG0 

120E7 

26JCC 

26JE0 


24SC9 

24SF9 

170C9 
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MlMT  U.  NPROttt  ( Af«2(ll* 

I • 1, 

19  > 

170C9 

ir  1 ML  J 972.  971,  971 

26JC0 

971 

PRUT  63 

26JE0 

PRUT  63 

06N09 

00  TO  973 

2«jro 

972 

PRUT  66,  KPROB 

26JE0 

PRINT  65 

6N09 

973 

CONTINUE 

26jro 

CALL  SLPlT  ( W,  eHX,  8HY, 

PSIGO, 

LI,  l2,  IPOP,  I ) 

24JE0 

9t3 

CONTINUE 

26jrp 

CALL  TIC  TOC  (6t 

25SE6 

00  TO  1010 

76AG3 

9990 

PRUT  99 

29N07 

9990 

CONTINUE 

19MR5 

9999 

CONTINUE 

06MY3 

PRINT  11 

•MY3 

PRINT  1 

21JL7 

PRINT  13,  ( ANKN),  N • 1, 

60  ) 

. 

170C9 

CALL  TIC  TOC  <2» 

26SF6 

PRINT  19 

26AG3 

CNO 
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c 

c 


c 

c 

c 

c 


c 

c 


c 

c 


c 

c 

c 

c 


c 


SUBR-^UTINE 

MATRIX  (LI* 

JJtMXP3*MV 

.AA.CP.CC.CO.rF.FF.A 

*ami. 

29N07 

1 

AM2 

.C*C.'ll*CM2.D.f  .ML  ) 

19APB 

UOUlii-E  PP.E 

CISION  AA.  B3 

* CCf  OD. 

EC.  FF*  A.  AMI.  AK2. 

P,  BMl.  0M2 

.160C9)3- 

1 

C*  CMl 

• CH2«  Dt 

E 

l60C9Ia- 

OIMCmSION 

AA  (Ll*li 

• 

BB  Il1*3>  * 

CC  (LI. 51. 

OTJAA 

2 

00  (L1.3) 

• 

EC  (Ll*l)  * 

FF  (Ll.l). 

19AP8 

3 

A U.l*l) 

AMKlI.I). 

AM2(L1 .1) . 

29N07 

A 

B (LI. LI) 

• 

BHKlI.LD  . 

eK2(Ll*Ll) 

.29N07 

5 

C (LI. LI) 

• 

CMKlI.LD  . 

CM2(L1.L1) 

.29N07 

6 

D(L1*L1)  . 

E(L1*L1) 

2&JA8 

TEST  FOR  .-luLTJPLE  LOAOlflC.t  IF  OFFSPRING*  RETRIEVE  D AND  E 
RECURSION  MULTIPUIEPS  FROM  TAPE  2 


1F(  ML  I SOO*  S20«  S20  200C7 

SUO  READ  (2)  ((  D(I*K)  « EIItKI  * t>  1*HXP3)  * K>  1,MXP3  ) 200C7 

• CO  TO  550  20OC7 

compute  RECURSION  hUETIPLlER  £ 

52U  CALL  MATMYl  (LI  . MXP3  * MXP3  • AA  » BK2  • E)  t»N07 

CALL  MATAl  (LI  « MXP3  « C » SB  * E>  29N07 


— --compute  recursion  multiplier  -1/0*  C USED  AS’ A TEMPORARY 

CALL  HATHPY  (LI  * HXP3  • MXP3  . E * BKl  * Dl  • 29M07 

CALL  MATMYl  (LI  * MXP3  * MXP3  » AA  • CM2  • C ) 29M07 

00  535  K • 1.  MXP3  21JL7 

00  530  I • i,  MXP3  21oL7 

OdtKI  • - OdtXI  - Cd»X)  01DE7 

530  • CONTINUE  AMY7 

535  CONTINUE  22JF7 

CALL  MATS2  (LI*  MXP3*  0.  CC.  0 I 29N07 

COMPUTE  RECURSION  MULTIPLIER  0 

CALL  INVR5  ( Oi  LI*  MXP3,  JJ  ) 19AP8 

compute  recursion  COEFFICIENT  C 

CALL  MATM2  (LI  . MXP3  . 0 * EE  * C » 29N07 

COMPUTE  RECURSION  COEFFICIENT  B*  NOW  USING  BM2  AS  A TEMPORARY 

CALL  MATMPY  (LI  , MXP3  * MXP3  * E * CMl  * B I 29N07 

CALL  MATAl  (LI  . MXP3  i B * DO  * BM2 » . 21Df7 

CALL  MATMPY  (LI  * MXP3  * MXP3  * D * BK2  • 6 > 210E7 

compute  recursion  COEFFICIENT  A*  EE  USEO  AS  A TEMPORARY 

550  CALL  MATMPY  (LI  , MXP3  * I * E * AMI,  A I 29N07 

CAkL  MATMYl  (LI  * MXP3  , I * AA  * AM2*  EC  > 02UA8 

OO  55v  I • 1*  MXP3  21JL7 

EE(I*1I  > Ad*n  * EC  (1*11  - FF(I*n  300C9 

560  CONTINUE  0AMY7 

CALL  MATMPY  (LI  , MXP3  * 1 * 0 * EE  *A  I 02JA8 

RETURN  BAGT 

ENO  6AC7 
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SUdMwUTIHE  HaTMVI  <H1  . L2  * L * X * Y * Z I 
C UOUWlt  PHECIStON  X*  V«  2 

01KE<«Slori  X(M1,1I  • 2(K1«C)  • YIKl.Li 
00  200  I * 1 * L2 
00  100  J ■ 1«L 

2(ItJ)  • Xdtlt  • 

JU  CONTINUE 

0 CONTINUE 

RETURN 
END 


2iJLT 

U0C9I";- 

03JAR 

19AP^ 

13JL7 

03JAA 

13JL7 

19AP8 

13JLT 

13JL7 


subroutine  MATAl  (HI  • t2  t 2 t X3  » Y)  29JL7 

C . DOUBLE  PRECISION  X3t  2*  Y ‘ )60C9la’^ 

01HE^S10N  X3(Ml,3>t  2(Ki.Mn.  YIHltMl)  13JL7 

MMi  « L2  - 1 25JL7 

00  60  1 » 1 * L2  19,\PR 

OO  SO  J « 1 t L2  25JL7 

• • YlltJI  « 2(I»J»  ISJI.7 

50  CONTINUE  19JL7 

60  CONTINUE  19AP6 

00  100  1 > 2i  MMI  . 15UL7 

Yll«l-1)  • Yll.l-ll  •»  Xin.n  13A7 

Ylt*!)  ••  YdtII  ♦ X3(U2)  13A7 

Yd.dlJ  • Vd.l*ll  ♦ X3d»3»  19JI.7 

100  CONTINUE  13X7 

Yd»l»  • YdilJ  ♦ X3d.2»  13JL7 

Yd.2>  • Ydi2I  ♦ X3d»3»  13X7 

Y(L2tL2-lJ  • Y(L2.L2-1»  ♦ X3(L2>d  25X7 

Y(L2«L2I  ■ Y(L2.L2J  ♦ X3(L2t2)  25X7 

RETURN  13X7 

ENO  13X7 


■ ■*  "1 


SUBROUTINE  MATMPy  (Ml  • L2  . L * X t Y » 2 ) 25X7 

C DOUBLE  PRECISION  Xt  Y.  2 160C9Ipr< 

dimension  XIHlfKld  Y(MlfL  I*  2IKlfL  ) 27Jr7 

DO  300  1 " 1 I L2  19AP8 

00  200  M • l.L  19APA 

2diMl  * 0.0  27JE7 

00  100  K • 1 • L2  25X7 

2d,M)  • XdtX)  • Y(X,M)  ♦ 2d,M)  27JE7 

luO  CONTINUE  27Jr7 

2wU  CONTINUE  19AP8 

300  CONTINUE  19AP8 

RETURN  27JE7 

END  27JE7 
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SUilHyjTlNE  KATS2  (LI  . • 2 » X3  . Yi  30N07 

C UOJoLt  ('RL.JSlUII  2*  X3t  Y 160C9IP” 

UIMt<»SlUN  X3«Llii>  . 2(LltLri  • YtLltLtl  25JL7 

KH2  • M»  - 2 13JL7 

00  iOO  I « 3t  W2  13JL7 

Y(l,t-2»  * 2(1. 1-2)  - X3(I.1I  30N07 

2(1. t-D  - X3(I.2I  30H07 

Vd.n  « 2(i.n  - X3(|.3t  30NO7 

Y(l.!*l)  • 2(l.Ul)  - X3ll.*t  30N07 

Ytl.I+2)  « 2(1. I+2)  - X31I.5J  30N07 

100  CONTINUE  13JL7 

Y(l.l)  • 2(1.1)  - X3(1.3)  30NO7 

Y(i.2)  * 2(1.21  - X3(1.4>  30N07 

Y(l.3)  . 2(1.3)  - X3(1.5)  30N07 

Y(2.1)  « 2(2.1)  - X3(2.2)  30N07 

Y(2.2)  > 2{2.2)  - X3(2.3)  30N07 

Y(2.3)  « 2(2.3)  - X3(2.4)  30N07 

Y(2.4)  ■ 2(2.4)  - X3(2*»)  . 30N07 

Y(M1-1,M1-3)  - ?{hl-l,Ml-3)  - X3(K1-1,1)  30N07 

Y(Ml-l,Ml-2)  « 2(hl-l.Ml-2)  - X3(M1-1.2)  30N07 

Y(M1-1,M1-1  ) * 2(Kl-l.Ml-l)  - X3(V1-1.3)  -^ONOT 

YlXl-l.Ml)  « 2(K1-1»K1)  - X3IH1-1.4)  30N07 

Y(Ml,Ml-2)  • 2(M1,M1-2»  - X3(M1.1)  30N07 

Y1M1,HI-1)  > 2(M1,M1-1)  - X3(H1,2)  30N07 

YlMl.Ml)  . 2(H1,M1)  - X3tKl,3)  30N07 

*|7U«N  • - rjJL7 

tNO  , i3JL7 


SUitKOUTlNE  1((VR3  ( A,  tl.  L2.  JJ  1 19APn 

C OOUUtE  PRECISION  A.  S.  S2  160C9IPM 

OIMEiYStUN  Adl.LlI  2SJA8 

20  PORMmT  <///2«H  no  INVEKSE  EXISTS  JJ«t  15  . lOH  A(I.J)  « ) 300C9 

30  FORHAT  ( 1X.1CC10.3  ) 100C7 

31  FOKmaT  (//  5CH  NOTE.  THIS  IS  USUALLY  DUE  TO  AN  ILL-DEFINED  OR  31AG0 


1 / 52H  1N5TASLE  STRUCTURE  - CHECK  yOuR  INPUT  FOR  TIAGO 

2 ■ / 52H  ERRORS.  (BEWARE  WHEN  MODELING  CRACKS  THAT  31400 

3 / 52H  THE  STIFFNESS  AT  THE  CRACK  IS  GREATER  THAN  31AG0 

4 ! 52M  2ER0)  - ALSO  THIS  WILL  HAPPEN  IF  THE  SLAB  31AG0 

3 / »1h  has  buckled  due  to  excessive  axial  thrust  )3]AG0 


EP  • X.OE-IO  210E7 

00  185  I - I . L2  20OC7 

•KK  • 1 ♦ 1 200C7 

IF  ( AuSIAll.l)  ) - EP  ) 990.  990.  150  19AP8C0C 

C IF  I OABSIAII.I)  ) - EP  I 990.  99C.  150  300C9IPM 

150  S • 1 / Ad.n  200C7 

DO  160  J ■ 1 . L2  25 JAl 

All.J)  • Ad.J)  • S 200C7 

160  CONTINUE  200C7 

Ad.ll  <•  S 25 JAR 

00  ISO  J • 1 . L2  200C7 

IF  I J-1  ) 170.  ISO.  170  200C7 

170  $2  . AlJ.d  200C7 
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AlJtll  ■ 0.0 
DO  175  tC  . 1,  L2 

A<J,IC»  • A<J.K)  - S2  • A(1»K) 
175  COWTIHUE 

TW COHTTnUE 

U5  CONTINUE 

RETURN 

990  print  20,  JJ 

print  3u,  (I  AII,J»,  J.1,L2I  , I.l,  L2  ) 

PRINT  31 

CNO 


29JA4 
25JAA 
200C7 
20OC7 
20OC7 
200C7 
20OC7 
10  JAR 
.25JAR 
31AGQ 
2COC7 


C 


100 

200 


SUOROUTINE  MATM2  <M1  , L2 
DOUBLE  PRECISION  X,  Y,  2 
dimension  XIMl.iJ  , Z(M1«MII 
DO  200  I • I , L2 
DO  100  J • I , L2 

2(  1 » J)  ■ X(  J.  1 ) • YU  , j) 
CONTINUE 
CONTINUE 
return 

END 


» 7 * X , 2) 

, YIMl.  Ml: 


25JL7 

1&0C9I'*^' 

03JA8 

19AP8 

25JL7 

03UA8 

1JJL7 

I9APR 

13JL7 

19JL7 


i 

f 

( 
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iiiMiia<»^ni«iTWTiirTiirT8aaaaMBaiaiaaa 


rkr\  n rt  r>  r»  r\  r\  r\  r\  n 


SUu«OUTtNi  TIC  TOC  (J)  240Cfr 

THIS  ROUTINE  IS  S^tCIF ICmLLV  FoR  THE  COC6600.  WHEN  USING  THE 
ieM3*0/S0  StSTEM  THE  INDICATED  lOM  CARDS  WILL  CALL  THE  SUBROUTINE 
PRTIHE  TO  printout  THE  REQUIRED  TIKE, 


TIC  TOC  <ll  » compile  tike  20DE7 

TIC  TOC  <21  • elapsed  TM  TIME  JOOC9 

TIC  TUC  <31  « TIME  FOR  THIS  PROBLEM'  200E7 

TIC  Toe  <41  » TIME  FOR  THIS  PROBLEM  AND  ELAPSED  TM  TIME  300C9 

lU  FORrtMT<///3uX19HELAPSE0  TIME  ■ I5.«H  HINUTESF9. 3 .BH  SECONDS  ) 31AGC 

11  fuArt«T<///3*Xi5HCO^-PlLE  TIME  » «I5.8H  KINUTES*F9.3.8h  SFCONDS  ) 2SSE6 

12  roR«A7<///3;X24HTl.vc  FOR  THIS  PROBLEM  » *15. «H  MINUTES. F9,3*  2SSE6 

1 8H  SECONDS  » 25SE6 

I • J - 2 21JY7 

IF  I I-l  > 40»  30.  30  2IJY7 

30  FU  - F 25SFN 

AO  CALL  second  <Fl  25SEACnC 

AO  CALL  PR TIME  30OC9IPY 

00  TO  990  30OC9IOY, 

III  « F 25SE6 

II  • III  / SO  255fS 

FI2  • F - ll»SO  ' 25SE6 

IF  ( I J 50.  70.  SO  24JL7 

50  PRINT  ll,  II,  FI2  21JY7 

00  TO  99d  255E6 

M FI3  • A-  JS.  FU  255E6 

12  • FI3  f SO  25SCS 

FIS  • FI3  - 12*60  25SF6 

PRINT  12,  12,  Ft3  25SES 

IF  < l-l  \ 99C,  990,  70  21JY7 

70  PRINT  lu,  ll,  F12  . 2UY7 

990  CONTINUE  OSSE? 

RETURN  25SFS 

END  25Srs 


SUbMUOTINE  PLOP  » J.  LI.  L2 . NCT6.  KCT7.  MX,  MY  I 

C - tHlS  SJBPOUTJf.E  ACCEPTS  SPECIAL  LOAD  H-AITERf.S  ANO  TmEIR 

f PLACE.'^LMTS.  Ai'<0  GEbCPAl^S  THf  APPHOPRIATE  LOAD  VALUES  FOR 

c Solution  oy  the  ■••ain  phoorak. 

UIME..SIU1N  lM6i  9.121.  JH16<  9,121  i 
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QP|  9,i2) . N0(9t . OtLl  .L2I  , 


JN17C 


DATA  OP 

60  ruRM«T  I 

61  format  ( 

64  format  I 
1 

2 /.  5X 

3 

4 / 26X 

65  format 

66  format 

67  format 
7U  format 

71  format 

72  format 
7a  Format 

1 
2 
3 

•0  format  (///50H 
50H 
50H 
50H 


1N17I  71 , 

/ 108  "MD.O  / 

12. 2X  12.  2X  2413  > 

8X  12F6.0  I 

5X50H  pattern  NuH  OF 
30H  CONCENTRATED  LOADS 
50H  NUM  LOADS 

X 


PATTERN  COORDINATES  .‘ND 


30H  Y XV 

iohreference  • 

/5X1S.  SX  U.  4X  61  15.  14  > ) 

23X  6(  3X  F6.0  ) > 

26X  6(  13. IX. 13  1 I 

//50H  TABLE  7,  PLACEMENTS  OF  SPECIAL  LOAD  PATTERNS  / 
1615  1 

2IIC.  IX  7(  15.14  I I 

50H  pattern  NUM  of  LOCATION  OF  REFE 

30HRENCE  LOAD  ISLAB  COORDINATES)  / 


50h 

35H 


NUM  PLACEMENTS  X Y X Y > 
Y XY  XY  XY/l 
NUM  OF  wheel  loads  APPLIED  TO  THE  SLAB 
SUM  OF  WHEEL  LOADS  APPLIED  TO  THE  SLAR 
NUM  OF  WHEEL  LOADS  PLACED  OUTSIDE  SLAB 
SUM  OF  WHEEL  LOADS  PLACED  OUTSIDE  SLAB 


lOU 

110 

115 


1 Ilu,/ 

2 CIO. 3, 

3 Ilu./ 

4 ElU.3  I 
PRINT  64 

IF  ( NCT6.E0.0  I 00  TO 
NCT602  « NCT6  / 2 
00  110  N * 1,  NCT602 

READ  60,  NP,  NL,  ( INloiNP.J),  JN16INP,J),  J 
NOINP)  ■ NL 

READ  61.  ( OPINP.II.  I < 1.  NL  I 

IF  ( NL  .GE.  12  I GO  TO  110 

NLPl  • NL  ♦ 1 
. 00  100  I - NLPl , 12 
OPINP.I ) < 0.0 
CONTINUE 


115 


NL  ) 


120 

122 


DO  150  NP  ■ 1.  9 
00  120  I » 1,  12 
IF  ( OPINP.I I .EO. 

NL  ••  NOINP  I 
GO  TO  122 
CONTINUE 
GO  TO  15C 
IF  I NL.GT.6  ) 


0.  ) GO  TO  120 


GO  TO  140 


76J£n 
18AG0 
I SAGO 
18AG? 
76JE0 
26JE0 
07SEP 
26JE0 
26JC0 
17AG0 
26JEr 
17AG0 
18AG0 
17AG0 
17AGP 
03JL0 
26JE0 
) 17AGC 
26JL0 
26JE0 
17.  GO 
17AG0 
17AG0 
26JE0 
17AG0 
26JEP 
17AG0 
26JE0 
26JE0 
26JE0 
07SE0 
26JE0 
07SE0 
26JE0 
25JE0 
26JE0 
07SE0 
07SFC 
07SE0 
07SE0 
07SE0 
07SE0 
07SEf’ 
07SE0 
07SEC 
07SEC 
07SE0 
07SE0 
07SE0 


PRI«T 

65, 

NP,  NL,  1 

IN16INP»J). 

JM6INP,J1  . J • 

1,  NL  1 

75Jro 

PRINT 

66. 

1 UPINP.I) 

, I ■ 1,  nl 

1 

26JE0 

GO  TO 

150 

26JfO 

140  PRIhT 

65, 

Up , NL , 1 

INI6INP,J1, 

JN16INP,J),  J - 

1,6  ) 

26JE0 

PRINT 

66  . 

1 OPINP.I I 

, I *•  1 , 6 

) 

26JE0 

PRINT 

67, 

1 {M6INP, 

Jl  , JM6CNP 

,J),  J . 7,  NL  1 

26JC0 

print 

66 . 

1 OPiNP.n 

. 1 « 7,  NL 

I 

26JE0 

69 


I 


I 


150 

CONTIKU'E 

26jrO 

155 

QPSUH  s 0.0 

31  AGO 

OPOSUM  • O.C 

PSJEO 

LOP  ■ 0 

26JF0 

LOPO  • 0 

26JF.0 

PRINT  70 

26JE0 

PRINT  74 

26JE0 

UO  300  N « If  NCT7 

• 

2SJE0 

READ  71f  NPf  NSPf  ( INHI  U *JN17( !)  f I 

» If  NSP  » 

26JE0 

PRINT  72tNPf  NSPf  < lNl7«HfJ.n7n)f  1 

« If  NSP  » 

26JE0 

NL  « NQ(NP> 

26JE0 

00  290  IP  « If  NSP 

29jrp 

DO  280  IL  < If  NL 

29JE0 

1 • INUtNPflL)  * IN17(IP) 

18AG0 

J • JNl6(NPfILI  ♦ JNl7fIP> 

18AG0 

IF  ( l.LT.O  *0.  I.GT.MX  *0.  J.LT.O 

.0.  J.CT.HY  I 

GO  TO  2fO  29JEC 

1 > I 4 2 

18A6P 

J • J ♦ 2 

. 18AG0 

OPSUH  • OPSUH  * OPtNPfILi 

29JEe 

LOP  • LQP  ♦ 1 

?sjro 

QtlfJ)  > O(lfJ)  4 OPINPfIL) 

29JE0 

GO  TO  260 

26JE0 

2«U 

OPOSUH  « QPOSUM  4 GPINPflLI 

29JE0 

LOPO  • LOPO  4 1 

. 

26JE0 

2tO 

CONTINUE 

26JE0 

290 

CONTINUE 

, ' 

2SJE0 

500 

CONTINUE 

2SJE0 

PRINT  UOf  LOPf  OPSU^if  LOPOf  CPOSUK 

26JE0 

RETURN 

26JE0 

END 

26JE0 
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SUtMuUTINC  COPS  I CH«  OX*  OY*  S*  Lit  L?  ) HJ'O 

COHMwft  ! CUPS  / EM*  SXK*  fOPO*  lOPC*  lOPS*  MX*  HY  ITAG" 

C - - - - IMIS  SU8P0UTIKE  GCNrPATCS  THE  OPTIONAL  PAVEMENT  STIEPNESS  18AOP 
C CONSTANTS.  APPROPRIATE  OUARTER  AND  MALP  VALUES  ARP  18AGC 

C GENERATED  AT  THE  EDGES*  1IAG9 


COMPM* 

/ 

SPLT  / MIS*  MY*  THK*  NCT9* 

NOE*  PR 

2AJEP 

DIMENSION  0X(Ll*L2t*  0YILl*L2t*  CHIL1.L2I*  SILli 

L2I 

26JE0 

FORMAT 

I 

/ 90H  table  a*  ( CONTO  > COMPUTED  OPTIONAL  PAVEMENT 

i8AG0 

19H  STIFFNESSES 

ISAGO 

// 

97m  prom  thru 

OX 

OY 

18AG0 

lOX  29H  $ 

C 

07SEC 

// 

TX 

9M1  1 I9*IA*2X2H9*PX2M0*2CX2ri0.8XPir*9* 

18AG0 

/ 

7X 

9H0  0 I9.IA*  2E11.9* 

nx  Ell. 9*  ah 

0. 

/ 

7X 

9Hl  0 I9.IA*  2Cn.9* 

IIX  Ell. 9*  AM 

0. 

18ACP 

/ 

7X 

9H0  1 I9*IA*  2E11.9* 

nx  Ell. 9*  ah 

c. 

18AC0 

/ 

7X 

9H1  1 I9.IA*  2C11.9* 

nx  Ell. 9*  AH 

0*  I 

18AG0 

OXGOA  « 0*0 

26JE0 

SGOA  • 0*0 

< 

2SJE0 

CMGEN  *0.0 

2,«JE0 

MYPI  « MY  ♦ 1 

, 

2SJC0 

MYP2  ■ MY  ♦ 2 

2SJE0 

MXP2  • MX  ♦ 2 

2SJro 

MXPl  • MX  4 1 

. 

2SJE0 

MXHl  ■ MX  •*  1 

2SJ»0 

MYMI  • MY  - 1 

2SJfO 

IP 

I 

lOPO.NE.l  ) GO  TO  no 

" - ■*** 

09Ae 

OXGEN  • ( EM  • THR4P9  1 / 

I 12.  ••  < 1*  - RR«PR  I » 

26JE0 

OXGOA  • OXGEN  • *29 

2SJE0 

loo 


IP  ( 


I 


tOPC.NE.l  > GO  TO  no 
CHGEN  • OXGEN  • ( 1*  • PR 
lOPS.NE.l  » GO  TO  200 
SGOA  ■ SXX  • MX  • MV  • *25 

200  PRINT  AO*  MX*  MY*  CMGEN*  MX*MY*  DX60A*  DXCOA*  SCOa 


no 


IP  I 


HY*  OXGOA*  OXGOA* 
MYHl*  OXGOA*  OXGOA* 


SGOA 

SGOA 


290 

900 


910 


920 


9A0 


1 • MXHl, 

2 ■ * MX* 

9 » MXMI,  MYMI*  OXGOA*  OXGOA.  SGOA 

IP  ( lOPC.NE.l  > CO  TO  900 
00  290  I • 9*  HXP2 
. DO  290  J « 9*  MYP2 

CMIl.JI  • CMII.JI  ■»  CMGEN 
CONTINUE 

OO  910  I • 2*  MXP2 

00  910  J • 2*  MYP2 

OXd.J)  • OXII.J) 

OYII.J)  ■ OY{|*U) 

$(l*Jt  • SII.JI  ♦ 

00  920  I • 9*  MXPl 
00  920  J • 2*  MYP2 

DXII.J)  • OXII.JI 
OY(l«J)  • OYd.JI 
Sd.JI  > Sd.JI  ♦ 

00  9A0  I ■ 2*  MXP2 
00  9A0  J • 9*  MYPI 

OXd*Ji  ■ OXd.JI 
OYdtJ)  > OYd.J) 

Sd«Ji  ■ Sd.JI  ♦ 


OXGOA 

OXGOA 


♦ 

♦ 

SGOA 


♦ OXGOA 

♦ OXGOA 
SGOA 


♦ OXGOA 

* OXGOA 
SGOA 


2SJE0 

2SJE0 

2SJC0 

•.2SJE0 

2SJE0 

2SJE0 

2SJE0 

2AJC0 

2SJE0 

2SJF0 

2SJE0 

2AJE0 

2SJC0 

2SJE0 

2SJE0 

2SJE0 

2iJE0 

2SJC0 

2SJE0 

2SJfO 

2SJE9 

PAJPO 

2AJP0 

2SJP0 

2SJF0 

2AJIC0 

UJi9 

2SJC0 
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w 


f 


00  330  I > 3t  HXPl  20JC0 

00  330  J • 3*  MVPl 

0X(l*J»  « OXIUJt  ♦ OXCO*  26JfO  • 

OYIttJt  • OY(i*Jt  * 0XG04  2«JeC 

330  $(t«JI  • * 5004  2OJE0 

RETURN  20JF0 

CNO  20JE0 


1 

f 

\ 

i 


i 

i 

1 

i 

( 


SUORUUTINE  ZOTl  ( XX*  YY*  NERO*  10  » 2iUE0 

C • • • THIS  ROUTINE  ( MHE^  ACTIVATED  I CRIVES  A CALCOMR  703 
C IF  YOU  HAVE  A 703  AND  DESIRE  THE  USE  OF  THIS  ROUT  INF 

C .CONTACT  FRANK  L CHORES  * AUSTIN  RESEARCH  CNCINCERS  INC* 

C 3120  MANOR  RO.  AUSTIN*  TEX  71723 

RETURN 

IRO  ^ . 


•UORUUTINC  SFLT02  ( KEY*  NCNO* 
OINCNSION  XINCNO)*  OUHd) 

COMMON  /PLOT/  II*  12*  Jit  J2 
IF  ( KCY.EO*!  I 
IF  ( XCY*C0.2  I 
return 
CNO 


X*  OUM  I 

CAU  SRL0T3  1 X*  NCNO*  OUM 
CALL  SRL0T4  ( X*  OUM*  NCNO 


02UC0 

02JC0 

23N09 

I0N09 

lfN09 

I30C9 

130C9 


I 
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SUI»AOUfiM£  SLPLT  I W*  BKX*  SMV*  SICO*  Ll«  L2«  IPOP*  IC  t 

CU>^«»  / ZOT  / LOPt  MC*  IPOLL*  MOP 

COMMON  / SPLT  / MX*  MV«  IHK.*  NCT3*  KOE«  PR*  SOM 

COMMUN  /PLOT/  !!•  12 f Jl*  J2 

OlHei^StON  iril3(  10  l<  tN2S(  10  !•  JNl3(  10  I*  JN23I  10  I* 

1 KASewdOl  f KASEXdOl  • KASPYtlOt  • KASCPdOl  * 

2 TEMX  Moot*  TEMY  I300t*  XXI 300 1 » 

S WIL1*L2I  « SMX(L1»L2I*  BMYILl*L2)*  StCOILULZt 

DATA  IDl*  102*  103*  1U4  /lOHOEELECTlON.  lOHEENO  MOM  X* 

1 lOHBENO  MOM  Y*  lOH  SIGO  / 

46  format  ISX*  2dX*i2*lX*I3l*22X*E11.3> 

34  format  I 41  2X>  13  i«  41  4X*  d t I 

37  FoRnAT  I 3.1*  21  IX*  12*  IX*  13  I*  4X*  12*  IIX*  12*  2dlX*l2}  I 


63  FuRMaT  (///«13X«26HX  MOMENTS  ONLY*  BETkEEN  I 

1 «H  ) AND  I *13*  1H*«I3*2H  ) * ///• 

2 30H  X * Y X MOMENT  « / 

44  format  I 13X*  12*  IX*  13*  ICX*  C1C*3  I 

44  FORMAT  I 15X*  12*  IX*  13*  tOX*  IPEIO.3  > 

47  format  <///*1SX*24HY  MOMENTS  ONLY*  BETWEEN  I 


1 IH  t AND  I *13*  lH*.t3*2H  t » ///» 

2 30H  ‘ X . Y V MOMENT  ./ 

41  Format  I 13X*  12*  IX*  13*  lOX*  C)0.3  ) 

41  FORMAT  I 13X*  12.  IX*  13*  lOX*  lPElO.3  ) 

43  FOR.MAT  l///•13X*32HB0TH  X ANO  Y MOMENTS*  BETWEEN  I 

1 1N**13*  IH  I AND  I *13*  lH.*13t  2H  1 *///* 

2 22H  X * Y X MOMENT  * 2RX*  lOH  Y MOMENT 

74  FORMAT  I 13X*  12*  IX*  13*  ICX*  C1C.3*  3X*  £10*3  I 

74  format  I 13X.  12*  |X*  I3*  ICX*  IPCIO.3*  3X*  lPeiO*3  I 

•4  format  I///.13X.24HOEFLECT10NS*  BETWEEN  I *t3 

1 BH  ) ANO  I *13*  1H**I3*2H  t » ///* 

2 3QH  X * Y DEFLECTION  */ 

ifFORMAT  I///*13X»1CHPR1NC1PAL  AB«14HES»  BETwrCN  . I 


*t3*lH«»t1* 


•t3*lH**t3* 


*t3»lH*»t3* 


90  format  I//  4SH 


1H**|3*  SH  I AND  I *13*  1H**13»  2M  » 
13H  X • Y A4  */ 


•///. 


••••  CAUTION,  total  number  OF  SPECIFIED 


1 lOH  POINTS  IS  • 13*  /• 

2 40H  SHOULD  NOT  BE  GREATER  THAN  300  *//) 

IF  I iC.NE.O  I GO  TO  BOO 

00  100  I ■ 1*  300 
XXdt  ■ 1 >1 
10  CONTINUE 

LOP  » -1 
NC  - 10 
IROLL  • 1 
MOP  • 0 

— IMB0TTABLI3  . . 

00  173  N • 1*  NCT3 

READ  34*  INlSiNt*  JN13INI*  IN23IN)*  JN23INI*  XASCWIN)*  KASEXiNl* 
1 KASCYINI*  KASEPINt 

PRINT  37*  IN13(N|*  JN13INI*  IN23(N)»  JN23INI*  KMEWtNt*  KABEXTMi* 
1 RASEYINI*  KASEPINt 

NCNO  • dN23INt-IN13INt4t}»IJN23INI’>JN13IN)4lt 
IF  I |N13INt  > iN23INI  » 134*  134*  133 
13  NOE  > NOE  4 1 


02JE0 

tfJCBCOC 

3iAG0 

2WY0 

19MYB 

12MYP 

13JE0 

IBMYO 

ISJCOCOC 

13JC0C0C 

23(>C7C0C 

20APO 

20APe 

0SNO9 

24N09 

23NO* 

23BPBCOC 

I70CB1PM 

0SNO9 

24N09 

25N03 

2Sse«coc 

ITOCBtpM 

03NOB 

24N03 

24NOf 

30OC3COC 

SOOCBI-M 

20AP0 

lOAPO 

20AP0 

31  AGO 

27AP0 

3XA60 

03MR0 

03MR0 

03MRB 

12MY0 

13JE0COC 

IfJFPCOC 

l3JECCf»r 

13JE0C0C 

13JF0CMC 

I3JEOCOC 

13JE0C0C 

12MY0 

12MY0 

IBMYO 

20AP0 

20AP0 

20AP0 

lOAPO 

eSMRO 

e20C« 

24Be« 
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WWW 


1»A 

If  1 J<<t3(MI  > JN2)IM  t 

190t 

150*  155 

2osr9 

M 

NOE  • NOE  ♦ 1 

• 

20SE9 

150 

IF  I 1N23IN>  - KX  1 

150 1 

150* 

157 

20Sf9 

I5T 

NOE  • NOE  1 

20Sf9 

I5t 

If  I JN23INI  - MY  1 

100 1 

100* 

159 

20SC9 

159 

NOE  • NOf  ♦ 1 

* 

20SE9 

100 

If  ( KASEXINI  - 1 1 

102.* 

102* 

101 

090C9 

lOI 

NOE  • Noe  ♦ 1 

090C9 

102 

If  1 KASEYINI  - 1 1 

100* 

loo* 

103 

090C9 

105 

NOE  • NOE  ♦ 1 

090C9 

100 

IF  ( RASEWINt  * 1 1 

100 1 

100* 

105 

fOAPO 

105 

NOE  • NOE  ♦ 1 

20AP0 

100 

If  1 KASCPINI  > 1 ) 

lOtt 

100* 

lOT 

20APe 

lOT 

NOE  > NOE  ♦ 1 

20AP0 

lot 

If  I NENO  - 300  1 ITOt  ITOt  ITS 

XOAPO 

IT5 

NOE  • NOE  ♦ 1 

20AP0 

MINT  90t  NENO 

05MN0 

ITO 

CONTINUE 

lOAPO 

1T5 

CONTINUE 

20AP0 

GO  TO  9000 

lIMYO 

— — OUmffTAMJt 

•Od  COMTIMUC 

If  ( I^.IO.S  > MO^  • 1 
IPOl*  • .1 W - I 
00  ft!  M • It  KC1S 
ir  ( KASCWIMI  • 1 I Oltt  tott  ffto  * 
tut  filM  Mt  iMISINIt  JMISINIt  INfflNIt  JNISIN> 
II  • INIKNI  ♦ 2 
--  II  • INISIMI  ♦ I 

J1  • JNtSfNI  ♦ I 
Jl  • JNIKNI  ♦ I 
NtNO  • 

R • 0 

If  < U2-J1I*(II-UI  I 99**  mt  000 
•M  00  006  J • Jit  Jl 
00  to*  I • lit  II 

' TCMXIKI  • WdtJI 
too  CONTINOC 

•0  TO  til 

tot  00  tlO  I • lit  12 

00  tio  J • Jit  Jl 

R • R ♦ I 
TIMXIRI  •WdtJI 
tit  CONTINOC 

ttl  if  ( IfOf  I ilOt  lUt  tit 

tl4  CALL  ifLTOI  ( It  NCNOt  TfMXt  It.  > 

■ 60  TO  «|tf 

tU  . CONTIKOe 

10  • 101 

CMt  lOT  I I XXt  TCNXt  NCMOt  10.  I 
Ofl  • NINO  ♦ t 
00  tlT  I • It  Nf| 
tlT  XXdl  ■ I ' I 


ItJfO 

ItJCOCOC 

lOHve 

lOMVO 

tOHYO 

lOHYO 

lOAfO 

lOAfO 

lOAM 

lOAfe 

ItAfO 

20AW) 

lONVO 

lONYC 

20NVO 

lOAfO 

lOAfO 

lOMVO 

I0MV9 

lOWVO 

lOWYC* 

lOAfO 

lOAfO 

lOMYO 

IOMYO 

IWO 

lOWYO 

llJft 

IfJCOCnC 

ISJttCOC 

ItJCt 

ItJft 

Ifjft 


ir  ( t^O^.CO.l  i CO  TO  •!« 

•It  CONTIKUe 

820  CUMTINUC 

DO  tlw  N « 1«  NCTS 

IF  I KASCXtNI  • KASCVINI  - 1 I tSAt  822*  tttO 
•22  FRlKT  49*  INISCMI*  tN23(N|*  JN23IN» 

CU  TO  432 

•24  IF  I KASEXiNI  • I > ii$*  »26»  99M 
•24  FftlM  43t  INISINI*  JN13(Nt*  tN23<Mi*  JN23(N) 

GO  TO  432 

•24  IF  ( KASEVtNI  * 1 I •40*  430*  9940 

•30  PRINT  47*  IN13(Nt*  •<N13(K|»  tN23iN}*  JN23«N> 

432  CONTINUE 

11  ■ IN13(NI  4 2 

• JNl3fNI  4 2 

12  > IN23(Nt  4 2 

J2  • JN23(Ni  4 2 

NENO  f <l2-ll.41l«IJ2-Jl4tl 
R • 0 

00  434  J • Jl»  J2 
00  434  I > 11*  12 
K ■ K 4 I 

TCNXIKI  • tMXIUJI 
tCMYIKl  • OMYIltJI 

•34  CONTINUE  

•34  CONTINUE 

IF  I IFOF  I •34*  449*  449 

•34  IF  I KASEXINI  • KASEVINI  • 1 I 442i  440*  9940 
•40  CALC  SFI>T02  I 2*  NENUt  TENX*  TEMV  I 
GO  TO  440 

•42  IF  I KASEXINI  • | I SAA*  443*  9940 

.443  CALI  SFI-T02  I I*  NENO*  TENX*  20.  I 
GO  TO  440 

•44  IF  I RASEVINI  • I I 9940*  444*  9940 

•44  IF  I IJ2>Jll-ll2>ni  I 494*  494*  444 

444  K > 0 

00  492  I • It*  12 

00  490  J • Jt*  J2 

K > K ♦ 1 

TENTIKI  • 41191 1*JI 
•10  CONTINUE 

442  CONTIMUt 

•94  CAU  SFLT02  I I*  NENO*  TENV*  20*  I 
•94  CONTINUE 

•49  COllTiNUC 

IF  I KASEXINI.EG*!  I 10  ••  102 

IF  I KASEXINI*C0.1  I CAU  EOT  I I XX*  TENX*  NENO*  10  ) 

NF2  • NENO  ♦ 2 
00  470  I • 1*  NF2 
•70  XXII I • I > I 

. IF  I KASEVINI*E0*1  I 10  • 103  ' ^ . « * 

IF  I KASEVINl.EO.l  I CALL  207  | I XX*  TVMT*  NENO*  10  I 
OO  472  I • I*  NF2 
.472  XXIIt  • I - 1 

IF  I IFOO«EO*t  I GO  TO  434  s 

440  CONTINUE  ’ 


fonyo 
fonyc 
20Nve 
09N09 
20NYC 
20NV0 
2IW70 
FO»iye 
20Nr0 
20X90 
20NVe 
20M90 
20NY0 
24SF9 
24419 
24Sf9 
24SE9 
ltN09 
S9N09 
20N90 
20X90 
29N09 
29N09 
29N09. 
20X90 
20M90 
21N90 
30N90 
29090 
29N90 
20090 
29090 
29090 
24090 
20X90 
20XYO 
'20OV0 
20OYO 
12029 
120E9  . 

20090 
20090 
29090 
20090 
tSJEO 
I9UE0CI>C 
I9JC0COC 
IfJEO 
IfJEO 
19UE0 
I9JE4COC 
ISJEOCflC 
ItJEO 
IfUCO 
24094 
9004 
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M 914  N • 1*  NCTI 

20MY(* 

IP  1 RASCPtNI  • t 1 914*  Ilf*  9910 

204ve 

Mt  PRINT  •?«  SOM*  INlStNI*  JNlSINl*  IN29tNI*  JMtStNI*.  SOM 

07SC0 

1 

11  • INlitMt  ♦ 2 

12Mve 

12  • IN2S<N>  *2 

12MV0 

i 

Jl  • JNlliNt  * 2 

12MY0 

( 

J2  • ^2flNI  4 2 

12NY0 

1 

NCNO  > lt2*>ll<»ll«<J2*Jmt 

12NY0 

R • 0 

12MY0 

i 

02JC9 

1 

••I 

00  IIS  J • Jl*  J2 

02J»e 

j 

00  IIS  1 • 11*  12 

02JC0 

( 

R • R 4 1 

02jeo 

1 

TCNRIRI  • llOOIltJl 

02JC0 

m 

CONTINUC  • 

OIJCO 

M TO  190 

02jeo 

( 

MT 

M Ml  1 • 11*  12 

•2Jf5 

! 

M Ml  J • Jl«  J2 

02JC0 

i 

R • R 4 i 

eijco 

1 

TCNXIRt  a 1100(1 *JI 

eiJco 

I 

Ml 

CONTINUE 

02JCO 

1 

•M 

CONTINUE 

exjEo 

! 

IP  1 IPOP  1 9II«  970*  970 

CONYO 

! 

Ml  CAU  SPLT02  1 1*  NCNO*  TfMX*  tO.  1 

2INY0 

i 

M to  914 

COMYO 

1 

fTI 

CONTINUE 

ISJfO 

i 

, 

10  • 104  — 

19JCOCOC 

CALt  XOT  1 ( XX*  TCMXt  NCNO*  tO  1 

ifjcocnc 

j 

NP2  • NCNO  4 2 

19JC0 

1 

M 972  1 • 1*  NP2 

19JC0 

! 

ft! 

XXIll  • 1 > 1 

19JC0 

! 

« 

IP  1 IPOP.CO.l  1 00  TO  901 

2IMY0 

1 

i 

MA 

CONTINUC 

IIMYO 

1 

Ml 

CONTIWIE 

20NYO 

j 

fVM  RITUIIN 

lONYO 

9fM 

CONTINUC 

19MY0 

INO 

j 

tUlAuUtlMC  S^LOT  SIX*  rei.?*  hIOTH  I ItNOt 

C • • • • THC  tATEST  REVISION  OAi'-  *i'^  THIS  ROUTINE  IS.  - - 11  DEC  S»  REVISEf) 


DIMENSION  XINENOI*  SPACCIISI*  Ttk 

)0OC9 

COMMON  /PLOT/  11*  12*  Jl«  J2 

29N09 

data  S^ACE  / 15*AH  /*  SYK8  / 

4H*  «4M  • *4H  • *4H  • / 

100C9 

THE  NEAt«  'MtUFS  OF  X 

100C9 

DCOINNINO  WITH  THE  INITIAL  VALUE  TRANSEERRED* 

100C9 

THE  PAPER  SHOULD  BE  POSITIONED  PROPERLY  At«0  ALL 

100C9 

HCAOINCS  PRINTED  BEFORE  CALLlNC.  F*  L*  E* 

lOOCt 

••••  INPUT  “ X*  THC  FUNCTION  TO 

BE  PLOTTED 

100C9 

NENO*  THE  NUMBCK  OF  X 

TO  BE  PLOTTED 

.140C9 

WIDTH*  WIDTH  OF  PLOTI  LESS  THAN  41  » 

t40C9 

••••  OUTPUT-  NO  ACTUAL  VALUES  ARE  RETURNED  TO  THE 

100C9 

CALLING  ROUTINE  - THE  VALUES  ARE  PRINTED 

I00C9 

AND  PLOTTED  VERTICALLY. 

tOOC9 

10  FORMAT  1 5X*  12*  IX*  IS*  IX*  £10. 

S*14A4  1 

ISflOf 

IS  format  1 5X*l2*lX*tS*lX*ElO.S  1 

ISM09 

IF  1 NENO  .LE.  0 1 

GO  TO  9*0 

IINOB 

IF  1 W10TH.GT.40.  .OR.  WIDTH. 

LT.l.  ) WIDTH  • 40. 

1S0C9 

ISXP  • WIDTH  / 4 ♦ 1 

* 

140C9 

SVHO  • SYHBIll 

100C9 

OMEGA  « Xtll 

looet 

THETA  • Xlll 

100C9 

IF  1 NENO  .EO*  1 i 

GO  TO  120 

29IM>9 

00  SO  1*2*  NENO 

• 

leoct 

IF  1 OMECA.LT.XIII  1 

OMEGA  • XfS) 

lOOCt 

IP  1 THETA.GT.XIII  1 

THETA  •XII) 

100C9 

SO  CONTINUE 

100C9 

IF  1 OMEGA.EO. THETA  1 

GO  TO  40 

29N09 

sigma  • 1 WIDTH  - 1.  } / 

1 OMEGA  • THETA  ) 

tOOCf 

SGT  CONTINUE 

tSMOf 

' IF  1 <U2-J1-I2-»I1I  .GT.  0 1 

GO  TO  ISO 

ISOfB 

1 • 0 

tSMOt 

00  110  US  • Jl*  J2 

• * 

I9N09 

JSTA  • JS  - 1 

2SN09 

00  100  IS  • 11*  11 

2SN09 

ISTA  • IS  - 2 

ISN09 

1 • 1 ♦ 1 ♦ 

ISN09 

IF  I OMEGA.EO. THETA  I 

GO  TO  BO 

100C9 

1 BETA  > SIGMA  • 1 XII)  > 

THETA  ) ♦ I. 

IG0C9 

IOTA  • SETA 

IG0C9 

iF  I ISETA  - lOTAl.GE.O.S  I 

IOTA  • IOTA  ♦ 1 

leoct 

ISKP  ■ 1 IOTA  - t ) / 4 

100C9 

IR  • IOTA  • 4 • ISRP 

100C9 

ISKF  • ISKP  ♦ 1 

100C9 

SYND  • SVMOIIR) 

100C9 

SO  PRINT  U*  ISTA*  USTA*  XID*  ! SPACEID*  L>1*  tSXP  )»  SYMO 

2SN09 

lUO  CONTINUE 

2SN09 

110  CONTINUE 

2SN09 

GO  TO  StO 

2SN09 

ISO  • CONTINUE 

— ^ 

120C9 

1 • 0 

120C9 

00  21U  .IS  • 11*  12 

120E9 

ISTA  • IS  « 2 

llOff 

00  200  US  a Ul*  U2 

II0C9 

• USTA  • US  - 2 

tiott 

I 

i 

t 


1 • t ♦ 1 2SN09 

ir  ( OMCOA.CO* THETA  I.  00  TO  100  120C9 

BETA  • StCMA  • ( Xdl  • THETA  ) ♦ 1.  100C9 

IOTA  « BETA  100C9 

tr  ( (BETA  - 10TAI*GC*0«5  t IOTA  • IOTA  ♦ 1 100C9 

ISKP  • ( IOTA  * 1 I / A T0OC9 

, IR  ■ IOTA  • A • ISKP  IOOC9 

ISKP  • ISKP  4 t 100C9 

SVMO  • SVMBdRi  100C9 

IBU  PRt.«T  1»«  ISTA.  JSTA*  Xdlt  I SRACCdl*  L>1*  tSKP  )*  SVHO  12DE9 

aOO  CONTINUE  12DE9 

210  CUNTINJE  12DE9 

GO  TO  990  . 120E9 

120  PRINT  15*  II*  Jl«  Xdl  29N39 

990  RETURN  11N09 

C INO  SRtOT  S lIMVO 

BND  ieoc9 
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HflMN 


i 
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SUlMtoUfiNC  SKOT  4i  t X*  V*  hCNO  I 

C • • • • THE  LATEST  IICVl$IO*i  OATC  FOK  THIS  KOUTIKE  IS  > • XS  NOV 

OtMEt%SlUN  XtNENOI*  VTNENDt*  SFACEt  SI*  SVHBtLt 
COMMON  /MLOT/  II*  12*  Jl*  JZ 

OATA  space  / ••AH  /«  SYM8  / AH*  *AH  • *AH  • *AH  • / 
••••••••  THIS  ROUTINE  ORINTS  AND  PLUTS  THE  NENt)  VALUES  OF  X ANO  V 

NCQINNINC  WITH  THE  INITIAL  VALUE  TRANSFERRED* 

THE  MAPER  SMOULU  SC  POSITIUNEO  PRO*ERLY  ANO  ALL 
HEADINGS  PRINTCU  BEFORE  CALLING.  F.  L.  C. 

••••  INPUT  - X*  THE  FUNCTION  TO  BE  PLOTTED 

Y*  THE  FUNCTION  TO  BE  PLOTTED 
NCNO*  THE  number  OF  X OR  Y TO  BE  PLOTTED 
••••  OUTPUT*  NO  ACTUAL  VALUES  ARE  RETURNED  TO  THE 

CALLING  ROUTINE  * THE  VALUES  ARE  PRINTED 
ANO  PLOTTED  VERTICALLY. 

10  FORMAT  I SX*  I2*1X*IS*1X*E10.3*7AA*C10.S*TAA  I 
IS  format  ( SX*I2*lX*I3<U*Clw.3*2SX*£10.3  I 

IF  ( NCNO  *LC.  0 t GO  TO  SSO 

IW  > 20 
IX  > IW  / A 
tSKP  > IW  / S 
ISRP2  > ISKP 
ISKPT  • a 
iSKPF  • ISKP 
SYMO  • SYMBIli 
SYM02  • SYMBIII 
WIDTH  • IW 
OMEGA  > X(ll 
TNETA  « Kill 
02  • YIll 
^ T2  • YIll 

IF  I NCNO  *C0*  II  GO  TO  120 

00  SO  1*2*  NCNO 

IF  ( OMEGA.LT.xni  I OMEGA  • Xllt 

IF  I THETA.GT.Xtll  I THETA  > Xill 

IF  ( 02.LT*Y(1I  I 02  • Yin 

IF  ( T2«GT*Y(|}  I T2  • YIll 

SO  CONTINUE 

. IF  ( OMEGA.CO. THETA  I GO  TO  SO 

I SIGMA  • I WIDTH  • 1.  I / I OMEGA  - THETA  I 

SO  IF  ( 02*C0«T2  I GO  TO  70 

. S2  • I width  - I I / I 02  - T2  I 
70  tONTINUE 

. I • 0 

00  110  US  > Jl*  J2 
USTA  < US  - 2 
■ 00  100  IS  • 11*  12 
ISTA  • IS  > 2 
1 • I ♦ 1 

IF  I 0MCGA*C0. THETA  » GO  TO  BO 


• IS  - 2 


GO  TO 


tCTA  ■ sigma  • ( Kill  > THETA  I ♦ 1*  .. 

IOTA  ■ BETA 

IF  1 ISETA  > IOTAI*CC.O*S  I IOTA  • IOTA  ♦ 1 
ISKP  • I IOTA  * 1 I / A 
|R  • IOTA  > A • ISKP 
ISKP  • ISKP  ♦ I 


ISNOt 

s«  revised 

ISOCO 

2SN0S 

ISOCf 

lOOCS 

lOOCt 

lOOCY 

lOOCY 

lOOCY 

lAOCY 

lAOCY 

lOOCY 

lOOCY 

lOOCY 

2SN0Y 

2SN0Y 

1|N0Y 

IBNOY 

lAOCY 

lAOCY 

IINOY 

lANOY 

lAOCY 

lAOCY 

lAOCY 

lAOCY 

lOOCY 

lOOCY 

lAOCY 

lAOCY 

SSNOY 

lOOCY 

lOOCY 

lOOCY 

lAOCY 

lAOCY 

lOOCY 

2SN0Y 

lOOCY 

lAOCY 

lAOCY 

2SN0Y 

2SN0Y 

2SN0Y 

2SN0Y 

2SN0Y 

2SN0Y 

2SN0Y 

lOOCY 

lOOCY 

lOOCY 

lOOCY 

lOOCY 

ISOCY 

lOOCY 
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(tMO  • IMC9 

•tt  If  ( ot.n.Ti  t to  to  «0  l«OCt 

•CU  • U * < V(ll  • t2  » ♦ 1*  IMCt 

IOTA  • KTA  lAOCt 

ir  t (tCTA  • IOtA).6C*0.»  I IOTA  • lOTX  ♦ 1 lAOCt 

lOKAt  • I IOTA  - I I / A lAOCt 

SO  • retA  > A • ISKP2  UOCt 

tVMOa  • SVNSdAl  lAOCt 

ISMI  • ISKA}  « 1 lAOCt 

ISAAt  • U - ISA»  ♦ I 150CA 

00  OOIM  U,  t&TA*  JSTAt  XUI*  ISPACetK)  •K»1«ISK9|  • SVMO  2SN09 

\ • ISOACCaitL-ltlSAATI.  V(tl«  IS»ACCIK)«K«ltlSRAa)»  SVM02  ttNOA 

two  CONTlNUI  aSAOt 

110  CONttMUC  aSNOt 

CO  TO  too  2SK09 

IIW  ORIMT  tS»  It*  Jl,  XUtt  Yd*  a»N09 

m RfTOIIN  11N09 

C CNO  IKOT  A lOMVO 

CNO  ItOCf 
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CO.'TTRACT  OACA  23>7C>C*CC7A  UMITS  AAC  LAS  ANO  INCHrs 


601 K 

BRIDCC 

approach  slar 

IFROM  R6F  3.  SIMILAR  PROR  IN  RtP 

1) 

0 

1 3 13 

2 6 2 

1 1 

0 

0 0 

12 

16  2 

.AOOEOl  2.ACCC4C1 

2.5006-01  4.C006406 

o 

4 

O 

o 

• 

0 

0 

15 

12 

15  1 1 

0 1 

A 

0 

A 

16  1 0 

1 1 

0 

12 

12 

16  1 0 

0 0 

• 

0 

0 

12 

16 

•1.5006402 

0 

1 

12 

15 

-1.5006402 

1 

0 

11 

16 

-1.5006402 

1 

1 

11 

15 

-1.5006402 

0 

7 

12 

7 -1.7T7C40t 

1 

7 

11 

7 -1.7776408 

6 

C 

6 

16-1.777E408 

• 

6 

1 

6 

15-1.7776408 

0 

0 

12 

1 

1.0006430 

0 

7 

12 

15 

1.AA06404 

1 

7 

11 

15 

1 .6606404 

* 

0 

7 

12 

16 

1.4406404 

1 

7 

11 

16 

1.4406404 

, • 

1 

12 

16 

• 

-6.0006404 

Jl 

1 

11 

16 

-0.0006404 

1 A 

0 

0 

3 0 0 3 3 

3 

•lOOOO-lbOOO-lQOOO-10000 


2 2 0 0 >3  0 

•>000  >$000 


1 

2 

1 

15  1 

-A 

2 

1 

A 

-8 

JR5 

6-2561 

SO  Taxi WAY 

SLABS*  2 C 

-5A  BOGIES  CENTCRCO  ON 

25  PERCENT 

JOINTS 

1^ 

1 A 

A 

0 

2 1 11 

1 

1 1 

20 

28 

3.0006401  2.1506401 

2.0006- 

01  5.0006406  1.2006401 

1.2506402 

3 

8 

17 

21  1 

0 

0 

1 

0 

20 

20 

21  1 

1 

0 

0 

0 

21 

20 

21  1 

1 

0 

0 

f 

0 

5 

28  1 

0 

1 

0 

10 

0 

10 

28-2.8126408 

10 

1 

10 

27-2.8126408 

, 

0 

16 

20 

16 

-2.8126408 

1 

16 

18 

16 

. -2.8126408 

1 7 0 0 -2  1 -1  1 1 1 2 1 -1  >2  1 -2 

0>30b00>30000-30bw0-s0000>s0000->0000 


1 2^  10  10  10  20 

JM«  A*2>FT  SO  TAXiWAY  SCAOS*  2 t>7A7  BOGUS*  2>  PfACINT  JOINTS 
-1  OAOO  '2  1.  000  11 


3 

8 17  21 

1 

0 

0 1 

0 

21  20  21 

1 

1 

0 0 

f 

0 f 28 

1 

0 

1 0 
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21-  -3. blit -02 

# 

. 

9 

21  -4,40<it-»,^ 

# 

10 

21  -*-.741t'.-02 

# 

11 

Cl  '-H.4'42E-tl2 

« 

12" 

2J  -3.67a»'>o2 

• 

t 

13 

21  -2.33'5E-02 

• 

./• 

16 

21  -1.34lt-02 

/ 

# 

15 

21  •-6.37rE”03 

• 

16 

71  -c.?ilr-.«5 

• 

17 

21  -1.40JE-OH 

• 

16 

21  7.2V»t-UH 

• 

19 

21  1.04»f>')3 

• 

20 

21  -l.l6J<-:-U3 

• 

0F.FLcCT10<i>*  0ei<*£CN  ( 0«  21  ) AnO  ( 2tt  21  ) 


* 

« Y 

OcFtCCTlO^ 

0 

21 

I • |43r -0  J 

1 

21 

1 o'*'*  1'  "'J J 

2 

21 

3 

21 

-1 .OO jt-04 

6 

21 

-2.2Slt-03 

5 

21 

-•i.37eR-uJ 

6 

21 

-l.Jili.-'JC 

7 

21 

-2.3h-^-0.2 

0 

21 

-J.«?7-r-oc 

9 

21 

-•*  .4 

10 

21 

-4, 74lK-'»< 

11 

21 

-H  1 -'iC 

12 

21 

• 

13 

21  -2.3"*i»“ii2 

• 

!*• 

21 

• 

IS 

21  -s.aT^f-oJ 

• 

16 

?1  -2.2S|r.-J3 

• 

IT 

21  -1  .*>0  jt-u-* 

• 

IS 

21 

• 

11 

21  1«Uh1£*U3 

• 

20 

21  l.l*.3r:-«J 

0€fl.ttTlu atlvtcN  C 0»  0 1 

AND  ( 

0* 

2t  » 

X 

» r 

OtFLCCTJON 

0 

6.32ot-')3 

• 

1 

1.663E>0i 

• 

Q * 

2 

-1.2S7t.-03 

• 

3 

-6.7-»7t-03 

• 

4 

— <.J33»  -Uj 

• 

5 

-l.Sl7f;-U2 

• 

6 

-2,23'.t-02 

• 

7 

«3«0f»S£*'02 

• 

S 

«3«S6SEI”02 

• 

1 

— 02 

• 

10 

66,'V«7  ,-02 

11 

•02 

• 

12 

— 6.|(;  -1  — V2 

• 

13 

-3«  ' )i 

« 

-2..6»>t’'»2 

• 

li 

-2 ,73jE-02 

• 

16 

-3  •OJci-02 

• 

IT 

-.■..S66E-02 

• 

IS 

-'!•  12«E-02 

• 

11 

-02 

9'- 

20 

-6.S62£-o2 

• 

21 

-6.612E-02 

• 

22 

-3.613t-02 

• 

23 

-2.S66E-02 

• 

26 

-2.002E-02 

• 

2$ 

-1 •33jE-U2 

• 

26 

-6.el3£-03 

# 

2T 

•1.122E-03 

% 

' 20 

6.1S1E-03 

• 

, 

• 

X vrt-tfc.Nrs 

V 

OMtY*  *9eT..JCN  i §•  1% 

I A*M  ( 21* 

21  1 

X 

« 1 

20 

-6.023E-11 

m 

20 

-1  »67'»t»0< 

m 

2 

20 

-4.70«t*02 

• 

3 

20 

»o«: 

• 

6 

20 

-1 

# 

S 

20 

-2.2  «Jc*0J 

• 

i 

I 

f 

1 


I 


tot 


6 

20 

# 

7 

20 

• 

d 

2U 

* » 

9 

20  S.2i»'*k*04 

# 

10 

20 

• 

li 

20- 

# 

12 

20  3.S72t«Oi 

* 

13 

20  ••l*o/l£*OJ 

• 

l*> 

20'  -2,*‘»eK*0j 

• 

‘ 

IS 

20  -4{.2<5jfc*t»3 

# 

16 

20  -l«o'*^£«Oj 

• 

17 

20  “•»«92Jt«02 

# 

• 

19 

2<i  “4.7-l!;t»oi; 

• 

• 

19 

20 

• 

20 

20 

• 

0 

21  “I  •0«)*4€“I  1 

• 

1 

21  <-l«**ost*02 

# 

2 

21  “4«Sh6£*02 

• 

3 

21  -9,*7lt*62  • 

# 

4- 

21  -l.S-«4t.»0J 

5 

2i  -2,<12r:>'M 

• 

6 

21  -2«4l5c*o3 

• 

f 

7 

21  •1.121£*0J 

« 

i 

■ 

21  S,«2Se.*«3 

* • 

t 

♦ 

21  4.4«*‘.e*03 

# 

j 

10 

21  1.416c«9i 

# ■ ' 

11 

21  a«64't£*0J 

• 

1 

It 

21 

• 

1 

13 

21  •1.12ib'«U3 

' 

14 

21  -t.alQt^oj 

« 

IS 

21  -2.2l2t**U3 

# 

* 

1 

14  21  -l.S»44e^03 

# 

j 

r? 

21  -'».471t*02 

# 

! 

21  -•♦tS*<4e*o2 

• 

1 

If 

21  •l.45nt'>«2 

• 

» 

tt 

2}  '•a«339€>12 

• 

•> 

t 

i 

1 

t 

. X n»9t;NTi  0Ni.it*  hctwccn  c 

t*  21  i AND  ( 20*  21  1 

j 

} 

K 

1 

$ 

21  -I.OOHK'tl 

• 

.21 

• 

* 

; 

2 

21  -4.4»4£*o2 

3 

21  -9.>7K*o2 

# 

4 

,21 

• 

S 

2l' 

• 

4 

21  -2.4t<i£*o3 

• 

1 

21  -l.t2i»>93 

• 

t 

21  f.429l*f3 

• 

• 

9 

21  4.4 

• 

It 

21  1.4la4*93  “ 

• 

» 

11 

21  4,4Ha£:*'«i 

• 

It 

21  *i."2VE*‘l4 

• 

• 

13 

21  -1.121r:«»i 

• 

- 

14 

21  -2.4l4f*#i 

4 

, 

If 

21 

• 

16  il'  ♦ 

17  ?l  • 

14  21  r^,obtK*‘J‘d 

14  21  '•I.**“i4t«u2 

20  21 


i 

» 


13 

-i.b?5t*0l 

• 

U 

a -l.CJ4k««2 

« 

15 

a -6.73R€*til 

l» 

16 

6 ‘-o.Obtrl^Ol 

• 

17 

a -3.b5l£«'Jl 

• 

3 

9 -J.'j5«»K*Ul 

« 

4 

9 -®.S3-Je*0l 

• 

5 

9 -y.25Ht*«l 

• 

6 

9 -l.03ffc*t»Z 

• 

7 

9 H«21«c*0l 

• 

a 

9 i.rife*o2 

• 

V 

• 

10 

9 l.®3?£*U2 

• 

11 

9 2.l50t»OZ 

• 

a 

9 l.l37t*«2 

• 

13 

9 8.212(i*Ul 

» 

14 

9 -l.036t*02 

• 

IS 

9 -9.2S^E*Ul 

• 

' 

16 

9 -6.S39k*Ul  . 

• 

17 

9 -3.ob9£*0l 

• 

3 

10  -H.04lt*0l 

• 

10  -0.b3£t*0l 

• 

5 

10  -9.5Ufc*01 

• 

6 

10  >1.0h1£*U2 

* 

7 

10  i.Hlc:»02 

# 

8 

.10  1.779t«ai: 

4 

10  2,32efc*u2  . 

o', 

If 

10  2.216E*02 

• 

It 

10  2«32bt*02 

• 

U 

to  1.77vt.»02 

• ♦ 

13 

10  l.Hie*U2 

•f 

IV'IO  •1.04lc*tf2 

• 

L§ 

10  >9.ba3E*Ut 

• 

16 

10  >4.«32E*01 

• 

IT 

10  -it.OHlE^Ol 

• 

.11  -H.OaTE*®! 

# 

4 

11  -6,V2be*0l 

• 

n 

11  -9.7H->e*0l 

• 

6 

11  -1.075e*02 

• 

r 

11  1.67®E*0«: 

■ • 

4 

11  *^.l53t*02 

♦ 

11  4.hO->£«U2 

If 

11  3.Q-iiE*02 

It 

ll''fc,‘*09t*O2 

*' 

l« 

11  4.15jE*i»2 

• 

13. 

11  1.67'3E*02 

14 

11  -l.Otit^'lZ 

# 

15 

11  -9.7vj>E*')l 

• 

16 

11  -t>t92ac*0i 

# 

17 

11 

• 

3 

12  -3,rr9<t>)l 

• 

4 

12  «'>«77'»t*Ol 

• 

5 

12  -9.b9ct*ul 

# . 

6 

12  -1.0/it*‘)2 

# 

.? 

12  l.047f.*04f 

• 

a 

12  l.-itT-iE’iiZ 

• 

9 

12  l.7lJ'‘*)2 

• 

10 

12 

• 

It 

12  l.7li»=*0< 

• 

12 

12  J«s73'»0<' 

• • 
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13  1^  1 tOHi 

“I 

15  *';i 

lb  12  — 

IT  12 
3 13 

b IJ  -e«<»Hif*01 

5 13  ••'f.C-<2c*Ul 

6 13  -1.03St»U2 

T IJ  -i.3lit*0l 
8 13  9.2l<.c»0l 

» IJ  9.R-)/£.»Ol 

10  13  b.<)7i)t»i)l 

11  13 

12  13  v;«;l-4c*Ul 

13  13  -9.31»t*0l 
lb  13  -1.053£»02 

15  13  -9,C92t»01 

16  13  >6.b3jL*Ul 

17  13  -3.7.b<»i:»ol 

3 1%  -3.2<J>it*ul 

b 1** 

5 U -7.7ij!E*ol 

6 lb  -a.753t»01 

7 lb  >a»6bJ£*Ul 

8 lb  •V.6‘Ut'*Ul 
6 lb  -KlO'tt^Uf 

10  lb  -1.0Hb£*u2 

11  lb  >1.109C»')2 

12  lb  ■•.9.63lt*0l 

13  lb  •J*6b3£*Ul 
lb  Lb  <rb.7^)a£*ijl 

15  lb  -7,79jf4Ul 

16  lb  -S.6bbC*al 
.17  lb  -3,29>t:*ul 

3 l5,-3.1bbE*Ul 
b’  15  >S*30>i£*Ul 
$ IS  >7.I13e*Ul 

6 IS  •7.b>J7£*oi 

7 IS  -a.db-JE^ul 

• IS  -1.230^*02  ♦. 

b IS  -l»b63t>02  ♦' 

10  IS  -l.S0bi:*u2  • 

11  IS  -l»b<i'»E»u2  • 

12  IS  -1.2JC£.*U2  • 

13  IS  -■j.aSoE^’Jl 
lb  IS  -7,bb7t»'il 

15  IS  -7«l|3c*ul 

16  IS  -5.*J34E»01 

17  IS  -3» I 
3 lt»  "J«3SJt.'«vi 
b IS  «S.63';*t>«i 

5 lb  -/.77vit,*\»l 

6 16  -9.<)31t*«l  . 

7 lb  -V,V7»t.‘»ul 

8 lb  >*1  • l'i<:£»il2 
8 lb  •! •♦112 

10  lb  >l.l-iSr.*M2 

11  lb  -1.2**i*g2 

12  lb  '’l*lS*f*02 
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W. 
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# 

15 
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lb 

lb 

• 

17 
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# 

3 
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« 
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17 

17  -:i.b3;t«ttl 

« 

3 

H - ji 

# 
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• 

5 
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• 

6 
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♦ 

7 
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• 

10 
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• 

11 
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• 

12 
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• 

13 
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• 
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• 

n 
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• 

u 
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• 

17 
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# 
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0 
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10 
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# 

12 
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4 
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• 

s 
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• 
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20  2.39U‘*92 

• 

10 

?o  2. i««it:-«^u^ 

. • •• 

11 

20  2..19tt*U2 

• 

12 

20  ^•''■*lt»92 

# 

106 


13 

2U  1.07tc.*u^ 

# 

il- 

20 “1  .'i2<s,»Ui 

• 

ls 

20  -■i.S2tff‘Ul 

• 

16 

20  '>6.ni^f.«Ul 

• 

17 

20  •H.O'j-t’Ol 

• 

3 

2 1 ••3«'»7rK  »oi 

• 

h 

21  “6.6 7 -jf.  *01 

• 

5 

21  ■•s,3ytfc*ol 

♦ 

6 

21  -1  ♦Oi»le.«u2 

♦ 

7 

21  1.7-<''>»;*02 

• 
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V 

21  »u<r 

10 
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# 

11 

21 

12 
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PR0G«»4;:  sLAhSOt  CK-'L  bt-t-ClAL  UECR-Awt-jjP.KLE 
SM-ii'i.E  fr'd-'H  -S  lo  ot  ^0'4Sr-<*fE  PROGRAM  0*,£ 

covr^Acr  oala  ^3-V'j-c-oo/6 


revision  date  07  SEP  70 

- Zm  SEPT  7u  - JJPffLEtJLR 
UNITS  ARE  LrJS  AnO  INCHES 


PROR 

JR6  A-25FT  SO  TAXIWAY  SLABS*  2 B-747  BOGIES.  25  PERCENT  JOINTS 


TABLE  1 . CONTROL  bfLTA, 

« 


HULTIPLt  LtiAO  option  IIP  BLANK  OR  2ENO*  PROB  IS  INDEPENDENT  — -1 

IE  *1.  PArENf  EUR  neat  PROb  — If  -1.  AN  OEESPnING  PROBI 


table  number 

2 3 4 S 

NO-«  CA-IOS  I'PUT  TmIS  PROhLEa  0 4 0 0 

COMPUTE  OP  1 10N«L  PAVEHtNT  STIEENESS  CONSTANTS  a«TES)  OX*l)Y 

0 

OPTION  TO  sOPRtSS  OETaILEO  OUTPUT  (IsYt'S) 

OPTION  lu  PRi-iI  PblN  STRESS  instead  OF  MUM  UbTES) 


6 

2 

C 

0 


7 

1 

S 

0 

1 

I 


TABLE  2.  CONSTANl'S-UMITS  MUST  BE  CONSISTENT. 

, NUi  INCPt'iP'iTS  IM  A UIpECTIOn 
■ NU*  lNC«c Stars  Ifl  r OIrECTION 
INCR  LE'nGIH  I'l  A OlRtCTION 
INCR  length  I-m  r OIRECIION 
POISSO  iS  ratio 
MOOOtUS  0'  tLAstlClTY 
SLAB  T**lCsatsS 
SUdORA.It  sviLOLoS 


20 

2S 

3.000E*01 

2.1S0£*0I 

2.000E-OI 

S.0U0t*06 

1.200E*0I 

1.2S0E*02 


TABU  3.  SPEOPIEO  AREAS  FOR  SELECTED  PLOTTED  OUTPUT. 


from  THrO 
3 B 17  «l 

0 21  Zj  cl 

9 0 y <0 

12  0 U <rt 


PLOT 

OEFL 

1 

I 

I 

1 


U«YES» 

A-MOMENT, 

ft 

I 

ft 

ft 


Y-MOmENT  . PRIN  STRESS 


0 I 

ft  ft 

1 0 

I ft 


TABLE  4.  LOAD  DATA-REPLACES  LOAD  IN  PREVIOUS  PROBLEM. 

ALL  allFf.iESs  ItR-S  ARt  RETAINED  PROH  PARENT- PM08L£<T  JRS 


'FROM  Tp-<l 


e 


M)\t 
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table  5.  AXIAL  THRUST  DATA. 
fl'OS  ImJc 


table  S.  SPECiAL  LOAD  PATTERNS. 


pattern  no-  uf 

NU<-'  LOtJb 


i^tFEMtwCt 


HATTEHw  COOPOI.xjAies  AND  CONCE.NTWATtD  LOADS 
'.rr  XT  X 


" “ ® 3 10  2 0 it 

>HlbOO  -Ali,oo  -AISOO  -<HS00  -ilSOO 


0 0 

-suou 


-3  0 

-5000 


TABLE  7.  PLACEMENTS  OF  SPECIAL  LOAD  PATTERNS. 


PATTEwn  no-  ijF 
NOH  PLACf.E„ 


2 3 

-A1500 


PLACt-EWS  A y OP  REFENEnCE  load  (SLAB  COOHOInaTES) 

* P ^ ^ Y X Y j 

^ 7 12  12  ia 


SUM  nr  “‘’Pl.lEO  10  ThE  SLAB 
SUM  OF  ^MtfcL  LOAOS  APPLltU  10  iHt  SIAM 
NUM  or  .«aL  LJAD,  plaCco  uursiot  Iu3 
SUM  OF  ,«tLL  LOAJS  PLALtu  uoTSlDt  SLAtf 


12 

-A,9dOE*OS 

0 

0. 
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Ct.'L  5>r'cCl«L  UEC(<-A.<£-JJ^.f LE 
Si  n^Lh  ■<>  Tu  Ot'jy^fSrWATt  f'KOGrtAM  USE 

CO-lI<tACI  <iJ- 7O-C-O076 


t?EVISIO<-4  DATE  07  SEP  7U 

- SEPT  70  - JJPtfLEtJLK 
' units  ape  lbs  and  inches 


PPO:»  ('■)  ID)  ■ - 

JK6  A-25FT  SQ  TAXI  WAY  SLABSt  2 3-747  BOGIES*  25  PERCENT  JOINTS 


TABLE  8.  BESULTS-USING  STIFFNESS  DATA  FROM  PARENT  PROBLEM  JR5. 

A .y,t..r  AVI)  A i-.isrifn>  .iu-!Eni  aci  In  ime  a oihection  (about  y axis». 

Y UISIHO  = -X  (-riSTlNti  MOAt'wl*  COUNTEkCLOCKwISE  BETA  ANGLES 

AJt  POoiriVt  fHU-(  X AXIS  r«)  THE  UlPtCTlON  OF  LAkGESI  principal  stress 


A 

LARGEST 

beta 

X 

Y 

twisting 

SUPPORT 

PRINCIPAL 

X TO 

X . Y 

oefl* 

MOMENT 

MOMENT 

MOMENT 

reaction 

1 

STRESS 

t 

LARGEST 

•••  the  OETaILEu  OuTPUI  mas  been  UELETEO  at  THE  option  in  table  1. 

A SINGLE  SET  OF  VALUES  IS  PRINTED  AT  OR  NEAi<  THE  SLAB  CENTER  FOR  REFERENCE 


10  14  -S.73jt-U2  r2»4<jOE*ul 


2.17«E*03  -?.16SE*03  4,622E»03 


I.461E*02  -Sa.5 


Y ■ OEFL 


A 

/ 

Y 

X 

TWISTING 

SUPPORT 

largest 

PRINCIPAL 

beta 

X TO 

MOMENT 

MOMENT 

MOMENT 

REACTION 

STRESS 

LARGEST 

statics  thECx. 


SUHRATION  OF  REACTIONS  « 4.9«OE*05 


* I 


1 


HQ 


PfiOGPui-"  3(.'*rJJt  Ct'L  i**£,Cl4l. 

SA-1f>Lt  “f' .'iliLtni  TU  USE 

COWTrtACI  jACA  £3- ^U-C“uII7«> 


REVIS!a>^  DATE  tf7  SEP  70 

- ih  StPT  70  “ JJP.PLEtAP 
UNITS  AAE  lbs  ANU  INCHES 


PPOtt  <C.)  TO) 
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IT  12  «9.6V<E*01 

3 13  •l*2Tl£*02 

♦ 13  *1.03le*u2 
5 13  -l»*^ne*02 

• I3'''l,5'i»£*a2 
T 13  O.l7i(.*02 
•'  13  5a31<r*u2 


>-12 
TO  12 


0 13 
5 13 
0 13- 

T 13 


♦ 

10 

11 

■ 1>~  o,oi<»t»o2 

13  1.00j£»02 

13  *'l«3&ii*v2 

12 

13  **i«2‘)9t*02 

• 

13 

13  ••l«l2/t;»iu 

• 

10 

13  **l««^c£*«i2 

• 

15 

13  •lt21er>02 

10 

13  **l«00<>t->*i2 

# 

TT 

13  rT>29lt:*at 

• 

3 

lo  ■>l*207t'*U2 

• 

0 

lo  -l«7-77t‘*02 

• 

5 

lo  •l»S2»!f*02 

• 

0 

T 

lo  1.3 17  V:  *02 
lo  «.20<«**y2 

♦ 

I|  1*» 

9 19 

10  |4  l»**il£*U^ 

11  l9 

12  19  ••l»oTJt*02 

13  l9  -».93il£*ul 

U U -1.1Ji2E*02 
15  19  -1.9Ht*02 

19  19 

IT  19  -1.01't*(/2 
3 15  -l.2ilc*i)2 

9  IS  -l«7lut*oi 

5 IS  -1.7Jtt*U«{ 

6 15  2.23-3E*u2 

T 15  S.3«IC*()< 
a 15  6«17et*U2 
9 IS  6.l9tC«U2 

10  15  3.01'jt.*a2 

11  IS  •>2.7<>cE»U2 

12  15  -l.V7Se*02 

13  15  -I.OOsC*^^ 

19  15  -1.12«.t*02 

IS  15  -USIOE^O^ 
15  15  -l.S-'jE*l)2 
17  IS  -1.2ShE*02 
. 3 15  -1.17.;t*U2 

9 15  •1.557i»a2 

5 15  •l*595t*02 

5 15  1*592C*02 

7 15  3.055E*02 

• 15  3.953E«V2 

, 9-^5  J.7?3L»t)2 
. 19  15  2.92bF.*il2 

11  15  -2.95««K*02 

12  15  2.305<*tt2 

13  15  U7-i9r*02 

19  15  1.65be*U2 

15  15  -l,392t*U2 

15  15  'l.^STE^Oa 

. 17  15  -l.350E*02 

3 17  >l.«55e*42 

9 17  "l»59v't»l)2 

5 .17  -l.«9»E*02 
5 17-1.710E*52 

7 17  1.550£*02 

5-  17  1.71»E*U2 

9 17  2.00«F*«2 
If  17  -2.7j7E*u2’ 

11  17  3.0iaE*u2 

12  17  9.02;t»U2 

13  17 

19  17  3.21'it.*()2 

15  17  2.0Ur.«*i< 

15  17  -l.5»>.f.*i;2 

I'T  17  -l.**2lE*«< 
3 19 

9 19  «1»353€*02 

5 H -l.72sh*«2 
5 It) 

7 19  -l.5jJ-:«ii2 


I/; 

- 


i 

i 

\ 

K 


i 


IN  -i.S)jL«0^ 

• 

9 

le  ~1  .r  1 

• 

10 

Is  -X.l-inc^U^ 

• 

11 

H a.n-iE*').: 

• 

la 

18 

13 

18 

I** 

Is  6.UVt.«l'2 

♦ * • 

IS 

la  2«THjt*<J2 

# 

16 

ll  -1.*i3/c.*02 

• 

IT 

18  -l.<*33i»02 

• 

3 

U -7.S>U»ol 

• 

N 

W -1.1-Vi *0^ 

• 

5 

i->  -i.si-;l«o2 

# 

6 

19  •‘i 

« 

7 

19  -i;Tl»r>02 

• 

8 

19  -1.55-.c*02 

« 

9 

19  -l,56<t»02 

• 

10 

19  -i«967t*u2 

• 

11 

19  l.9a0r>v>2 

• 

1? 

19  9.3H3t*02 

# 

13 

H S,37-P*02 

' 

19 

19 

# 

IS 

19  2.l9St*02 

# 

16 

19  -l,«9-it>02 

• 

17 

li*  -1.90aE«02 

• 

3 

20  -5.91cf*0l 

• 

6 

20  -9,l43t*01 

5 

20  -1.2H2t*02 

• 

6 

20 .-l^NoaE^O^ 

• 

7 

20  -l.S.iie«0<' 

• 

• 

20  -1,3S=E*>J2 

• 

/ 

1 

-1.191t*U^ 

• 

j 

10 

20  -9.0lSF.*ol 

• 

/ 

11 

20  1.692r1»02 

/ 

■ 12 

20  9.3S<E*02 

# 

13 

20  5.*.?SE««2 

19 

2tf  ,‘*.2«^f*02 

• 

IS 

2U  2.173P.*02 

• 

16 

20  -l.2Vie*02 

« 

17 

20  -1.39£E»u2 

• 

' 

3 

21  -9,333l>Ul 

9 

21  -6.b95t«0l 

# 

5 

21  -9.6/2f*01 

# 

6 

21 'il.21lL*U2 

• 

7 

21  -l.3lfit*02 

• 

8 

21  -1.222c*02 

« 

9 

21  -l.0-<jt«02 

• 

10 

21  ”rt«203£*Ui 

# 

11 

21  2.1«*iE«o2 

• 

12 

21  b»''i-<r'« '•? 

13 

21  9.vj«i:*02 

• 

19 

21  6.3-5E*'J2 

IS 

21  2,»>-:)E*02 

• 

16 

21  -l.37jt.«02- 

• 

17 

21  -l.3ae£*02 

• 

flMC  FQM  This  PtTUHtCH 


fLAPSU)  THE  « 


0 HINUTCS  4.72»  SECONOS 

%% 

Q MINUTES  seconds 


♦ 
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Cf-'L  Sf'rLi*!.  l)tC^-Aw£- LE  hEVISION  0*Tt  OT  Stf*  TO 
b*-.Vi.h  *1  f ut*<0"»STw»it  i^mOGkAM  (JS£  “ 2a  SEkT  7o  - JJP*rLE«JLM 

CiiMl--Ct  j>*Ca  <;j-/u-C-«i0  7t>  UNitS  A«E  LBS  AND  INCHtb 


SACI  A » d hEO 


TAMJE 1.  CONTROL  DATA. 


HULIIPLt  LOAD  uhTIo.4  tlf  ttLArtK  Oa  2EW0»  PXOB  IS  INDEPENDEMT  — -0 

If  4l.  PAAENf  FO-i  NEXT  PilOB  — IF  -I*  AN  OFFSPHING  PftOB) 


NUA  CA^OS  npor  TnU  PROBLEM 
CONP07E  O^riUNAL  PAwEMt.'iT  STIFFNESS  CONSTANTS  (l*YtS> 
option  lu  SOPAtaS  OETaILEO  OUTPUT  Il»ttS> 

OPTION  to  pAioT  paIn  itotss  Instead  of  mom  ii«yes> 


TABLE  NUMBEP 
2 3 4 5 6 

1 0 6 0 0 


OXtOY 

0 


c 

0 


7 

0 

s 

0 

0 

0 


TABLE  2.  OONBTANTS-UNITS  MUST  BE  CONSISTENT. 


NU4  INCACmEhTS  in  A OIMECTION 
NU4  INC^EntoTS  M Y OIMECTION 
INCM  length  in  X UIMECTION 
INCP  LEnGIh  in  Y UIPCCTION 
PulSbO.NS  MAT  10 
• MOOUcUS  OF  ELASTICITY 
' SLAn  TfilCX'iESi 
SUOOMAOt  SOOULUS 


8 

6 

6.000E*00 

6.000E*00 

2.500E>01 

0. 

9.7SBE-01 

0. 


TABLE  3.  SPEOFtEO  ANEAS  FOB  SELECTED  PLOTTED  OUTPUT.  

PLOT  (i>ves> 

FROM^.  ThMU  OEFL  X-HOMEKT  Y-HOHENT  paiN  MOMENT 


NONE 


TABU  4.  STIFFNESS  AND  LOAD  DATA. 

FROM  rNPO  OX  OY  0 S 

• 0 4 

■t  0 I 

• 1 tf 

1*  1 7 

4 6 N 

1 1 • 


b OoZSOE^OS  6.2S0E*V5  >0.  1 

9 6.2iOE*Oa  6«2S6E*0S  -0»  -0 

7 o,2Si)e*U3  6.2S0E*05  -0.  -0 

7 6*2>i)E*vS  6»250E*0p  -1 

•»  •(>.  -O.  •l.tOOEAfS  -9 

6 -0«  -0*  •§•  «• 


.0&0^*20  >0. 

-0. 

•0. 

••0E*20  -0. 

-s* 

I.S7SE«96 
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^ tables.  AXIAL  THRUST  DATA. 

FHi)'l  Th-o 

» 

NOHb*  « 

TABLE  6.  SPECIAL  LOAD  PATTERNS. 


PX  tx 


PHOb^A.i  sL4-30t  Ct-t  vrtCl-L  OtCK-HW£-JJ*».»^LE  HCVISlOU  D*l£  OT  SEP  70 

Sa  Vut  Tu  utMyibfxAlt  P><OG*<AH  USE  - 2**  SEPT  70  - JJPfFLE*JtP 

COxTxACT  JACa  23-/U-C-00/O  ‘ UNITS  AHE  LHS  AND  INCHES 


P90d  (COf lOl 
SACl  4 A 3 A£(i 


TABLE  8.  RESULTS. 

• A v.(i.it...a  anO  a TxISTInO  H0m£NT  ACt  IN  TmE  X DIRECTION  (ABOUT  T AXIS)* 
t '•‘O' It.'*  I = -A  I.iISTInO  mo-jEnI.  COUNTEhCLOCKkISE  beta  angles 

.APE  Posnive  F«OA  A AXIS  to  the  direction  of  largest  principal  homent 


X 

LARGEST 

seta 

X 

T 

TUISTING 

SUPPORT  PRINCIPAL 

X TO 

X • 

T OEFL 

MOMENT 

moment 

MOMENT 

REACTION  moment  LARGEST 

0 

0 1.267E-16 

7,lvue-<fb 

-1.23JE-2S 

-1.S3AE*03 

-I.267E*0A  -I.586E*03 

6S.0 

I 

0 -2.273t-(7 

-3.3y0t-lc 

-2.S'7oe-il 

-3.U13E*03 

2.275E*03  -3.0I3E*03 

AS.O 

2 

0 -A.77'.t-17 

-1  .SAHt-l 1 

-6.ioaE-u 

-2.AV0E*03 

A.779£*03  -2.690E»03 

65.0 

3 

0 -7.33)6-17 

-j.0‘*7e-u 

-1.233E-10 

-l.A(iOE«03 

7.3S0£*03  -l.A80E*0J 

65.0 

A 

0 -'i.38-tt-l7 

-2.V/'»C“(1 

-1.231E-10 

-6.129E-10 

8.866E*03  -6*912E-10 

67.2 

S 

0 -7.35)t-lV 

-3.0a  7t-i 1 

-1.23JE-I0 

I.Ad0E«03 

7,350E*03  -I.680E*03 

-65.0 

8 

0 -A, 77^6-17 

-1 .3aaE-1 1 

-6.1b8E-ll 

2.AV0E«03 

6.r79E*03  -2.A90E*03 

-65.0 

7 

0 -2.27SFM7 

-3.3U0£-12 

-2.976£-lI 

3.U13£*03 

2.27bE*03  -3.0I3E*0J 

-65.0 

8 

0 1.247E-10 

7.SJJE-CO 

-1.0V6E-25 

l-.SMAE^Oa 

-I.267E*06  -I,586E*03 

-65,0 

r* 

I -2.2736-17 

-6.0b0E-U 

-i.iio£-n 

-3.U13E*03 

2.27SE*03  -3.0I3t**03 

65.0 

1 

1 -l.2I6F.-uI 

1.0JAE«03 

1 .0Ja£*03 

-5.773£*03 

0.  6*807E*0J 

-65.0 

2 

I -2.3IA6-II1 

2.2'»7£*03 

».%3VE*03 

-4.bd2E*03 

0.  6 >S6E*03 

-63. 7 

3 

I -3.12'.t-ul 

3. 7boE»0J 

2.20V£*03 

-3.016£*03 

0.  6..42E«0J 

-37.8 

A ' 

I -3.‘.Slt")l 

A.7ehE*u3 

2.lbS£*03 

-1.203E-09 

0,  •A.720E*O3 

-.0 

5 

1 -3.I2'»t-ol 

3.7ot)E*OJ 

2.209E+03 

3.Ult>E*03 

U.  6.I02E*03 

37.8 

6 

1 -2.3iHe-ai 

2.2u7E«u3 

1.8JVE*03 

A.rtB2E*03 

0.  6.9S6E«0J 

63,7 

T 

1 -I.2l6c-Ul 

r.ojAt*uj 

1.0JhE*03 

5.773E»03 

0.  6.807£*U3 

65.0 

8 

1 -2.275F.-17 

a.7jaE-1 I 

1 .SdVE- 1 1 

3.0I3E*U3 

2.27SE*03  3.013E*03 

65.0 

0 

2 -A.775t-l7 

-6. loUE-ll 

-l.SAAE-ll 

-2.‘*90E»03 

A.779E*03  -2.490E»03 

65.0 

1 

2-  ’•d.3  l-«F-o  1 

1 .FiJ>,t*03 

2.2V7£*03 

-♦.d«2E*03 

0.  6.9S6E*03 

-66.3 

2 

2 

A.23lt»03 

A.2blE*03 

-4.A49E*03 

0.  S.700C*03 

-65.0 

3 

2 -6. 0636-01 

7 . AooE  *03 

5.3d2e*03 

-3.155E*03 

0.  9.7S9E*03 

-35.8 

A 

2 -6.7Si't-'il 

1 .OtfoE^uA 

S.BmTE^OS 

-1.203E-09 

U.  1.0ZBE*OA 

-.0 

5 

2 -6.0«>jF-iil 

7 . Ajot  *0 J 

3.3d2E*03 

3.15bE*03 

0.  9.7S9E*03 

35.8 

6 

2 -A.A.TJt-'jl 

-.23l6*uJ 

a.23I£*03 

6.AA9E*03 

0.  8.700E*03 

65.0 

7 

2 -^.il^F-i'l 

1 

2.2y7£*03 

A.(582E*03 

0.  6.956E*03 

66.3 

8 

2 -A.77'ir -1 7 

J.Ocjc-l 1 

7>,tooot-12 

2.h90£*03 

6.779E*03  2.A90E«03 

65.0 

0 

1 -7.33)t-l7 

-1  23JE-IU 

-3.0A7e-ll 

-1.A80E*03 

7.350E*0:i,  -I.ABC1*03 

65.0 

I 

3 -3.1<'tt-0l 

2.?OVL*u3 

3.7t)*jE*«3 

-3.Ul6£*03' 

~0;  6.I02£*0J 

-52.2 

2 

3 -6.'?63''.-iM 

s.33<:t*uJ 

7.Aa*»E*03 

-3.1SSE*03 

0.'  9.7S9E*03 

-56.2 

3 

3 -'i.AA'.r -ol 

1 .OOfJ*.  ♦IIA 

1 .OooC*OA 

-3.016£*03 

0.  l.370E*0A 

-65.0 

A 

3 -O.**?’*;- Jl 

i . A 1 uE ‘OA 

l.l 77E*«A 

-4,7b9£-10 

0.  I.ei0E«06 

-.0 

5 

3 

1 .obhE*!;** 

1 .0«J»>t*06 

J.01bt*03 

0.  1.370C*0A 

65.0 

6 

3 -6  , Oh  >r  -U  1 

3.3o.tf' *03 

7,‘iOh£*0) 

3.1S5E*03 

0.  9.7S9E*0J 

56.2 

7 

3 -3.1A,-  -.;1 

/ *<1 J 

J«7tjc>6*u3 

3.0l6e*03 

0.  6.I02E*03 

52.2 

8 

3 -7.3-'  /f  -1  7 

■j.l  776-11 

l.S7yE-U 

i.'tboE^ba 

7.350e*03  I.A60E«0J 

65.0 

124 


0 

4 / 

2 ibCrtK-l J 

1.0blE>12  >2.J13E>n 

a.644£*03 

I 

4 -j.4Sl‘--ol 

2.  l'>b' .yj 

4.72OE.0J  -l.JH.-iE>l» 

0. 

2 

4 -o.  7ti- -1)  1 

■j.  /t»o3 

I.U<Mt.*04  -5.U8VC-1U 

0. 

3 

4 -'i.6^3r-yl 

1.17  /r  .y*. 

l.ol0L*04  -V.2S2E-I0 

0. 

4 

4 -I.iJ“t*UU 

J.034t*U4 

3.0bb£*04  -4.626E-10 

0. 

5 

4 -9.62j»=-'.tl 

1.177c»04 

1.01UE.U4 

2.J1JC-10 

0. 

6 

4 -6.7^'it-Ol 

b.34/E.u3 

1.02bE*U4 

6./ObE-10 

0. 

7 

4 >3.4bl6-<Jl 

£•  lo^E*03 

4.72o£»b3 

7.064C-IO 

0 . 

8 

4 -a.8f)4E-l/ 

2.47ot-10 

b.273£>ll 

2.b44E>10 

8.864E.03 

0 

5 -7.35'J'-I7 

V.lalt-14 

3.672E-13 

I.480E*03 

7.3S0E*03 

1 

5 

a.^yyr.yj 

J.7bbb*u3 

3.ylfc£*03 

0. 

2 

5 “6.0r»JK“i)I 

:3.3a2t*y3 

7.4ObE*03 

J.155£»03 

0. 

3 

5 -3 .43-». -!)  1 

1 .0o?r>U** 

l.0bdE»04 

3.01bE*0J 

0. 

4 

5 -■».c?3r-0l 

1 .biuL.y*. 

l.i77E»04 

3./0i£>t0 

0. 

5 

5 -d.4S'»f.-Ul 

l.Ob0E«U4 

l.UbrtE*U4  >3.ylb£»03 

0. 

6 

5 -O.OSjE-Ol 

5.3e2c*y3 

7.40b£»03  >3.IS5E*03 

0. 

7 

5 -3.l2'.t>0l 

2.2U'4E«U3 

3.76b£*03  -3.Ui6E«U3 

0. 

a 

5 -r.35*)e-I7 

-i.844£-10 

>4.5a9E-ll  >i.4B0E»03 

7.350E»03 

0 

6 -4.77st-l7 

-0. lod£-I 1 

>1.544£>n 

2.4V0E*03 

4,77VE*03 

i 

t>  “2*3l-»t“yl 

1 .nji(t«u3 

2.2V/t»03 

4.d82£»03 

0. 

2 

6 -4.4J3t'"0l 

4.2bl£*y3 

4.2’»iE»ll>3 

4.449E*03 

0. 

3 

6 -6.043t-0l 

7 «4ebE  *03 

5..a2E»03 

3.I55E»03 

0. 

4 

6 >6.750t>0l 

1.028E*U4 

5.347E*U3 

l.Jb5E>09. 

0. 

5 

6 >6.04je-Ul 

V.4boE«03 

5.3a2E«03  >3.IS5£*03 

0. 

6 

6 -4.43jc.*‘Jl 

4,2sli;*y3 

4.25l£*03  >4.449£*03 

0. 

7 

. 6 -2.J14C-01 

l.OiVE.yb 

2.2V7£*03  >4.ay2E»03 

0. 

0 

6 -A.TTVE-l/ 

>1.2j4E>10 

-3.0O6C-11  -2.490E403 

4.77VE»03 

0 

7 >2,275fc-iy 

-b.ObyE-tl 

-I.liO£>ll 

3.V13£»03 

2.275E»03 

I 

7,-l,2lft£-Ul 

1 .0J4C*U3 

1.034E.03 

5.773E*63 

0. 

2 

'T -2.3l4E-yi 

2.2y7£*03 

1.03VC»03 

4.802E*03 

0. 

3 

7 -J,l2-.e-yl 

3.7bbt*0J 

2.2UVE.03 

3.U16E*03 

0. 

* 

7 -3.45lt-yl 

4,728£*y3 

2.I65E»03 

1.711E-09 

0. 

5 

7 -.>.l2^t-Ul 

J.7bot*y3 

2.2UVE»V3  >3.UlbE*03 

0. 

6 

7 >2.3l4r;-ir 

■2.2v7t»0J 

1.03Ve*03  >4.»H2£.03 

0. 

7 

7 >1.2lSt-ol 

1 .0  j‘<E.y3 

I.034£*U3  >5.773£»03 

0. 

8 

7 -2.27'nfc-l7 

-7.6y2£-U 

>1.4'#bE>ll  >3.0l3E*03 

2.27SE.03 

0 

• 8 "l.2t7e-lb 

-4, I0v£-2o 

-2.414£-2S 

1.584C*03 

-1.267E»04 

1 

8 -2.275t-17 

-J,yj7e-1 1 

>l.J/7£-lO 

3.013£»03 

2.275E*03 

2 

8 -4.77'it-l7 

-4,4<:Ht-i  I 

>I.ObOe>IO 

2.*V0E*03 

4.77VE*03 

3 

8 -7.350r”I7 

-3.047c-1 1 

-1.2J3E-I0 

1.4iioe*03 

7.3S0E»03 

4 

8 -8,644t-l7 

-2.V7-»i>li 

-1.2JIE-I0 

8.V05C-10 

8.064£*03 

5 

8 -7.35f.E-l7 

-I  .SosE-ll 

-b.I5V£>U  >1.4«0E*03 

7.35U£*03 

6 

8 -..77.t-I7 

>2.347c“14 

>S.«b7£-l5  >2.4V0E*O3 

4.77yE*03 

7 

8 >2.27ot-l7 

l.?yJt-U 

3.lV2b>ll  >3.013E»03 

2.275£*03 

8 

8 l.267f.-I& 

A 

-6.048E-26  >l.584E*03 

X 

r tmistinc 

-l.267E*04 

suppohY, 

X 

. Y OEfL 

hO.hEnT 

moment 

moment 

react  ION 

«.72<tL*UJ  -VO.O 
1.024Ct0<»  >90.0 

l.tiio{;*o<»  'VO.O 

3.0S6C»0<> 

i.aiee«04»  vo.o 

1.02llE*0<»  90.0 

4.728E*03  90.0 

9.2S4E-10  35.0 

i.«»«OE*OJ  ^5.0 
6.102E«0J  52.2 

9.759E*0J  54.2 

1.370E«0h  45.0 

1.810E«04  .0 

I. 370E*04  >45.0 
9.759E*0*i  >54.2 
«*102E«0J  >52.2 

>1.480E*03  , 45.0 

>2.490E*0J  >45.0 
6.956E»03  46.3 
•.700E«03  45.0 
9*7S9£^03  35.0 
i*020E*04  .0 

9*759E*03  -35.0 
0.700E»03  >45.0 
6.956E^03  -46.3 
-2.490E«03  45.0 

-3.0|3£*OJ  -45.0 
6.007E«03  45.0 

6.9S6E*03  43.7 

6.102E*03  37.8 

4.720E»03  .0 

6.I02E»03  -37.0 
6.9S6E«03  -43.7 
6.007E*OJ  -45.0 
>3.013.<’*03  45.0 

-1.504E*03  -45.0 
-J.013E*03  >45.0 
-2.490E*03  >45.0 
>i.400E*03  -45.0 
-9.602E-10  >46.5 
>l.400C*03  45.0 

>2.490£*03  45.0 

J. 01je*03  >45.0 

-l.S0<*E*03  45.0 


LAOGEST  OETA 
PRINCIPAL  X TO 
moment  LAKbEST 


I 
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